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A system at filling factor 2/3 could be a candidate for a quantum Hall ferromagnet at integer filling factor 
of composite fermions. Using exact diagonalization with electrons on a torus we study the transition from 
the singlet ground state to the polarized ground state at this filling and look for signatures of quantum Hall 
ferromagnetism. Differences between the fractional and corresponding integer systems are emphasised. 
Most interestingly, we find around the transition a low excited half-polarized state which might become the 
ground state in the thermodynamical limit. We study its structure and compare it to the singlet and polarized 
ground states. A new interpretation of the singlet state is suggested and comparison of the filling factors 
2/3 and 2/5 is presented. Adding magnetic inhomogeneities into the system we investigate the stability of 
all the three involved states and the tendency to build up domains like in conventional ferromagnets. 
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K. Vybomy: Spin in fractional quantum Hall systems 



Similar to superconductivity, the fractional quantum Hall effect [14, 79] is a unique field, where correla- 
tions between electrons give rise to macroscopically well observable ground states which we would not 
expect on the level of a Hartree-Fock approximation. The correlations are introduced by interelectronic 
interaction and, contrary to atomic physics for instance, the quantization of single-electron energy levels is 
a consequence of the strong magnetic field (Landau levels) and of the suppressed motion in the direction of 
the field (quasi two-dimensional systems). The latter phenomenon leads to another unusual feature of the 
fractional quantum Hall systems: the many-electron states in a non-interacting system are highly (macro- 
scopically) degenerate, since all electrons within one Landau level have the same energy. In particular, for 
filling factors below one, where it is useful to be restricted to the lowest Landau level, all many-electron 
states have the same energy. Now, the effect of interactions between electrons cannot be investigated by 
perturbation theory, as there is no single ground state to start with or, in other words, there is no small 
parameter in which we could expand the perturbation series. Since energy spacing between the many-body 
states is zero, the interaction is never a small perturbation, regardless of how weak it is. This fact ren- 
ders the fractional quantum Hall systems unique from the theoretical point of view and makes completely 
novel types of quantum-mechanical ground states possible. The best known of these are the incompressible 
quantum liquids. 

Quantum Hall ferromagnetism was one of the companions of the integer quantum Hall effect (Subsect. 
1.6). The observed long-range spin order can be explained by exchange energy gain in the ferromagnetic 
state and hence Hartree-Fock models are sufficient to describe the ongoing physics. However, new ex- 
perimental publications appeared in late nineties. Phenomena reminiscent of ferromagnetism have also 
been observed in Ihe fractional quantum Hall regime, being most noticeable at filling factors 2/3 and 2/5 
([45, 15]). In this situation, the Hartree-Fock approximation is no longer acceptable, the spin-ordered states 
are highly correlated. This area is not very well explored. Instead of a lattice of spins which are all point- 
ing in the same direction, here, we are dealing with itinerant electrons which are either in a fully polarized 
or in a spin singlet state (Subsect. 2.1). Although both states are incompressible, their structure is quite 
different [13]. 

How far can we extend the analogy between an Ising spin-lattice ferromagnet and fractional quantum Hall 
systems where two ground states with different spin order compete with each other? This was the central 
question of this thesis. There are several fundamental differences between these two systems. The latter 
one is itinerant and the liquid-like ground state is stable only owing to correlations while, in a spin-lattice, 
the electrons are spatially fixed and the ferromagnetism occurs also in classical systems with suitable 
site-to-site coupling. By observing e.g. hysteresis in magnetotransport, experimentators have provided a 
lot of evidence that the two phenomena are indeed very closely related [71, 70, 43], on the other hand, 
observations without an analogy to usual Ising systems have been reported too [46]. 
Before we start the theoretical introduction, let us summarize the basic experimental facts. At filling 
factors 2/3 and 2/5 two different ground states may appear. Depending on the ratio between the Zee- 
man to Coulomb energy, EzjEc, it is the fully spin polarized {Ez/Ec —* oo) or the spin-singlet one 
(Ez/Ec 0) [46]. This transition can be accomplished (a) by varying the electron density at a fixed 
filling factor [36], (b) by tilting the magnetic field [16, 21] or (c) by applying hydrostatic pressure which 
modifies the bulk g-factor [50]. 

When the two ground states are brought to degeneracy, transport experiments show hysteresis, time- 
dependent resistance with Barkhausen jumps [11] (see Refs. above) and huge longitudinal magnetoresis- 
tance [45] which is related to the spin polarization of the ion lattice of the hosting GaAs (NMR experiments 
[44, 19]). These could be related to formation of spatial domains of the two ferromagnetic ground states, 
even though surface acoustic wave experiments could not confirm this [20]. On the other hand, optical 
experiments [46, 24], suggest that a half -polarized ground state occurs near the transition. 
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1 Theoretical basics 

1 . 1 One electron in magnetic field 

When the mutual interactions are left aside, electrons in a plane subject to homogeneous perpendicular 
magnetic field B fill the macroscopically degenerate equidistant Landau levels (LLs) with energies E = 
{n + ■7;)huj, n = 0, 1, 2, . . . The degeneracy of all the levels is the same, and it increases proportionally 
to the magnetic field. Therefore, occupancy of the Landau levels, the filling factor, depends both on the 
number of electrons A^e (per area L^) and on B: 

" eB/h Ly{27r£l) {BL^)/{h/e) ^/^^ N^' 

Note that this v, i.e. number of Landau levels occupied in the ground state is equal to the inverse number 

of magnetic flux quanta <i>o = /i/e per electron in the system (the second last expression in (1). 

These facts can easily be obtained by solving the single-electron Schrodinger equation with Hamiltonian 

Ho^:^{p+eA)\ Vx/\ =(0,0,B), (2) 

a particularly nice and understandable example of this calculation is given by Murthy and Shankar [56]. 
Suitable energy and length units are the cyclotron energy fku — HeB /m and the magnetic length yjTijeB 
denoted by Iq. 

Let us now focus on the lowest Landau level. There are eB/h- states per unit area having the same energy 
and infinitely many possibilities of constructing a basis of this space. Choosing the Landau (sym- 
metric) gauge in (2) which is translationally symmetric in one direction (rotationally symmetric around the 
origin) we are lead to the following bases 

Landau: A = (0, Bx, 0) , i/ife^ (x, y) = exp{-ikyy) cxp[-(a; + ky)'^/2£l] , (3) 

ky/{2TT/L) =0,1,2,... 

symmetric: /\ = i(j/,-a;,0), Vml^) = z'" exp(-|2|V4^o) > to = 0,1,2,... (4) 

Especially for the latter basis, formulae are often more transparent if we use a complex variable z = x + iy 
rather than x, y separately to address the points in the plane. 



1.1.1 Magnetic translations 

A plane with perpendicular homogeneous magnetic field is obviously translationally invariant. However, 
spatial translations applied to the Hamiltonian may alter the gauge even though they leave the magnetic 
field unchanged. Operators which conserve also the gauge (and which therefore commute with Hq) are 
the magnetic translations [83, 84]. These operators will thus replace the ordinary translations applicable 
to systems without magnetic field. 

Magnetic translation operators depend on the choice of the gauge, in particular for the Landau gauge (3) 

u={ui,U2): T{u) = exp{-iu^y/el)t{u) , (5) 

where t{u ) is the ordinary translation operator cxp(iu • p /ft). General explicit formula for any gauge can 
be found e.g. in the article of Haldane and Rezayi [32]. 

Note, that ordinary and magnetic translations (5) coincide for = 0, exactly as a wavefunction in the 
form (3) remains unchanged up to a constant phase under the replacement y ^ y + Uy. For u = {u^, 0) 
this is not the case and that is why we must resort to magnetic translations. 
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1.2 What to do when Coulomb interaction comes into play 

The quantum mechanical solution of one electron - or many non-interacting electrons in a plane subject to 
a perpendicular magnetic field is at the root of the integer quantum Hall effect. The basic fact is that for 
integer filling factors, any, even arbitrarily small excitation costs at least the energy huj. This gap renders 
the ground state incompressible. The fractional quantum Hall effect cannot be explained in this picture. 
For instance at filling factor = 1/3, a non-interacting system has a manyfold degenerate ground state, or, 
some excitations cost zero energy (those which involve only rearrangement of electrons within the lowest 
LL), and the ground state should be compressible. Today it is well established that the effect is due to 
electron-electron interactions which select among those states one special ground state and separate it by a 
gap from the excitations. 

Now, the Hamiltonian of the many-electron system consists of two terms; the kinetic energy (leading to 
Landau level quantization) and the electron-electron interaction. 

g = y^ + ^ly, ' , (6) 

^ 2m 47re 2 ^ rj-r, 

Consider some particular filling factor, u = 1/3 for example, and let us vary the magnetic field. Since 
ly — n/{2Tii^) — n/{eB/K), this implies changing the electron density n simultaneously. The kinetic 
energy will change proportionally to hu) cx B. The interaction energy on the other hand scales with 
1/a oc ^/B, as the mean electron-electron distance a is proportional to the magnetic length a — \/lJn = 
\J\l(yeBlti) cx (for a more thorough discussion see Yoshioka [79], Chap. 4) 

In the high field limit we can therefore expect the Coulomb interaction to be a small perturbation which 
lifts the degeneracy of Landau levels. Looking for a (high-_B) ground state at some particular ;y < 1 we can 
therefore omit the higher Landau levels and study only states within the lowest Landau level. This model 
gives qualitatively coiTect predictions (for high B) and the inclusion of the Landau level mixing leads only 
to quantitative coiTections (e.g. in the ground state energy, see Chakraborty and Pietilainen [14]). 



1.2.1 Ground states: analytical many -body wavefunctions 

It is very surprising that even though we now handle a many-body Hamiltonian (6), there are still analytic 
(correlated) wavefunctions which describe the ground state at some special filling factors. The best known 
example was suggested by R. B. Laughlin [48] earning him the Nobel Prize: 

*i(zi,...,^„) = exp(-(|zip + ... + |z„|2)/4^2)TT (^^_^.)3. (7) 

There are several beautiful physical arguments why this wavefunction must be the ground state at ;^ = 1 /3 
(in fact, an excellent approximation to it, see Subsec. 1.3). These are explained in other more detailed 
pubhcations [79] (Chap. 4), [14, 30]. Let us mention here only two basic ideas about the interpretation of 

First, the (zj — Zj)^ term makes the Laughlin wavefunction isotropic and translationally invariant. More 
detailed studies (density-density correlations, see Subsec. 2.1) suggest that it resembles a liquid. Second, 
resembles 'i'l, the wavefunction of completely occupied lowest Landau level which is the same as (7), 
only with {zi — zj)^ replaced by {zi — Zj). In fact, 

1/ = 1, GS: *i oc l[{z, - Zj) ^n.<i(f;--.)' ^ ^ q^. ^ -Q^^^ _ 

i<j i<j 

Now consider a single electron in a state i/^miz) oc z™ (4). If we pierce the system by an infinitely 
thin solenoid at zq and pass two magnetic flux quanta adiabatically through it, this state will evolve into 
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ipm{z)(z — zq}^. This leads to the following interpretation. The Laughlin state is just the completely filled 
lowest Landau level, but the constituent particles are electrons with two attached magnetic flux quanta 
rather than bare electrons. 

So far we have only spoken about the filling factor = 1/3 and fully spin-polarized electrons. Generaliza- 
tions of these concepts are possible also to systems where electrons are not fully spin polarized. Halperin 
proposed the following WFs [35] 

<^m^'n[z]= n (^^-^.r n {4-4r' U^'^-''krI[^M-\^^\v<)eM-w//<) 

i<j<N^ k<l<Ni i<N-\ i,j 

(8) 

with TO, to' odd. The state assumes N-^ (N^) particles with spin up (down) and Zi (zp describe their 
positions. The filling factors of the two components are 

m' — n m — n 

= — ; 2 ' = — 1 2 ■ 

mm — toto — n 

Thus, they describe a state at filling v = v-^ + vi and polarization p ^ {v-^ — vi)/v. For example, the 
choice m ^ m' = 3 and n = 2 leads to the total filling factor 2/5 and zero spin polarization = I'd'n). 
These analytical results will always be a good starting point for investigations going to regions where only 
numerical methods are possible. Before continuing, however, it must be emphasised, that wavefunctions 
in (7,8) are not the only analytical trial wavefunctions known in the lowest Landau level. A more detailed 
review can be found in [17] (MacDonald and Girvin). 

1 .3 Other types of electron-electron interactions 

A model of short-range interaction (SRI) in fractional quantum Hall systems is the main issue of this 
section. We will explain how a general interaction V{r) between two particles within the lowest Landau 
level (LL) can be represented by a set of Haldane pseudopotentials {Vm] [30] and show how this concept 
makes it easier to study different classes of interaction. In particular, this discussion will unveil under what 
conditions the Laughlin wavefunction (7) becomes the exact many-body ground state. 

1.3.1 Two particles, magnetic field and a general isotropic interaction 

Let us consider two negatively charged particles in a plane subjected to a perpendicular magnetic field B. 
Assume that their interaction is described by a potential (energy) V{r) which depends only on their mutual 
distance. Classically, when starting from rest, the particles would move along a straight line towards or 
away from each other were it not for the magnetic field. The Lorentz force bends their trajectories and 
makes them orbit around their centre-of-mass on a circular trajectory. In quantum mechanics, this circular 
motion is quantized just as in case of an electron orbiting around a hydrogen nucleus. Roughly speaking, 
only discrete separations between the two particles are allowed. Interaction energies V{r„T) = Vm 
rather than the full form 1/ (r), r G (0; oo) fully determine the spectrum of a many-body system of particles 
interacting via V{r). 

Let us now follow this idea in more detail and derive the precise claims. The Hamiltonian for two particles 
reads 

H = -^(Pi + |e|/\i)2 + -L(p2 + + V{\r, - ral) . (9) 

ZTO 2to 

Following Laughlin [49], we write it as a sum of the centre-of-mass (CM) and relative parts [rcM = 
{I'l + ''2)/2, rrei — {ri — r2)/\/2]- Without V{rrei), both parts will be equivalent to a single particle in 
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—mhuj 

m= 1 2 3 4 5 6 7 m= 1 2 3 4 5 6 7 

LLL n=0 ■_~-^'-.4-'--^'-V-4"~-4-"'-,- ~ 

■a \ \ ^ ^ ^ 
i=0 1 2 3 4 5 6 7 

Fig. 1: Spectrum of one particle confined to a plane subject to a perpendicular magnetic field or, equivalently, the 
spectrum corresponding to the relative motion of two particles (Hrei)- Left: without particle-particle interaction, the 
two terms, 'harmonic oscillator' (quantum number i) and 'angular momentum' (quantum number m) combine into 
degenerate Landau levels (quantum number n). Right: interaction lifts the degeneracy. If the interaction potential is 
small compared to the harmonic oscillator term (i.e. (V) ^ hu). Landau levels are roughly preserved. The energy 
levels Vm within the lowest Landau level (sorted according to m = (Lz/h)) are then the Haldane pseudopotentials. 



a plane in magnetic field. Hence, a Landau level index (A^ = ncM, n — n-rei) and an angular momentum 
(M = racM, m = rrirei) will be attributed to both parts. Out of these, N is fixed to zero (two-particle state 
within the lowest LL) and M is unimportant as it is merely tantamount to fixing the CM to some particular 
position in the plane. Eigenstates of the relative part will be sorted as shown in Fig. 1 on the left. 
With V{rrei) included, the relative part reads 

Hrel = ^ + lf>loirrel/eo)'' ~ \^K^l +V rd\) (10) 
^ ^ 

with L^^i denoting the (z-component of relative) angular momentum. This is the well known Fock-Darwin 
form, the standard usage of which is to describe one particle in a magnetic field and confining potential 
V. However, we will use it in a different way here: we consider V{r) to be weak, and for example a 
repulsive oc 1/r potential, compared to the parabolic term in Hrei,kin- It is only our initial assumption that 
the interaction is much weaker than the cyclotron energy, Ec ^ fu^ (Subsec. 1.2). 

Owing to [Hrel: ^reil ~ cigenstatcs of Hrel can still be classified by angular momentum. Moreover, 
assuming the states of the lowest LL to have no admixtures from higher LLs (Ec <C fu^), the eigenstates 

will not depend on V{rrei) at all. This is a combined effect of the Lz symmetry and the requirement of 
analyticity (confinement to the lowest LL): angular part of the form cxp(—im(/3) implies ■i/;'"(z) oc z™. 
On the other hand, energies of these states will shift differently for different 7ti's (Fig. 1, right). The 
energy shift due to the interaction between particles is Vm — (V"™; I^IV'rei) ^ state with relative angular 
momentum to. Since r„i ~ (V'^z \A'^Tei) — ^oVSto + 1 we can roughly estimate Vm to be V{rm)- 
The operator of any weak {Eq ^ huf) interaction (in the lowest Landau level) can be then written in terms 
of its spectral decomposition 

oo 

V{rrel)=Y.\'i^rel)VMll\. 
m=0 
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The spectrum and eigenstates of a many-body system confined to the lowest Landau level and interacting 
by V{rrei) is thus completely determined by the discrete set of numbers {Kn}- 

The quantities {Vm] are called Haldane pseudopotentials. They were first introduced in [30] in the context 
of interacting electrons on a sphere. Finally, we add two remarks. 

Fermions and bosons. A careful reader may have noticed that we have spoken just about two particles so 
far An additional constraint that e.g. (spatial part of the) wavefunction should be antisymmetric implies 

Therefore, only the states with m odd (11) are allowed in the case of two electrons with the same spin 
(where a symmetric spinor part implies an antisymmetric orbital part of the wavefunction). In other words: 
only the values of Vi, V3, . . . are needed when we describe motion of fully spin polarized electrons. 
Uniqueness. If V{r) is given, the pseudopotentials Vm are determined uniquely. However, the opposite is 
not true: knowing only the values of Vm, we cannot reconstruct the full form of V{r). 

1.3.2 Particular values of Haldane pseudopotentials for the Coulomb interaction 
The ideas above are not valid exclusively for the lowest Landau level. 

Let us consider a numerical example for electrons in a plane, one of them located in arbitrary Landau level 
rii and another in 712. Their relative angular momentum be m. Given these three numbers, the state is 
uniquely defined, up to the center-of-mass part of the wavefunction, as we have already stated. Assuming 
interaction of the form V{q) (in the Fourier space), their interaction energy can be written as [60] 

/>OC 

ym'-'''= qdqV{q)Lr.,{qy2)L„,{qy2)Lm{q^)exp{~q'). (12) 
Jo 

The Laguerre polynomials are defined by Ln{x) — (l/n!)[a;"e^^]*^"^e^ . For the case of Coulomb inter- 
action, V{q) ~ a/l^l, the integrals in (12) can be evaluated (easily and) analytically. Figure 2 shows their 
values for the cases (a) both particles in the Lowest Landau level (n = 0), (b) both particles in the first 
Landau level (n = 1) and (c) one in the lowest and one in the first Landau level. 

For ni = 7i2 = the coefficients Vm decay monotonically with increasing m, exactly as the Coulomb 
energy does with increasing distance. The non-monotonic structure of Vm for the case of particles in the 
first Landau level is due to the additional structure of wavefunctions in higher Landau levels (e.g. a node 

at ?■ = for n = 1). 

1.3.3 Model interactions: hard core, hollow core 

Why is a hard-core interaction {short-range interaction) important for the physics of the lowest Landau 
level? 

There are three reasons: (i) it is the strongest part of the Coulomb interaction, (ii) the Laughlin wavefunc- 
tion is an exact (zero energy) gapped ground state for this interaction and (iii) the ground state changes 
only little if the other terms of the Coulomb interaction are considered. 

Let us discuss this in more detail. For the purposes of this paper, the short-range interaction (SRI) for spin 
polarized electrons is defined by the Haldane pseudopotentials 

short-range int. (spin polarized electrons): {Vi, V3, V5, . . . } = {1, 0, 0, . . . } . (13) 

For the Coulomb interaction. Fig. 2a, V\ is indeed the strongest pseudopotential. 

On the other hand, considering the Laughlin wavefunction of N particles (7), any pair of electrons in it 
is in a state with relative angular momentum m = 3: owing to factors {zi — Zjf" ■ As there are no pairs 
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(b) 



0123456789 



(c) 



0123456789 



(d) 



0123456789 



Fig. 2: Values of Haldane pseudopotentials Vm for Coulomb-interacting electrons (a-c): between two particles (a) both 
in the lowest Landau level, , 712) = (0, 0), (b) both in the first Landau level, (ni ,712) = (1, 1)> and (c) one in the 
lowest and the second in the first Landau level, (m ,712) = (1,0). The pseudopotentials referring to electrons in a plane 
(a,b,c) are 'universal', those related to electrons on a torus (a,d) depend on its size (here Nm = 30, Subsec. 1.5.1). 
Values of {Vm} chosen for as a model for short-range interaction in this work are shown in (d). 



with angular momentum m = 1, the total energy of this state will be Vi • + V3 • N{N — l)/2 = 
for SRL Also, this state is rigid: any excitation of this state must remove the triple zero from some of the 
electrons (leaving only a single zero required by antisymmetry) thereby creating some pairs with m ~ 1. 
This implies a finite excitation gap. 

These are analytical results. A surprising numerical result is that the many-body ground state changes only 
slightly if other pseudopotentials V3, V5, . . . are 'turned on' up to their Coulomb values (Fig. 2a). This 
has been confirmed by Haldane and Rezayi [33] (later also by others, e.g. [23]) by calculating the overlap 
between the real ground state and the Laughlin state for different sets of Vm- It is also shown in [33] that if 
Vi is lowered beyond some critical value (while keeping other pseudopotentials on their Coulomb values), 
the gap collapses rendering the ground state compressible. These observations have been systematised 
by Wqjs and Quinn [75]: they argued that both Coulomb and short-range (13) interactions belong to the 
same class of superharmonic pseudopotentials where particles try to avoid low m pair states, implying 
the Laughlin ground state. A drastic change in the ground state occurs first when we leave the mentioned 
class of interactions. The 'superharmonicity' means roughly that V„i decays fast enough with growing m 
[75]. The aforementioned collapse of the gap is then no longer surprising as the decreased value of Vi will 
eventually violate potentials superharmonicity between to = 1 and 3. 

An advantageous property of the interaction (13) is that it is effectively non-parametric, the only present 
parameter Vi determines only the overall scaling of the energy scale within the lowest Landau level. 
These results can be summarized by stating that the short-range interaction is the component of a realistic 
interaction which determines almost completely the properties of the v = 1/m spin-polarized ground 
states. 

Obviously, for non-fully spin polarized systems it is not possible to keep Vi 7^ only. Electrons with equal 
spin are still closest in the state to = 1 (with energy Vi), electrons of unlike spin however are closest in 
the state to = with energy Vq. Such a model is obviously not as elegant as in the former case, it contains 
two parameters Vq, Vi whose ratio cannot be factored out of the Hamiltonian. An alternative might be the 
potential 

{yo,Vi,\/2,...} = {oo,l,0,0,0,...}. (14) 

Another typical model potential presented by Haldane and Rezayi [34] was inspired by the low value of 
Vf) in the first Landau level compared to the lowest Landau level (Fig. 2). They suggested the hollow-core 
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Table 1: Values of Haldane pseudopotentials corresponding to particle-particle interactions of the type V{q) — q^''. 
These values are additive, e.g. V{q) — — + 1 corresponds to the 'hollow core interaction': {Vm} = 
{0,1,0,0,0,...}. 

potential 

hollow-core interaction: {Vq, Vi, V2, ■ • ■ } = {0, 1, 0, 0, 0, . . .} . 
and tried to explain the even-denominator fractional quantum Hall effect at = | using this interaction. 



1.3.4 An alternative definition of Haldane pseudopotentials 

Haldane introduced the quantities Vm originally for interacting electrons on a sphere [30]. In that case, or 
for electrons in a plane, m can be identified with the relative angular momentum of the electron pair, having 
in mind that m is closely related to the average separation between the particles increases. In contrast to 
that, rotational symmetry of the configuration space is lost on a torus and angular momentum is no longer a 
good quantum number. Here we will introduce an alternative definition of Haldane pseudopotentials which 
is applicable also for particles on a torus [74]. 

First, recall that matrix elements of the Coulomb interaction on a torus can be conveniently evaluated in 
Fourier space (Subsec. 1.5.3) where 

Vir)^ey\r\ ^ Viq) = ey\q\. 

Second, consider a general (radial, bounded) interaction with its Fourier transforms V = V^(|<7|) and 
expand V^d^l) into a Taylor series. Owing toV^(r) = V{—r) = V{\r\), the series will be free of odd 
powers (72'=+!. Now, go back to the direct space and use .F[/(r)]''^' = [iq^ T f [r) 

V{q)=vo + V2q^ + VAq^ + --- => V {r) =VQ5{r) -V2^^5{r) +ViVH{r) - . . . . {\5) 

The coefficients Vi now fully characterize the particle-particle interaction. In an extension to the article of 
Trugman et al. [74], let us show how to translate them into y^'s, i.e. interaction energy of the two-particle 
state in a plane (or on a sphere) with relative angular momentum m. 

For the evaluation of V„i = {ijjm\V{r)\ipm), let us take the functions from the planar system, (11) 
plus normalization. If V{q) = q^^ then 



2ku.. (-1)" 



1 d d 
r dr dr 



(16) 



This is a unique prescription of how an interaction of the type V{q) = q^'^ can be transcribed into the terms 
of Vm- Table 1 contains these coefficients for several lowest powers of q. Note that Vm = for m > k. 

In conclusion, an interaction potential defined by some particular set of values of Haldane pseudopotentials 
Vm can be recalculated into the coefficients Vi in (15) (Taylor series of V{q)) using Table 1 or, more 
generally using (16). 
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Again, several remarks should be made. 

(1) The expansion in (15), being first suggested by Trugman and Kivelson [74], looks unusual. In the 
distributional sense, we say that a non-zero ranged potential V{r) can be written as a sum of terms with 
'zero range'. 

Instead of a (5-function imagine rather a sharp peaked function 6^, a Lorentzian of width b, for instance. 
Functions V^'^5f,(r) will then have 'the longer the range the higher the k is': it is instructive to draw a 
sketch of the first few functions ((5;,)'^'^'-' and consider as 'range' the position of the local extreme which is 
the most distant from the origin. In this sense, (15) is an expansion of V{q) in terms of 'increasing ranges'. 

(2) When calculating the Coulomb matrix elements for particles on a torus (Subsec. 1.5.3), we do not use 
the full function V{q) but only its values in discrete 'lattice' points q . This is obviously due to the periodic 
boundary conditions. In particular, q = is missing among these points. 

Thus, we need not worry about the long-rangedness of the Coulomb potential, V{q ^ 0) ^ oo which 
renders it unexpandable into a power series of q. Instead of l/g we may imagine to have considered any 
other polynomial in q which matches the values of l/q at the 'lattice' points. Both interactions must lead 
to the same results. 

(3) Example: consider two electrons in the lowest Landau level interacting via V{q) = aq^. Eigenstates 
sorted according the to increasing value of the particle-particle distance (r) may be indexed by an integer, 
say TO. The state m = will have an energy of —a, the state m = 1 will have an energy of a and all other 
states (with larger interparticle separation) will have zero energy. 

The state with to = will have a symmetric wavefunction and will be thus prohibited for electrons with 
equal spins. Thus there will be only one state with non-zero energy for this case and it is the state with the 
lowest interparticle separation. The potential V{q) = aq^ defines therefore a hard-core interaction. 

1 .3.5 Short-range interaction on a torus 

The decomposition of the Coulomb interaction (in the lowest LL) into the set of Haldane pseudopotentials 
has already been shown in Fig. 2. This is also the spectrum of two Coulomb-interacting particles on a 
sphere. 

Let us now consider a pair of particles on a torus. Fig. 2a. The index m is no longer the angular momentum 
of the pair as this is not a good quantum number. The wavefunctions ip^li{z) oc z™ exp(— |zp/4^§) in 
(11) must be modified, in order to comply with the periodic boundary conditions. 

In Fig. 3 we show some of the wavefunctions corresponding to the relative motion on a torus of size 
N,n = 30 (= a6/27r^Q, Subsec. 1.5.1). We will denote them simply by -0™, m = 0, 1, ... , m,nax and 
skip all other indices which would be appropriate, e.g. to indicate that they depend on the size of the torus 
(a X a). Even though these states are more complicated than those in (11), they can still be sorted according 
to growing values of r,„ = {ij}"^\r\ip^"'). It is not surprising that the states -0™ for low m (<C Nm), Fig. 3b, 
look very similar to the eigenstates of angular momentum to for infinite systems (11). First when r™ 
becomes comparable to the system size, deviations from the circular form of occur (middle column of 
Fig. 3b). It is an intriguing property of the periodic boundary conditions that the states with very high to 
look very similar to those with very low m. If we fix one electron to r = (0, 0), then the second electron 
orbits around (0,0) at a distance in the state ijj™', whereas in the state ■0™™"^-™ it orbits around 
(a/2, a/2) at the same distance. This can be seen by comparing the left and right columns of Fig. 3b. 
Now, with 0™, as a substitute for the relative angular momentum eigenstates, we can define Haldane 
pseudopotentials on a torus by Vm ~ ('0™l^(''rei)IV''")- Their values (Fig. 3a) are almost equal to Vm 
in a plane, as long as ip"'- is not affected by the periodic boundary conditions (Fig. 2a) i.e. for small 
values of to. 

A reasonable model mimicking the short-range interaction keeps the first two energies of the spectrum in 
Fig. 2, i.e. the pseudopotentials Vq, Vi at their 'Coulomb' values while setting the other ones to zero. 
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Fig. 3: Two-particle eigenstates of Coulomb-interacting particles in the lowest LL confined to a torus. States are 
shown irrespective of the symmetry or antisymmetry of the wavefunction. Left: spectrum (horizontal axis has no 
meaning: levels are 'randomly' distributed into four groups in order to show degeneracies). Right: relative part of the 
wavefunction, j?/)'"(r,.e;)p, for several states. Here, r^ei ~ corresponds to the corner of the square (the four corners 
are identical due to periodic boundary conditions). Some of these states resemble the eigenstates of relative angular 
momentum, see text. 



Table 1 gives a prescription how to encode such an interaction into V{q). We thus arrive at an interaction 
potential defined by V{q) = 0.34(j^ — 1.51 which is used throughout this work to model a short-range 
interaction unless something else is explicitly stated. 



1 .4 Composite fermion theories 

Let us recall the observation from Subsec. 1.2: three zeroes are bound to each electron in the Laughlin 
state 'I'l. A good way to see this is to fix the positions of Z2, • ■ • , z„ and use the last 'free' coordinate zi 
to inspect (the zeros of) the wavefunction. One zero is required by the Pauli principle (when zi = Z2, the 
wavefunction must vanish), the others are 'voluntary'. 

Whenever an electron goes once around a zero in 'i'l, the wavefunction acquires a phase equal to the 
Aharonov-Bohm phase corresponding to one magnetic flux quantum. From this point of view, the Laughlin 
state can be interpreted as the v = 1 state where two magnetic flux quanta are attached to each electron. 
These objects (electron dressed by two flux quanta) are called composite fermions (CF). Note however that 
the precise definition of a composite fermion may vary in different theories as will be explained below. 
Intuitively, this concept explains the existence of a gapped ground state at filling factor v = l/i. Originally, 
there are three flux quanta per electron (1) and the huge Hilbert space of many-electron states in the lowest 
Landau level is completely degenerate without interaction. In other words, we expect no gap without 
interaction. If we now assume, that the Coulomb interaction leads to the formation of composite objects, 
an electron and two flux quanta, then there remains only one free flux quantum per CF. This in turn implies 
the filling factor of vcf = 1 for CF (1). We know that in this case the ground state of particles obeys 
Fermi statistics (see comment [1]) and is gapped. If a Landau level is completely filled, then any, even 
infinitesimal, excitation requires promoting at least one CF into a higher CF Landau level which costs the 
finite energy > Hojcf- Now let us describe some of the current composite fermion theories in a little more 
detail. 
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1 .4. 1 Chern-Simons transformation 

Looking at the Laughlin wavefunction in the way sketched above, we might find it reasonable to incorporate 
the flux attachment into the Hamiltonian. 

The Chern-Simons (CS) transformation is just a gauge transformation of the magnetic field 

acs{r) = / d?r^^-^^^^!—^-^\r,)^{r^) . (17) 

It does not change the magnetic field ('gauge transformation') felt by the electrons only owing to the fact 
that two electrons cannot be simultaneously on the same place (see comment [2]). The price for this is 
that the transformation is singular, a cs diverges for r = r,. The objects '^^{r) are the one-electron field 
operators and a is the number of attached magnetic fluxes. 
After this transformation the full Hamiltonian 



i y"d'r*t(r)[-iftV, + e/^(r)-e3cs(0]'*(0 



(18) 



contains - apart from one-particle terms - two-particle terms (those containing a cs) and also three-particle 
terms ^'t(r)5'(r)5'^(ri)^'(ri)^'t(r2)^'(/'2) (they originate from a ^g). The CS transformation alone thus 
does not really simplify the original Hamiltonian. 

A mean field approximation can be made at this point where the density operator '^^ {r i)'i> {r i) in acs 
is replaced by the mean value ns- We arrive at a single particle problem with an effective magnetic field 
BcF = B — a(t>ons- In illustrative terms: 

CS transf. mean tield 

a many-body system — * a very complicated many- — * a simple one-particle 
at J/ = 1 /3 body problem at = 1 /3 problem at = 1 

The final one-particle problem alv = \ has a non-degenerate ground state, the lowest Landau level fully 
occupied by CF. By means of this procedure we thus circumvented the original problem that the Coulomb 
interaction must select the ground state out of the vast number of degenerate many -body v — l/i states 
within the lowest LL. 

A mean field approximation is not the only possible treatment of the Hamiltonian (18). However, theories 
beyond the mean field i.e. those treating fluctuations of the gauge field, are very complex [51]. 
Using the CS transformation we attach 2s ^ a vortices (not zeroes) to each electron. In the mean field 
approximation the problem is equivalent to non-interacting particles in reduced magnetic field Be f which 
then corresponds to a filling factor vcf- It turns out that many of the experimentally observed fractions v 
(exceptions see in [59]) correspond to integer vcf- Let us conclude with an overview of relations between 
quantities referring to electrons and to CF, cf. (1). 

BcF = B{\ - 2su) = B- 2sns$o , t = tcF - ^° , (19) 

Vl — 2si^ 

1 1 „ P ■ . 

= 2s , 1/ = , (p, s integer). 



VCF V 2sp + 1 



1.4.2 Composite fermions a la Jain 

Compared to the Chern-Simons transformation, Jain chooses to go in some sense the same path but in the 
opposite direction [38, 39]. It starts with a wavefunction of particles (fermions) at integer filling vcf = P, 
attaches s zeroes {not vortices) to each particle and, after projection into the lowest Landau level, it presents 
the result as a trial wavefunction for the ground state at filling ly = p/{2sp + 1) (19). This procedure 
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Table 2: The scheme of construction of Jain's wavefunctions for CF with two flux quanta attached: examples of 
composite fermion filling factors (p| , pi are numbers of fully occupied spin up and spin down CF-Landau levels) and 
corresponding electronic filling factors. 

reproduces exactly the Laughlin wavefunction and at other fractions it gives wavefunctions with very high 
overlap with ground states calculated numerically by exact diagonalization. 

There are two central reasons why this approach is very popular. On one hand, it gives a simple single- 
particle picture of what is going on in the highly correlated many-body problem. On the other hand, it 
offers explicit formulae to work with since it is easy to write down a wavefunction of p full Landau levels. 
A very pleasant feature of this approach is that it allows to incorporate the spin of electrons easily [77]. 
Take p-\ of full Landau levels with spin up and pi of full Landau levels with spin down. These Landau 
levels are then called composite fermion Landau levels. The magnetic field felt by the CF, i.e. the field 
corresponding to filling factor vqf = p is called effective magnetic field B^ff. It is weaker than magnetic 
field B corresponding to the electronic state at (19). 

Note, that the filling factors in (19) are all in range i' < ^. For i < ly < 1, Jain et al. [77] suggest the 
idea of antiparallel flux attachment: the effective field B^ff is antiparallel to the real field B, however, 
the additional flux quanta are added in parallel to B i.e. antiparallel to B^ff- In terms of (19) this means 

p —p 01 V = p/ {2sp — 1). 

An example of candidates for ground states and their polarization provided by Jain's composite fermion 
theory is given in Tab. 2 (see Chakraborty [13] for a review regarding ground states with various spins). 



1.4.3 Composite fermions a la Shankar and Murthy (Hamiltonian theory) 

The Hamiltonian theory of FQHE (Shankar and Murthy [56]) builds on previous works of Jain and those 
concerning the CS transformation, quoting words of its authors, it combines the strengths of the both 
theories. 

It provides a projected Hamiltonian of the lowest Landau level which scales only with the Coulomb inter- 
action. In addition to each electron a new independent object is introduced; a pseudovortex. Its definition 
on the level of commutation relations (Eq. 129 in [56]) assures, that if an electron goes around a pseu- 
dovortex, it picks up the phase of 27r 2s i.e. it has the same effect as an insertion of 2s flux quanta. Note 
however that it is yiot a zero of the wavefunction. The projected Hamiltonian is written in coordinates 
which are a combination of the electron and pseudovortex position (Eq. 138 in [56]). This combination is 
then called composite fermion coordinate. 

For this Hamiltonian an ansatz for a ground state can be written down. At filling v = p/ {2sp + 1), it is p 
Landau levels filled with CF. It is then possible to evaluate their Hartree-Fock energies. 
The first substantial success of this theory is that it produces the correct scaling of spectra within the lowest 
Landau level (cx ^/B). Compared to Jain's theory, it keeps track of the fact that the two fluxes (which sit 
exactly at each electron in the Laughlin state) can be only loosely bound to electrons. This is owing to the 
dynamical degree of freedom given to the pseudovortices. On the other hand, the electronic coordinates 
are actually the only really independent ones, for instance in the Laughlin wavefunction, all the zeroes 
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{zi — Zj) are expressed in terms of electronic coordinates. Thus, the price we must pay for the extension 
of the Hilbert space is that we must perform a projection to the space of physical states at the end. 
Nevertheless, this does not seem to be a substantial problem and thus the Hamiltonian theory is probably 
the most advanced achievement in an effort to understand the many-body physics in the FQHE. 

1 .5 How to test the CF theory? 

The concept of flux attachment (Sec. 1 .4) provides a well-understandable model of the FQHE. However 
transparent it seems at the first look, predictions based on it must be tested against a model which contains 
less approximations. Exact diagonalization (ED) is a good choice for this purpose. In exchange for exten- 
sive numerics to perform, the only substantial approximation of the method is to take a finite instead of an 
infinite system. 

The main part of this Section concerns the exact diagonalization (ED) [14]. It is definitely not the only 
numerical method used in the context of the FQHE. Some numerics is at the end of nearly any method 
as soon as many-body problems are concerned, be it a Hartree-Fock treatment of CFs or Monte Carlo 
simulations of the Laughlin state mapped onto a one-component plasma. O 

We take the complete many-body Schrodinger equation but confine the interacting electrons moving ac- 
tually in an infinite plane onto a compact (i.e. finite-sized) surface, possibly without edges. The standard 
choices are a sphere [30], a torus (square with periodic boundary conditions) [80] and a disc [48] (see 
Yoshioka [79] for an overview). Although these manifolds are locally flat and therefore with growing 
system size a convergence towards infinite-plane results can be expected, they all break some of the sym- 
metries of the infinite plane. For instance, the sphere keeps the angular momentum while the torus retains 
the translational symmetry. In any case, the hope is that effects inflicted by the finite size can be separated 
from those generic to a two-dimensional electron gas. Another usual yet not necessary approximation is 
to neglect Landau level mixing, i.e. restriction to the lowest Landau level only. Also note, that there is a 
long way from an ideal 2D system which study here, to the experimental reality (impurities, effective mass 
approximation, finite thickness of the 2D electron gas etc.). 

1.5.1 Torus boundary conditions 

One possibility to model an infinite plane by a finite manifold without edges is a rectangle with area a ■ b 
with periodic boundary conditions (PBC). Topologically, this is the same as a torus, although it is better to 
stay with the former picture for the sake of twisted PBC, even if we sometimes use the word 'torus' as a 
shortcut for this model. 

What are the single particle states of the lowest Landau level in this case? Recall (3) where single-particle 
states complying with translational symmetry along y are given 

^o,fc; {x', y') = cM-^Ky') cxp[-(x' + k'y f/2] , 

primed variables are in units of magnetic length, x' — x/Iq, k' = fc^o- Periodic boundary conditions 
along y admit only discrete values of k'y = {2Tr£o /b)j with j integer The wavefunction is centered in the 
a--direction around Xj — kyi^ and if we require Xj to lie within [0; a), we have < —k'y < a/£o. Thus, 
up to a sign, 

0<,<^.m. (20) 

Equation (1) with = ab implies that a6/27r^Q is equal to the number of magnetic flux quanta (<i>/$o) 
which pass through the rectangle and by virtue of (20) it must be an integer. This brings us to the central 
insight that there is only a finite number m = Nm of states in a square with periodic boundary conditions 
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(subject to magnetic field and discarding all but the lowest Landau level) and that the size of the torus 
(area in units of £q) can be measured by the number of magnetic flux quanta iV,„ penetrating the torus: 

ab = 2Tr£lN„i. 

States i/^o.fc' {x',y') shown above are not periodic in the x direction and this can be accomplished by 
periodic continuation: tf'ix, y) fpix, y) + i!{x + a, i/) + . . . . The (non-normalized) single particle states 
we will be dealing with are thus [80, 78, 81] 



fj{x',y') = exp 

A:— — OO 

j = 0,l,2,...,m-l. (21) 
These states constitute the single-particle basis of the lowest Landau level. 



,y'(^L+k)C-l(x'-i^+k)C 

m 2 \ m 



C= \/- •27rm, 



Twisted boundary conditions 

Consider what happens if we require the modulus of ip rather than ip itself to be periodic, similar to 
Bloch's theorem. Thus, the wavefunction may acquire a non-trivial phase when going once around the 
torus. Mathematically, this can be described using the magnetic translation operators (5) 



T{aex)ii = cxp{i4'x)ip , T{hey)ii = cxp(i0j^)?/! . (22) 
Fixing phases <t)x, (py, the correct (non-normalized) periodic single particle states are 



k— — oo 



(Pj{x,y) = Y exp{ik(l)x)t{kaej,)exp{-iXjy/il + i<f)yy/b)cxp[-{x- Xjf/2£l], 

j = 0, 1, . . . ,TO - 1 , 



■' m 



(23) 



where t{£_ex) is an ordinary translation, i.e. an operator transforming -0(2;, y) into 7/'(-'J;+Ci y)- For (f>x,4>y = 
the original result (21) is recovered. This choice of cpx, (py is also used throughout this work. 
Interpretation of (px, <py By imposing the PRC we arrived at the statement that wavefunctions must be 
centered (along x) at Xj ~ (a/m) ■ j, j ~ 0,1, . . . There is no a priori reason for the point x = Xq = to 
be more important than x = Xq,^ = (a/m) ■ 0.5 which is not among the just mentioned X/s. By varying 
(px, the set {Xq,Xi, . . .} — {a/m){0, 1, . . .} is transformed into {a/m){0 + (px/2'K, 1 + (px/2Ti, . . .}. 
Thus, sweeping (px from to 27r, we probe all points between and a in the x direction. Independently 
on this, we may sweep through all ky points in the interval [0; 27r/6] by changing (py. Thus, (px and (py are 
analogous to lattice wavevectors within the first Brillouin zone in an ordinary periodic system defined by 
ordinary rather than magnetic translations. 

In summary, by considering only a finite system, we have only m states to probe the whole plane (i.e. [0, a] 
in X and [0, 27r/6] in y). Sweeping (px, (py from to 2ti we can access an arbitrary point in the plane. 
Another interpretation of (px , (py was given by Tao and Haldane [73] in terms of additional magnetic fluxes. 
These come from two ideal anuloids (closed solenoids): one goes inside the torus and another around the 
torus outside. {h/e){(px.y/2Tr). It was also shown [32, 37] that (px increasing linearly in time acts as a 
homogeneous electric field in x direction. 



General basis of single-particle states on a torus: complex coordinates 
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A precise discussion of one-particle states on a torus including the phases 4>x , (j^y was first given by Haldane 
and Rezayi [32]. They showed that the most arbitrary state is 

m 

V'(a;, = exp(— ix^) • exp(ifcz) ?9i (7r(z — , z = x + iy (24) 

1=1 

analytic 

where 'i?i(w|r) is an elliptic theta function ([29], p. 921), fc is a real number in the range < irm/b and 
Zi are some fixed complex numbers within the rectangle [0, a] x [0, b]. In the terminology of (23), these 
states correspond to any j and any (fi^, <j>y. The most important things to know about the theta functions are 
that it is analytic, that 7?i (z — z/ |i) cx z — zi for |z — z; | and that z/ is its only zero in the rectangle. In 
this form, it is also clear that m is equal to the number of flux quanta in the elementary cell (the rectangle). 
Going once around the rectangle, the wavefunction gathers a phase of 2ttx number of zero points inside. 
That number is just m, each factor in (24) contributes by a single zero. 

By choosing fixed 4>x, (py, there arise m possible choices for the values of k and zq = J2i J ~ 

0, 1, . . . , TO — 1. For each pair (fc, zq) we can construct one function of the form (24) and the resulting to 
functions will constitute a basis of the lowest Landau level, just as the basis in (23). There is naturally a 
large freedom in choosing one particular basis. This happens by choosing some particular position of the 
zero points z;'s while observing the constraint on zq. The basis in (23) can be obtained from (24) by putting 
the zeroes on a line, zi = i ■ bl/m + j /ma and choosing fc = {2tt /b)j for the state with (px = 4>y ~ 0. 
Even though it is by far not obvious in (23), Fig. 4 shows a 2D plot of one of such functions. 
In principle, the wavefunctions in (24) are very similar to those obtained in the circular gauge (4) except 
for substituting z by ■di{z\i). This is a manifestation of the fact, that even on a torus, circular symmetry 
is approximately preserved at short distances and deviations occur first when '&\[z\i) deviates from z at 
larger distances. One could say, ■di{z\i) is the function /(z) = z adapted to the torus i.e. deformed to 
comply with periodic boundary conditions. 

On the other hand contrary to the infinite plane, each single-electron wavefunction on a torus has as many 
zeroes as there are flux quanta passing through the torus. 

1.5.2 Many-body symmetries on a torus 

Center-of-mass 

What changes if we consider n-body states instead of single -particle ones [32]? Given the considered 
Hamiltonian (33), the most obvious symmetry is the separation of center-of-mass and relative part of the 
wavefunction 

^'(Zi, . . . , Z„) = *CA/(^)V'rei , Z = Zi + . . . + Z„ . (25) 

The center-of-mass part is just a one-particle wavefunction. Hence it must have the form shown in (24). 
Haldane and Rezayi [32] showed that it has q zeroes in the region [0; a] x [0; b] for filling factor v = 
Ne/{qNe). Again as for single-particle states, there are q basis states for Since the energy does not 

depend on the center-of-mass position in a homogeneous system, these three states will lead to degenerate 
many-body states, provided iprei remains the same. 

This introduces a delicate topic. The electron density in a given state depends on the center-of-mass part 
of the wavefunction. Different choices of bases in the g-fold i.e. threefold for ly ~ 1/3, degenerate space 
of center-of-mass wavefunctions may lead to a g-tuple of states with practically homogeneous density in 
some cases or with quite strongly varying density in other cases (Fig. 5). This is true in spite of that we 
always describe the same ground state subspace. Even worse, in homogeneous systems we often want to 
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Fig. 4: One possible one-particle state on a torus pierced by three flux quanta (i.e. m = 3). 



study only the relative part of the wavefunction, which must be the same in all cases. If it is for example 
the Laughlin wavefunction, we know that it leads to a homogeneous density. The central trouble is then 
that the Hamiltonian eigenstates obtained by exact diagonalization contain 'i'cM- 



Relative-motion part of the wavefunction 

The discussion in the previous paragraph is based on (magnetic) translations of the center-of-mass Tcm (u). 
In an n-body state, the translation of a single (i-th) particle (Subsect. 7.2 in [14]), ti{v), can be split into a 
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(a) A basis leading to more inhomogeneous densities. 
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(b) A basis leading to less inhomogeneous densities. 

Fig. 5: Two different bases for 'i/cM, (a) and (b). At filling = 1/3 there are three allowed CM states on a torus. 
They are labeled a, /?, 7 in this figure. For each element of each basis we show the modulus of '^cm, the density of 
the corresponding Laughlin state with six electrons, i.e. the state *I'ca/^l, and section of the density along x and 
along y. Note the positions of the three zeroes in different 'I'ca/ 's (marked by the red arrows). 
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Fig. 6: The first Brillouin zone for relative translations on the torus (square with periodic boundary conditions). Its 
size depends on the number of particles in the system, cf. (27). At filling factor u — p/q with A^e = Np particles 
(and Nm. ~ Nq fluxes), number of allowed /c'^-points is A^^ and the upper right corner has /c'' = 7(71,71), 7 = 
N/2-Kq /£o . Different colours indicate points of different symmetry (or degeneracy number of a state with this k ^ 
in a homogeneous system), the rightmost figure shows the limit of large A''. 



translation of the center of mass Tcm{}^v) and a relative translation Trei,i{v — ^v). Owing to the indis- 
tinguishability of particles, the effect of the relative translation ti{v) on a particular many -particle state is 
the same for any i. We may thus omit the index and imagine i = 1, for instance. 
Again, as in Bloch's theorem, wavevector can be attributed to these relative translations [31] 

Treliv)^ = CXp{tk'- ■ V)ij . (26) 

Since Trei{v) commutes with the Hamiltonian (33), the Hamiltonian eigenstates can be sorted according 
to values of k^. In Bloch's theorem, the allowed translations are given by an arbitrary lattice vector v. Not 
all of them are allowed for Trei though [31]. 

This concept is very similar to a single particle in a periodic potential. However, there is no real periodic 
potential in an infinite plane and we introduced one particular period artificially. The largest period possible 
within our model is the size of the rectangle. 

The Brillouin zone for k"^ is rectangular (Fig. 6) and its size grows with the size of the elementary cell. For 
filling factor v ^ p/q {p,q with no common divisor > 1) and number of flux quanta per cell Ng = Nq, 
the allowed values of /( are 

k'^eo = t) ' \t\ < N/2 and integer (27) 

The quantity A is the aspect ratio. For the sake of comparison between systems of different sizes we will 
sometimes use size-independent units for k^, where k^ = (tt, tt) means the upper right corner of the 
Brillouin zone, i.e. s ^ t = N/2. 

It can be verified ([14], p. 169), that application of the operator 

CDW: cxp(ic7 • r,) , or SDW: ^ 5*+ cxp(iq • r,) (28) 

i 3 

to an arbitrary state corresponding to k^ increases its wavevector /c '' by q . On the other hand, the operator 
(28) generates a charge-density wave (spin-density wave) with wavevector q , as can be best verified by the 
simple example of the Fermi gas. Isotropic states are supposed to have k^ ^ Q. 
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The wavevector for states on a torus is also related to the angular momentum of the corresponding 
states on a sphere or on a disc, \k'''\ = {\L \/h)/R, where R is radius of the sphere [30, 14]. The direction 
of L (or alternatively L^, for instance) is related to the direction of k^: for example a plane wave going 
around the equator specified by the wavevector k"^ will have L pointing to the pole and L, = h\k^\R. 
This correspondece allows to directly compare spectra for FQH states obtained for different boundary 
conditions [33] and this, in turn, helps to sort out the finite size effects. 

Momentum 

So far, we have introduced two sorts of translational symmetries of states on a torus. One of the center-of- 
mass part of the wavefunction and another of the relative part. Since the corresponding magnetic translation 
operators commute with the homogeneous Hamiltonian, it would, in principle, be possible to spUt the basis 
of the whole lowest Landau level into several smaller bases and diagonalize in the subspaces separately. 
Each basis would be characterized by a particular value of kcM and k^. 

This procedure can help to treat larger systems but it costs some extra effort to implement it and moreover 
it is only possible in homogeneous systems. We will now discuss another of Hamiltonian's symmetries, 
described by a new quantum number J, which is a combination of the previous two. This symmetry is 
preserved with a certain class of inhomogeneities and it can be implemented straightforwardly. When 
constructing the basis for a particular value of J, we only have to select the matching Slater determinants 
(31) rather than to construct linear combinations of them. 

The homogeneous Hamiltonian in the Landau gauge (i is the particle index, Vint is the Coulomb interaction 
between particles) 

i 

conserves the total momentum in y direction. Due to the periodic boundary conditions allowed values of ky 
are (27r/&)j, j = 0, 1, . . . , m — 1. In an n-body state constructed as a Slater determinant of single-electron 
states (fj. (21), the total momentum along y is 

n 

{b/2TT)Ky = (V27r) ^ fcj, = ji + . . . + j„(mod m) = J . (29) 

i=l 

Values of J thus range for instance from to m — L It is useful to keep in mind, that ji is (up to the factor) 
the point in x-direction at which tpj. is centered, Xj. = {ji/m)a. Thus, J can also be interpreted as the 
x-coordinate of the center-of-mass of the n-electron state. 

Without proof, let us now present the precise connection between J and the wavevectors following from 
TcM and Trei (i.e. kcM, k^). Let u ~ p/q (p, q with no common divisor > 1) and m = Nq the number 
of flux quanta per cell. An arbitrary J can be decomposed into two parts 

J = JcM ■ N + Jrei , \Jrei\<N/2, JcM integer, (30) 

i.e. Jrei is J modulo N and Jcm is J divided by N . The quantity J^ei is directly the y-component of k"^ , 
more precisely, J^ei = t or J^ei = N/2 — t in (27), the former for — 1) even, the latter for pq{n — 1) 
odd [14]. 

JcM distinguishes states which differ only in the center-of-mass coordinate. By a successive application 
of TcM to one state 5* we can go through all possible values of Jcm = 0,1, 1. 
Each subspace with definite ,/ contains states of all different k^.. Since Trei by allowed translation vectors 
commute with the total momentum along y, it is in principle possible to split a basis corresponding to a 
particular J into subspaces with k^^ioy^ N/2t: = —N/2, . . . N/2. However, the basis state will no longer 
have the simple form of antisymmetrized product states of ipj (21). 



{x',y') = {x/£o,y/io) 
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1.5.3 Exact diagonalization 

Many interacting electrons in a rectangle with periodic boundary conditions can be described in the fol- 
lowing way. 

• Choose the number of flux quanta penetrating the rectangle (m). All allowed single-particle states ipj 
are those written in (21) or (23) for nontrivial boundary-condition phases (p^, <t>y Their number is m. 

• Construct all possible linearly independent n-particle states (for the given number of flux quanta 
m). Most conveniently, these can be antisymmetrized products (Slater determinants) of n states (pj. , 
denote them by 

=4---«LlO>- (31) 



• The filhng factor is then v = n/ra,cf. (1). 

• Take an arbitrary many-body Hamilton operator and calculate its matrix elements in the basis \ ji . . . jn)k, 
k = 1, . . . ,N. The dimension of the matrix is = ( " ) . 

• Diagonalize the Hamilton matrix. Eigenvalues are the total energies Ei, eigenvectors Vi ~ {vj, . . . , vf^) 
are related to the many-body eigenstates by 

N 

H\^,) = E,\^,) , 1^.) (32) 

fc=i 



This procedure is exact if we consider a system where electrons in the lowest Landau level form a periodic 
system. The approximation rests therefore in representing an infinite system by a periodic repetition of 
a representative finite cell, a procedure which has been very successfully applied in condensed matter 
theory. Formulated in other words: the Hamiltonian is exact and all approximations are implemented by 
the choice of the basis. The dimension of the matrix is finite by construction, no cutoff for one-particle 
states is needed. 

In the rest of this Subsection we present the particular form of the Coulomb matrix elements ([80] or [14], 
Sect. 5.1). 

The exact Hamilton operator in first and in second quantization is 



H 

H 



e 

Aire 



(33) 



= + -^Ji J2 'h ji aji , 



Jl ,32 
J3'J4 



where create single-electron states. The latter expression assumes already periodic boundary conditions. 
The first sum is the Madelung-type energy of the electron interacting with its own periodic images [12] 



W 



2- J2 ^-iH'^ + ^'/A)) 



(^„(z) = / dte-'V. (34) 



If only the electrons were considered, this energy would diverge at least as J^n 1/'^- ^'^ keep it finite, a 
neutralizing positive background must be considered [12]. 
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Choosing the single-electron basis according to (21) or (23), the interaction matrix elements are 



1 ^ Viq 



Aireio 2m 



+j2,j3+j4, 



■ exp 



gx=(27r/a)3 
<jH = (27r/6)t 



4 <-0 



(35) 



X exp [-2TTit{ji - x a{ji + .72 - js - jiAy)- 



with both integrals taken over the rectangle [0; a] x [0; 6]. Primed Kronecker 5 compares the two arguments 
modulo m. The last factor a is solely due to the boundary condition phase (j)y 

a(A J, (j)y) = ^,7,0 + SAj,m exp{i4)y) + dAj-m cxp(-i0y) , 

and the matrix elements do not depend on (j)x, only the basis vectors | ji . . . jn)k do. 
The periodic continuation of the Coulomb interaction in two dimensions is given by 



c2 y > \r\ 



1 27r , , /27r 27r \ ^ 



a5 ^ |q 



(36) 



hence the Fourier series of y(r) used in (35)hasy(qf) = l/|c7|- 

The Hamiltonian (33) assumes spin-polarized particles. Its extension to particles which may have different 
spin is straightforward, since the Coulomb interaction conserves spin [85], 



(37) 



31 ,32 
33,34 



operators aj^ must be extended appropriately. They create a particle in state (pj either with spin up or spin 



down. 



1.5.4 Symmetries and choices of bases 

Regarding the structure of the basis of our choice (31) there are two Hamiltonian symmetries which are 
easy to use: conservation of J, total momentum along y (29) and conservation of the z-component of the 
total spin 5*2. 

'Easy to use' means here that the basis of the whole lowest Landau level in the form of Slater determinants 
I (jiCTi . . . jncrn)k) Can be sorted into groups corresponding to particular values of J and S^. 
Sorting according to J splits the basis into m subspaces of approximately the same size w ( ) /m. Util- 
isation of Sz brings a smaller profit, since the Sz = subspace is larger than the Sz ^ n/2 subspace by 
a substantial factor of about (^" j) ~ 2"^^/\/27m. The size of the largest group is then not simply the 
number of all states divided by the number of subspaces. 

Other symmetries of the homogeneous Hamiltonian would correspond to conservation of the total spin S"^ 
and conservation of (Subsect. 1.5.2). The eigenstates of these operators, however, are generally not 
of the simple product form (31), but they are linear combinations of such states. More importantly, these 
symmetries are gone if inhomogeneous systems are considered. However, suitably chosen inhomogeneities 
can preserve the 'easy-to-use' symmetries mentioned previously (Subsect. 3.1) while still lowering the total 
symmetry of the Hamiltonian. 

If the aim is to choose n as high as possible, then the largest accessible systems have about ten electrons. At 
filling I' = 1/3 with J symmetry employed and Sz = n/2, the basis counts 1 001 603 elements for n = 10 
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and J = 5. The largest bases used in this work contained 5 x 10^ elements, extremely elaborate programs 
can handle bases up to sizes about an order of magnitude larger [52]. An alternative to the classical exact 
diagonalization is presented in Subsect. 1.5.5 

Particle-hole symmetry 

Particle-hole symmetry provides a mapping between systems at fillings v and 1 — ly (spinless electrons) 
or ly and 2 ~ ly (spinful electrons). The mapping is exact provided Landau level mixing is absent. As 
an illustrative example consider fully polarized electrons i.e. only the lowest LL with spin up is relevant, 
lowest LL spin down and all higher LLs are so far in energy that they can be neglected. The spectra of 
a ly = 1/3 and 2/3 systems are identical up to a constant shift and the corresponding wavefunctions are 
related by a simple transformation. 

Think of the lowest Landau level as of a ID chain. The Landau gauge is particularly illustrative for this as 
the one-electron states (3) are localised along x. The basic idea of the particle-hole symmetry is that two 
electrons at positions i and j feel the same repulsive force as two holes at the same positions, i.e. when the 
whole ID chain is full and only at i and j electrons are missing. 

Let us put this into mathematical terms. Let a| be a creation operator of a single-electron state with 
momentum ky ~ 2Trj/b, being therefore localised in x-direction around Xj = kyt^ (3). Assuming that j 
can take values 0, . . . , m — 1, particle-hole conjugated n-body states are 

. . . |0) (particles) < — > Uj^ . . . |1) (holes), (38) 

where |0) is an empty Landau level (vacuum) while |1) = . . . a|„_i|0) is a completely filled Landau 

level. For example aJajjO) = | • • • • • •) < — > | ••••••) = aoa2|l) = aiC^aalaslO). 

A straightforward calculation shows that matrices of a translationally invariant two-body operator A are 
the same (up to a multiple of identity matrix and complex conjugation) in an arbitrary n-particle basis and 
its conjugated (m — n)-electron basis. The only approximation we must concede is to neglect the Landau 
level mixing. 

Result of the calculation is the following. The diagonal terms of an operator A in the particle basis and in 
the hole basis fulfil 

(l|a]^ . . . al A a,„ . . . a,Jl) = A |1) + (0|a,, . . . A a]„ . . . a]^ |0) (39) 

and the off-diagonal terms remain the same up to the complex conjugation. 
Two cases are worth of special attention: 

The spectra of (fully polarized) systems at — n/m and v — [m — n) / m art the same up to a shift 

El = E\_^ + Ei[m-2n)/m, (40) 

where Ef is the energy of a completely filled (lowest) Landau level. This result does not depend on the 
form of the interaction V{r). A nice demonstration of this formula is shown in Fig. 41b (see the comment 
[3]). 

Conjugated states (38) may have different values of J (30). For instance: for m = 4, consider a three- 
electron state \j1i2i7i) ~ |013) and its particle-hole conjugate ~ \2). The former has J = while the 
latter has J = 2. 

For the density-density correlation function Q'sir) = (X]i<j ^i'' ^ ''i + ''j))* we get 

5*(r) = (l-cxp(-rV2^g)) + (41) 
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where ^ and are arbitrary particle-hole conjugated states. Note that g^i/ refers to electrons in the 'hole' 
state. Correlations between holes in ^P' are the same as those between electrons in 'f. 
Note, that g(r) in (41) is «of defined in the normalized form S{r — ri + rj)/{n{n~l)). Also g(r) of a full 
Landau level may depend on system (finite-size) parameters, e.g. in a rectangle with periodic boundary 
conditions, it depends on aspect ratio. 

Let us mention that densities of particle-hole conjugated states are related by nii,(r) = m — n^^/l{r), exactly 
as we expect from the picture of a hole as a missing particle. The plus sign in (41) might look puzzling. At 
the second glance, however, gijf = n • ri (schematically) and therefore g*/ = (1 — — n) = 1 — 2n + g*. 

1 .5.5 Density matrix renormalization group 

Exact diagonalization as it has just been presented, boasts of taking the complete basis of the lowest Landau 
level on a torus. As long as the low-energy states are considered, many of the basis states will be almost 
absent in the product-state expansion (32). Especially those which place many electrons close to each other 
and thus contribute with a large Coulomb energy. Leaving out such states from the basis will not affect the 
calculated ground state noticeably while it reduces the matrix sizes. 

Density matrix renormalization group (DMRG) is a systematic method to leave out irrelevant basis states. 
Roughly, its basic idea is to successively enlarge the considered system and to use only the most important 
n-particle states for calculating the (n + l)-particle ground state. 

The idea was used originally for one-dimensional systems (a review in [64, 65]). Shibata and Yoshioka 
[66, 67, 68, 69, 82] noticed that the single-electron basis of the lowest Landau level is in principle one- 
dimensional (23) and adapted this method as an extension of the exact diagonalization for studies of the 
lowest Landau level. They were thus able to study systems with up to about 20 particles at fillings close to 
V = 1/3. 

1 .6 Quantum Hall Ferromagnets 

Consider the situation u ^ 1 and vanishing Zeeman energy [26]. What is the ground state? 

In the absence of Zeeman splitting, the lowest Landau levels {n = 0) for spin up and for spin down have 

the same energy, thus, without interaction, there are 2eB/h single-electron states available with energy 

huj/2, which is the lowest energy an electron can have in the presence of a magnetic field B. Filling factor 

one means that only eB/h states (per unit area) are occupied. Hence there is a vast number of degenerate 

many-electron ground states without interaction. 

One of these states has the form 

vfH = <J'(;^i,...,z„)|TT •■•?>• 

Antisymmetry of implies antisymmetry of $, or in other words $ vanishes when any z,; approaches 
any zj. Each particle is surrounded by a correlation hole, cf. (45). Moreover, the state 5'// is the only 
one (with v = 1 within the lowest Landau level) whose orbital part is fully antisymmetric, up to SU (2) 
spin rotations. If we do not neglect the repulsive interaction between electrons, the 'optimal correlation 
hole' of the state h will make its Coulomb energy lower than for any other v ^ 1 state and h becomes 
the absolute ground state even at zero Zeeman energy. The long-range order in spins (all are pointing in 
the same direction) which are not localised at fixed positions, e.g. as it is in a spin lattice, renders h 
an itinerant ferromagnetic state. In the absence of the Zeeman spUtting the v = \ quantum Hall system 
constitutes an example of a Heisenberg ferromagnet. 

For Coulomb interaction, the energy cost of a single electron flip, which implies a violation of the anti- 
symmetry of $, can be evaluated analytically: E ~ (e^/e£o)-\/7r/8 [26]. Quantitatively, this number is 
comparable to the cyclotron energy hio at magnetic fields in the range of few tesla in GaAs. The fully 
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(a) The u = I QHF: (0, t) and (0, J.) levels are active, the rest 
is empty. A Heisenberg ferromagnet. The SU (2) spin sym- 
metry is manifest in the degeneracy of all possible Sz states 
for one given S. 



(b) The 1/ = 2 QHF: (0, T) and (1, J,) levels are active, (0, J,) 
is full and treated as inert, the rest is empty. An Ising fer- 
romagnet. The Z2 symmetry implies that only Sz and —Sz 
levels are degenerate; S is no good quantum number Ground 
state energy was shifted to zero. 



Fig. 7: Spectra of two examples of quantum Hall ferromagnets (exact diagonalization. Coulomb interaction, eight 
electrons). The ferromagnetic ground state occurs in both cases at /("" = which is a necessary condition for the 
state to be non-degenerate (in orbital degrees of freedom). See Fig. 8 for an explanation of the choice of Landau level 
indices and spins of involved levels (n, a). 



polarized state then becomes the ground state stabihzed by the huge gain in exchange energy. A spec- 
trum obtained by the exact diagonaUzation in a small system is shown in Fig. 7a. In agreement with the 
argumentation above, the ground state has S = n/2 and it is well separated from excited states. 
For the Pauli principle to apply ($ vanishes as Zi zj), it is only important that all spins have the 
same direction, not that they are all pointing upwards. Thus, the v = \ ground state is characterized by 
full spin polarization {S = n/2) and arbitrary Sz- All states {S^Y'^ h, k ~ 0, 1, . . .rt are degenerate 
ground states. A finite Zeeman energy will lift this degeneracy and the v ^ \ system will then have a 
nondegenerate ground state h, e-g- Sz = n/2 for B pointing in the z-direction. 

Other types of integer quantum Hall ferromagnets are possible, but they all share the common scheme: two 
degenerated Landau levels which provide 2eB/h 'free places' and only eB/h of them should be occupied. 
Depending on which two Landau levels are degenerate, different types of ferromagnets can follow. A 
classification of possible cases was given by Jungwirth and MacDonald [40]. 

Let us introduce one more example, the ly = 2 QHF which turns out to be an Ising type ferromagnet (see 
also Jungwirth et al. [41]). By changing the ratio between Zeeman and cyclotron energy, (n, a) = (0, "f ) 
and (1, 1) Landau levels can be brought to coincidence (Fig. 8). Experimentally, this can be accomplished 
either by changing the 51-factor (it decreases with pressure [15] or by tilting the magnetic field [18] (cy- 
clotron energy depends only on the component perpendicular to the 2DEG plane, Zeeman energy depends 
on the total field). The low lying (0, i) Landau level is fully occupied (eB/h states) and can be taken as 
inert. The remaining eB/h states (giving in total = 2) can be distributed among the 2eB/h available 
places of the two crossing Landau levels (Fig. 8). Contrary to the v = 1 QHF, there are only two ground 
states now; either (0, t) is full or (1, 1) is full (Fig. 7b). To obtain this result we should use the exact 
diagonalization because of the large degeneracy present when interaction is switched off. However, the 
fact that distributing the electrons between the (0, t) level and the (1,1) level costs extra energy (compared 
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Fig. 8: Integer quantum Hall ferromagnets occur when two 
crossing Landau levels should be only half-filled (2eB/h 
free states, eB/h electrons to occupy them). Thick (thin) 
grey levels indicate completely filled (empty) Landau lev- 
els, the pair of black levels are the 'active' ones. Despite 
the large number of possibilities of how to distribute elec- 
trons in the two levels, the Coulomb interaction selects (up 
to degeneracy in S,) one state. Depending on which Lan- 
dau level crossing is active, different types of ferromagnets 
occur: Heisenberg type for v — 1, Ising type for u = 2,3. 



to placing all electrons into one of the levels), is probably a consequence of the fact that spin up orbitals 
are not the same as spin down orbitals [40] as they lie in different Landau levels. 

For a more detailed discussion of spectra of a Heisenberg and an Ising QHF (Fig. 7), see in Subsec. 2.3.3. 
Quantum Hall ferromagnets which occur at integer filling factor have the advantage that they can often 
be well described by Hartree-Fock models, at least as far as the ground state is considered. Even here, 
exact diagonalization studies can sometimes unveil unexpected ground states, as shown by Nomura [58] in 
bilayer systems (spin degree of freedom is substituted by pseudospin which refers to the two layers). 
The principial question which is addressed in this thesis is, whether quantum Hall feiTomagnetism can also 
occur at fractional filling factors. Naively, one may expect that phenomena occuring for electrons (integer 
z/) would also occur for the CF (fractional v). The Coulomb energy (of CF cyclotron energy) would take 
over the role of the cyclotron energy within the integer QHF. A pleasing fact is that now the ratio of the 
CF cyclotron energy and the Zeeman energy is B-dependent, so that the coincidence of the CF LL can be 
induced just by adjusting B ata given filling factor Experimentally, there are strong hints on the existence 
of ferromagnetism even at fractional fillings [70, 22] and this work should contribute to the understanding 
of these phenomena from the side of theory. 

2 Structure of the incompressible states and of the half-polarized states 

2. 1 Basic characteristics of the incompressible ground states 

Being interested in phenomena occurring at the transition between two incompressible ground states, the 
spin-polarized and the singlet one, it is reasonable to get acquainted with these two ground states first. 
In the very illustrative model of non-interacting composite fermions (NICE), introduced in Subsect. 1 .4.2, 
the ground state at electronic filling factor i/ = 2/3 = 2/(2 -2 — 1) corresponds to two completely filled 
composite fermion Landau levels (LL). If, in some particular situation, the CF cyclotron energy is smaller 
than the Zeeman splitting, these will be the = 0, f, n = 1, | CF Landau levels and the ground state will 
be fully spin polarized. Fig. 9a. If the ratio between Zeeman and CF cyclotron energies is reversed, the 
ground state has n = 0,1, n — 0, 1 CF Landau levels filled and is therefore a spin singlet, cf. comment 
[4]. Here, the CEs are electrons with two flux quanta attached antipamllel to the effective magnetic field 
Beff [77], which leads to a minus sign in the denominator of the CF filling factor (19). 
A similar situation, i.e. occurrence of two incompressible ground states, the singlet and the polarized one, 
occurs also at filling factor v — 2/5. Here, the ground state can be interpreted as two filled CF LLs 
where the two flux quanta were attached parallel to B^ff. Thus, these ground states should be completely 
equivalent to the ground states at ;^ = 2/3 within the NICE approximation. 
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(a) Systems at filling factors v = 1/3, 2/3 and 2/5 correspond to vcF = 
1, 2 and 2 within the non-interacting CF picture. The composite fermions 
(CF) are electrons with two magnetic flux quanta attached parallel (for u = 
1/3, 2/5) or antiparallel (for u = 2/3) to the effective magnetic field B^ff 
(but always parallel to the real external field B). When Zeeman splitting is 
increased, crossings between CF Landau levels occur and spin polarization 
of the ground state changes. 



(b) Energy levels of 8 particles on a torus at 
filhng factors 2/3 and 2/5 without Zeeman 
splitting (short-range interaction, see Subsec. 
1.3). Note the large excitation energies (re- 
lated to the incompressibility gaps) for the 
ground states at 5 = and S = 4, as com- 
pared to other inter-level separations. 



Fig. 9: Ground states at filling factors 1/3, 2/5 and 2/3 from the point of view of a composite fermion picture and the 
exact diagonalization. 



Let us compare this picture of an infinite two-dimensional system with a finite system treated exactly. 
Looking at the exact spectra of a ly = 2/3 and a i/ = 2/5 finite system. Fig. 9b, we readily recognize 
ground states in the S* = and the 5* = N/2 sector which are well separated from excited states, as 
compared to the typical level separation within the excitation spectrum or in subspaces with other values 
of the total spin. Also, as the NICF model predicts, the spin singlet ground state (n = 0, t, = 0, |) has 
a lower energy £^(5* = 0) than the polarized one {n = 0,1^, n ^ 1, f), £'(5' = N/2) if the Zeeman energy 
is set to zero. Both ground states have k '' ^ (0, 0) which corresponds to L in a system with circular 
symmetry, Subsect. 1.5.2. Angular momentum equal to zero is in turn a property inevitable in any state 
with no partially filled Landau levels, corresponding argumentation is analogous to the comment [4]. 
It should be emphasised at this place that however strong support for the NICF model these findings pro- 
vide, they cannot be taken as a proof of its complete correctness. The interacting electrons cannot be exactly 
mapped to non- interacting CFs and even the quality of the approximation is hard to control. Although the 
NICF model gives correct answers to questions indicated above, there is no guarantee of correct answers 
in other cases, especially at other filling factors. In the following, we will continue discussing properties 
of both incompressible states at v = 2/3 and of those at ly = 2/5 as calculated by exact diagonalization 
and we will occasionally mention links to composite fermion theories. 

2. L 1 Densities and correlation functions 

Having computed a many-particle wavefunction numerically usually does not automatically mean that we 
can say much about the nature of the state it describes. Very often, the only statement to be made is that the 
state is highly correlated, or entangled. By this we mean that the state cannot be written as a single Slater 
determinant [5], not even approximately, and thus its description goes far beyond any Hartree-Fock model. 
In order to learn more about the state it is apt to evaluate expectation values of observables such as den- 
sity or density-density correlation functions. In the first quantization formalism these are the following 



28 



K. Vybomy: Spin in fractional quantum Hall systems 



operators 

n{r) = 5: <5(r - r,) , g{r) ^ ^ ^ ^(^ " ^^)) ' ^42) 

summations running over all particles in the system. For inhomogeneous systems it is also useful to 
consider an unaveraged density-density correlation operator 



-'^'"•'■)-mh^T.oi''-'-)»i'"-'A 



where an average over r" gives (42). This is the probability density of finding a particle at place r' 
provided there is a particle at place r". The function g{r) is just g{r + r', r') averaged over all r', hence 
g{r) (X g{r + r', r') for homogeneous systems, i.e. both quantities are the same up to a proportionality 
constant. 

For not fully spin polarized states it is also useful to watch quantities 7i|(r) or g^^l{r) and its analogues 
with other spin indices. For example 

The normalization of density and density-density correlation functions we chose in (42,43) is the following: 

'drn{r)=N,, jdrg[r) = l, y drg,,(r) = ^|^^ , a e {t, i} , (44) 
where integrals are taken over the whole system i.e. elementary cell. 

As long as homogeneous systems are concerned we naturally expect density and also polarization to remain 
constant. For the incompressible states this is true only up to finite size effects. The density shows a slight 
modulation which decays rapidly as the system size is increased. Discussion of these effects which have 
no relevance for the real infinite 2D system will be presented later, Subsec. 2.1.4. 

In the following, by g{r) wemean5(r) with r — |r| for isotropic and homogeneous systems. Also, when- 
ever we will speak about 'correlation functions' we mean equal time density-density correlation functions. 

Fully occupied Landau levels 

The density-density correlation function can be analytically evaluated for a state with v — n fully occupied 
Landau levels [42]. This is the ground state of non-interacting electrons at integer filling factor. For the 
spin polarized case. 



g{r) = 1 - i-cxp(-[(rfcF)VH) 



ri 



{rk 



F) 



^2^ -12 



An 

Here L'^{x) are the associated Laguerre polynomials [6, 29]. In particular. 



(45) 



1^ = 1 ■■ 9.=i{r) = l-cxp(-rV2£g) , y = 2 : g^^^ir) = l-exp(-rV4^g)-^[2-rV4€^]' . 

(46) 

The Fermi wavevector kp for a system subjected to a perpendicular magnetic field is defined as kp in 
exactly the same system (i.e. the same areal density of electrons) just with magnetic field switched off. In 
this scheme 

{kpiof = , or kp = V2^ig^ . (47) 
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Correlation function, integer filling factors Correlation function, free 2D Fermi gas (B=0) 




Fig. 10: Density-density correlation in a free 2D electron gas and in magnetic field at integer filling factors (spin- 
polarized electrons). 



It is a pleasant news that by taking the limit v = n 

n n 2 



oo in (45) we obtain 



9Fs{r) = 1 - 



-; Jiikpr) 



(48) 



which is the correlation function of free electrons in two dimensions (Fermi sea). It should not be anything 
else because ly ^ oo with kp kept constant means that B is decreased to zero at a given areal density of 
electrons. 



Filling factor 1/^1/3 

Provided Landau level mixing is absent and considering only the short-range interaction between parti- 
cles (Sec. 1.3), the ground state at filling factor ly = 1/3 is described by the Laughlin wavefunction 
. . . , z„), (7). Up to my knowledge, no closed [7] analytical expression of the correlation function 
in this state is available. Only the short range behaviour can be determined analytically. For (zi — 22) ^ 0, 
I'I'lP vanishes proportional to {z^ — 22)^ {zi — 22)^ = \zi ~ Z2\^, hence g{r) = cr^ + o(r^) for r ^ 0. 
Numerically, {'^ L\9{r)\'^ l) can be evaluated by various Monte Carlo techniques. Fig. 11. These results 
are closer to the thermodynamic limit, referring to larger numbers of particles, than g{r) which can be 
obtained from exact diagonalization. Fig. 12. This is however only because we know an analytic WF 
of the GS for arbitrarily large systems in this case, ^P^. Exact diagonalization can be performed only 
for systems with < 10 electrons, but it is not necessary to know anything about the ground state in 
advance apart of that it Ues in the lowest Landau level. Therefore, exact diagonalization provides us a way 
to confirm that ^'l is indeed the ground state or a good approximation to it, e.g. for Coulomb-interacting 
electrons. Note also that Figs. 12 refer to electrons on torus whereas Fig. 11 refers to the disc geometry. 
Indeed, the correlation functions are very similar in both geometries, compare Fig. 12(b) and Fig. 11. This 
fact supports the hypothesis that the corresponding states, ^fi, on a disc and those on tori, are universal and 
hence basically the same as the ground state in an infinite 2D system. 

The correlation function g(r) in Fig. 12 is rather isotropic, at least on distances smaller than a/2. This 
distinguishes the Laughlin state from a Wigner crystal (Subsect. 2.4.1) or a unidirectional charge density 
wave in which some special directions exist, Subsect. 2.1.3. This fact motivates also the incompressible 
liquid terminology. Regarding the distinction between liquids and gases, the difference is the strength 
of interparticle interaction. Whereas negligible in gases, the interaction in liquids is strong compared to 
kinetic energy. In the lowest Landau level, kinetic energy is zero, or better a constant huj/2. 

The first maximum in g{r) occurs at ri w 4.4fo, Fig. 12b, and this separation can be taken as a typical 
interparticle distance in the Laughlin state. This distance lies close to the mean interparticle distance 
determined by the filling factor, rmean/^o = \/2Ti jv k, 4.35 (1). 
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Fig. 1 1 : Left and middle: Correlation functions of the ground states of 50-60 particles at filling factors 1 /3 and 
2/5 of the principal Jain's sequence, v — p/{2p + 1) (cf. Subsect. 1.4.2). The wavefunctions (WF) predicted by 
composite fermion theory were taken (for u — 1/3 this is identical with the Laughlin WF) and g{r) was calculated by 
a Monte Carlo method. Taken from Ref. [42]. Right: correlation function between oxygen atoms in liquid water as an 
example of a density-density correlation function in a well-known liquid (see text on p. 30). Results of both numerical 
simulation and experiments are shown, see the original paper by AUesch et al. [9] for details. 

After ri, oscillations in g{r) decay rapidly. The overall form of g{r) in the Laughlin state clearly differs 
from the correlation function of a free 2D Fermi gas (48). We will emphasise three aspects. 

(i) Laughlin state, Fig. 12b: the first peak of g(r) is relatively high, measured for instance by ratio 
g{ri)/g{a/ V2) > 1.1. Here a/-\/2 « lO^o is the maximum interparticle distance in the considered fi- 
nite system. 

2D Fermi gas. Fig. 10: the first structure of g{r) is about ten times weaker. Here, it is more appropriate to 
watch the depth of the first minimum, see (ii). 

(ii) 2D Fermi gas: all maxima (at Vpg) of g{r) have the same value, g{rpg) — 1. Laughlin state: the first 
maximum g{ri) w 1.1 (for oo) is much higher than other maxima. 

(iii) 2D Fermi gas: g(r) cx for r ^ 0. This is purely the effect of Pauli exclusion principle. Mathemat- 
ically, it comes from the antisymmetry of the wavefunction ^P, in other words, is a Slater determinant. 
Laughlin state: g{r) (x r^. This is a manifestation of correlations in the state, i.e. of the fact that 
cannot be written as a single Slater determinant. g(r) cx 7'^ also means that any two electrons avoid being 
close to each other very efficiently and this helps to minimize the Coulomb energy which is high at short 
inter-particle distances [33], Subsect. 1.3.3. 

Just as an illustration, a correlation function g{r) of liquid water is shown in Fig. 11, right. Of course, 
it is not possible to directly compare water and a 2D electron gas in the fractional quantum Hall regime. 
Nevertheless, the pronounced structures in g{T) beyond the correlation hole in the Laughlin state. Fig. 1 1 
left, are definitely more similar to g{r) of liquid water. Fig. 11, right, rather than to g{r) of a 2D Fermi gas, 
Fig. 10. 

The Laughlin state. Fig. 12b, also differs from integer filling factor states apparently. Fig. 10. The latter 
ones (i = 1,2,...) namely have always g,y~i{r) oc at r ^ 0. Also g,y=i{r) has exactly i — 1 maxima, 
i.e. (7^=1 (r) is free of maxima. 

This demonstrates the fact, that in the vcf — ^ composite fermion (CF) state, which is the model of the 
= 1/3 electronic ground state (Sect. 1 .4), the electron-electron correlations are different to those in a 
V = 1 electronic state. This is a bit counterintuitive, since the CFs were created by adding two zeroes 
to electrons in the v ~ \ state and we could have therefore expected that the electrons 'remained at their 
original positions' under this transformation. Figures 12b and 10 however show that even though the CF 
density equals the electronic one the electron-electron correlations are different in both states. 
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(a) Ne = 10. Due to the absence of circular 
symmetry on the torus, (;(r) is in general not 
only a function of r = | r | . For \ r\ <^ a, 
g(r) is however quite isotropic, Fig. 12b in- 
set. The first electron is sitting at the cor- 
ner, the four corners are identical owing to 
the periodicity. The lower plot differs from 
the upper one only by a finer z-scale which 
highlights the structures in (;(r) at larger dis- 
tances. 



(b) Section of g{r) for Ne = 5,8,9, 10-electron ground states along 
r = (x, x); the perfectness of match to (j(r) in Fig. 1 1 gives us a feeling 
how little the ground state is affected by the finiteness of the system. It 
is noteworthy that g(r) can be astonishingly well fitted by the [gFs('')]^ 
(48) up to distances beyond the first maximum (up to vertical scaling, only 
kp must be fitted, see the text). Inset: sections along diagonal and side 
of the square, i.e. g{x/V^,x/V^) and g{x, 0) for the A^e = 10 system. 
The good match of the two curves as far as well beyond the first maximum 
(~ 6^0 ) indicates that the isotropy on length scales < 6^o is not much 
affected by the rectangular geometry (periodic boundary conditions). 



Fig. 12: Correlation functions in the ground state of A^e electrons on a torus (square, length of sides a = 6, with 
periodic boundary conditions) at filling factor v = 1/3. The function g{r) gives the probability of finding an electron 
at position r = {x, y) provided there is an electron sitting at r' — (0, 0). For a homogeneous system the choice of r' 
does not influence the probability distribution of finding the second electron. 



On 'intermediate length scales' (1 to 5 magnetic lengths), the correlation function of the Laughlin state 
g{r) in Fig. 12b can be strikingly well fitted by 

c-[9Fsir)]\ (49) 

where gpsii") is the correlation function of a free 2D Fermi gas, (48). Herefore, we put kp w 0.874£q ^ 
which is only by about 7% more than what we would expect for filling factor = 1/3, (47). 
The quality of the match relies on the choice of m = 3 for the exponent in Expr. 49 (for r ^ 0) and 
on the fitting constants c and kp (around r « ri). The surprising fact is therefore only the good match 
between r = and r = ri. Also note that long-range (r ^ ri) behaviour of expression (49) and of g(r) of 
the Laughlin state are different. This again emphasises the differences between the Laughlin state and the 
Fermi gas. Expression (49) provides therefore only another representation of the exchange hole, parallel 
to approximate formulae given e.g. by Girvin [25]. 

In conclusion, we have seen that the correlation function of the correlated = 1/3 ground state (Fig. 
11) has a strong first maximum (near to 4.4^o) and an unusual exchange hole g{r) cx r^. These features 
distinguish the 1 /3 state from both free 2D Fermi gas and completely filled Landau levels and indicate the 
liquid-like and correlated nature of the Laughlin state. 
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(a) The polarized state. 
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(b) The spin-singlet state: correlation between unlike spins (left) and like 
spins (right). 



Fig. 13: The u — 2/3 ground states shown in their correlation functions. Section along the diagonal of the square 
elementary cell is shown, i.e. g{r) — g{x/\/2,x/\/2) and systems of different sizes are compared. 
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(a) The polarized state. 
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(b) The spin-singlet state: congelation between unlike spins (left) and like spins (right). 



Fig. 14: The f = 2/5 ground states, correlation functions g{r) = g{x/V2,x/V2). Note the good match of peak 
positions in systems of different sizes. 



Filling factor v ~ 2/3 

Provided the Landau level mixing is absent, the particle-hole symmetry in one Landau level gives a direct 
relation (isomorphism) between Hilbert subspaces of fully polarized states at;^ = 2/3 = 1 — 1/3 and 
1/ = 1/3, Subsect. 1.5.4. Owing to this relation eigenvectors of any radial two-particle interaction are 
exactly the same' in both spaces and corresponding eigenvalues are identical up to a constant shift. 
The correlation function in the/M//y polarized i/ ^ 2/3 ground state, Fig. 13a, is thus linked to the one of 
the Laughlin WF by an analytical formula (41). For a system with Nm flux quanta, i.e. having an area of 
27r^giV„i, it reads 



2/3iV™(2/3/V„ - l)5.=2/3(0 = l/3iV,„(l/3/V,„ - l).gi/3(r) + l/3N^,,g,=i{r) . 
The g{r) oc short range behaviour is thus obscured by the second term. 



(50) 



In the following sense: Take an eigenvector for u = 1/3. This is a linear combination of Slater determinants from the u = 1/3 
space. Replace each of them by its pailicle-hole counterpart and the resulting state from the u = 2/3 space is an eigenstate. 
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Fig. 15: Different density-density corre- 
lation functions in the v = 2/3 singlet 
state seem to be related to each other: 
5Ti.('') +5tt('') (properly normalized, see 
text) is very similar to the function 1 — 
exp(— r^/4^o)' the density-density corre- 
lation in a full Landau level with ^o\/2 in 
the place of la. 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 
x/lo 

The spin singlet ground state at i' = 2/3 has a different character. Here, we can distinguish between 
correlation functions for electrons of like spin, (7||(r), and for electrons of opposite spin, g-\i{r). Fig. 13. 
Neither of them bears any apparent resemblance to either the v — 2/3 or i/ ^ 1/3 polarized ground states. 
We should like to point out some of their particular features. 

(i) The ring-like form of g^i{r) suggests that the state consists of pairs of particles with opposite spin with 
average separation r^i sa 3.3^o- 

(ii) There is a deep hole in g^i{r) around zero. This cannot be due to Pauli exclusion principle which 
applies only to electrons of like spin, but rather solely due to Coulomb repulsion. As a check (not presented 
here), a comparison between g-iiir) in Fig. 13 and the 'lowest LL Pauli hole' g^=i{r) (46) reveals that 
their forms are indeed different. Also note that the value of g-^i (0) is not exactly zero, it is several percent 
of the maximal value of guir), Subsect. 2.1.2. 

(iii) There is a well pronounced shoulder in around r w 2£o. It is very suggestive, how well this 
shoulder can be fitted by the correlation function of a full lowest LL, g^=i{r), i.e. the lowest LL exchange 
hole (46). This is shown in Fig. 16a. 

This feature reminds of the relation between 1/3 and 1 — 1/3 systems (50). This is also supported by 
the fact, that after the shoulder is subtracted [g,^=i{r) times a constant], the remaining part of 5||(r) is 
g{r) oc at short distances (Fig. 16a), just as it is the case in the 1/3 Laughlin state. However, particle- 
hole conjugation between filling factors 1/3 and 2/3 is applicable only for spin-polarized states. 

(iv) The sum of g||(r) and g^iir) properly scaled for A^e — > oo lies very close to g^=i{r) with (q sub- 
stituted by ioV^, Fig. 15. Proper scaling means that g||(r) and g^i{r) should have the same norm, e.g. 
equal to one, in sense of (44). With the current notation (43) this is true only for oo. 
Therefore, if spin is disregarded, the singlet ground state at i/ = 2/3, created by magnetic field B, strongly 
resembles the state of a completely filled lowest LL at magnetic field B/2. 

Summary: the polarized ground state at = 2/3 is the particle-hole conjugate of the Laughlin state at 
1/ = 1/3. The electronic coiTelation function of the 2/3 state reproduces the liquid-like maximum at 
ri sa 4.4^0 but the = 1/3 broad exchange hole with g{r) oc is hidden behind the lowest LL exchange 
hole, gu=i{r). 

The singlet GS seems to consist of pairs of spin up and spin down electrons with characteristic size of 
3.3£q. Together with the sum rule, point (iv) above, this could be interpreted as that Ne electrons in the 
singlet GS form Ne/2 pairs, each with total Sz ~ and these pairs form the same state as Ne/2 fermions 
at = 1 in the ground state. 

In particular, it should be emphasised that the singlet state cannof be described as a mixture of two mutually 
uncoiTelated = 1/3 Laughlin liquids, one with spin up, another with spin down, as we could wrongly 
infer from the picture of non-interacting composite fermions, see comment [8]. 




Filling factor ly ^ 2/5 
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This filling factor should be the counterpart to = 2/3 within the CF picture. The two magnetic fluxes are 
attached parallel rather than antiparallel to the effective magnetic field and in both cases the CF filling is 
two (Sect. 1.4). In spite of this relation the density-density correlations between electrons show significant 
differences. 

The correlation hole of the polarized ground state (Fig. 14a or Fig. 11, middle) is much broader for 
V = 2/5. The first maximum occurs in both systems (2/5 and 2/3) at about the same distance « 4.1£o, it 
is however much better pronounced in the 2/5 system and also more structure is present beyond the first 
maximum here. Around r = both systems follow g{r) oc r^. However, whereas g{r) for v = 2/3 
is dominated by the 'exchange hole', i.e. g,y=i{r), cf. (50), the 2/5 state has a much broader minimum 
around r = 0. 

These findings are not unexpected. Consider two systems of the same area 2n(.QNm, one at fillings 2/3 
and 2/5, respectively. The latter will be more diluted ('emptier'), since it contains only 2/5N„i electrons, 
compared to 2/3N„i in the v ~ 2/3 system (1). Therefore, the correlation hole in g(r) can be broader in 
the 2/5 system. This conclusion is not a controversy of the CF picture, rather, it is a warning. By far not 
all claims which are true for electronic Landau levels (e.g. ly = 2 state remains the same regardless of the 
direction of the magnetic field) are true for composite fermions as well (it matters whether B^jf is pointing 
parallel and antiparallel to the attached flux quanta). 

There is also a close relation between the polarized 2/5 GS and the Laughlin 1/3 state according to the CF 
picture. The latter one corresponds to filling factor one, the former one to filling factor two of composite 
fermions. Comparing these two states, we find a bit stronger structures in the density-density correlation 
of the ly = 2/5 GS and also the first maximum shifts to smaller distances {AA£q at i/ = 1/3 and 4.1^o at 
1/ = 2/5). Both effects are quite similar to what happens when going from ly = 1 to i/ = 2, cf. Figure 
next to (46). Comparing the jy ^ 2 and 2/5 systems, we again (cf. = 1 and 1/3) find much stronger 
structures of g{r) in the latter case, just as we expect for a liquid state. 

Some marked differences occur also in the singlet ground states at both filling factors. At 2/5, correlation 
functions g^^{r) as well as g^i{r) seem to be quite flat beyond r„i w 6£o. We may speculate that the 
same is true for the filling 2/3, Fig. 14 (the shoulder in g^^ would probably have to be subtracted first), 
but then the plateau would occur first beyond some larger distance which is not accessible by exact 
diagonalization. 

It is remarkable that after subtracting the shoulder from g^^{r) of the 2/3 singlet state (point (iii) in the 
discussion of 2/3), the rest g{r) is oc near to r = 0. This is the same behaviour as we find in g^i{r) of 
the 2/5 singlet state, see Fig. 16. 

Correlations of unlike spins exhibit one clear maximum which is, as compared to 2/3, slightly but percep- 
tibly shifted to a bit larger r-^i sa 3.7£o. This agrees with the above argument that 2/5 systems are more 
diluted than the 2/3 ones, but quantitatively this shift is too small. It is only « 30% of what we would 
naively expect from comparing the areal electron densities. 

Finally, the r ^ behaviour of the 2/5 singlet state, g^i{r) oc and g^^{r) oc r^, matches the behaviour 
of the {3,3, 2}-Halperin wavefunction (8) and this $332 [z] is in turn identical with the ground state wave- 
function proposed by Jain's theory, Subsect. 1.4.2. This is because $„„'m[^] lies completely in the lowest 
LL and thus the last step of Jain's procedure, namely the projection to the LLL, is out of effect. Seen from 
the opposite direction: the singlet 2/5 state corresponds to filling only the lowest CF LL spin up and spin 
down. 

Summary: From the viewpoint of composite fermion theories, the polarized 2/5 state {p = 2), Tab. 2, is 
related both to the 1/3 Laughlin state {p = 1) and 2/3 polarized ground state {p = —2). The electron- 
electron correlations in exactly diagonalized systems clearly support the former relation, the latter one (2/5 
with 2 /3) is however far from being obvious in this way. 
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Short-range behaviour of g^-|-(r), 2/3, S=0 

N^/N^ =10/15, S=0 ground state, section g(x,x) 
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(a) Filling factor 2/3. The shoulder at r 2£q is apparently 
caused by a term proportional to 1 — exp(— r^/2^g), i.e. 
g^—i[r), cf. (46), which contributes to the total g-jfir). Af- 
ter this term was subtracted, a local power analysis of gif{r) 
has been performed. 



Short-range behaviour of g-|"f^(r), 2/5, S=0 

N /N =8/20, S=0 ground state, section g(x,0) 




r/a 



(b) FiUing factor 2/5. Local power analysis near to r = 0. 
Noise at very small distances is purely due to numerical inac- 
curacies: values oi g{r) are already very small there. 



Fig. 16: Correlation of like spins, gT\{r), of the singlet ground states at filling factors 2/3 and 2/5. Local power 
analysis (51) shows, that both correlation functions are oc for r 0; however, the shoulder in the state at filling 
factor 2/3 has to be subtracted first. 



Neither is the analogy between 2/5 and 2/3 apparent for the singlet ground state. Although similarities ex- 
ist, perhaps most importantly pairing between electrons of unlike spin, short range behaviour of correlation 
functions is very different. 



2.1.2 Ground state for Coulomb interaction and for a short-range interaction 

Short-range interactions as they were introduced in Section 1.3 have a special significance for the FQHE. 
It has been repeatedly emphasised that the Laughlin WF is on one hand an extremely good approximation 
of the ground state of a Coulomb-interacting (CI) system while on the other hand, it is the exact ground 
state of electrons feeling only a short-range mutual interaction (SRI) as it was defined in section 1.3. 
Consequently, it is very popular to say that a short-range interaction Hamiltonian captures the essential 
physics of the FQHE by inducing the correct correlations in the ground state. By the correct correlations 
we mean the 5* cx (z^ — ZjY behaviour when two particles approach each other. 

The SRI was used in most of the calculations presented in this work. This choice has been made for two 
reasons. It brings better chances in finding analytical results like the Laughlin WF. Moreover we may 
hope that the results in finite systems converge faster to the thermodynamical limit {N oo) because the 
electrons see only as far as their interaction reaches and thus - sooner than for a long-range interaction - 
they will not 'realize' anymore that they live on a torus and not in an infinite plane. Aim of the following 
section is to show and discuss how the ground states at = 2/3 change if the character of the interaction 
changes. 

The ground state energies for CI and SRI are naturally quite different. This is however for the largest part 
only an unessential shift, a part of it is the missing Madelung constant (34). Under SRI an electron of 
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course cannot interact with its own image in the neighbouring primitive cell. More importantly, the gap 
energies are quite similar in both cases, Subsect. 3.2.1. 

Since the density of the incompressible ground states should always be constant, up to finite size effects 
to be discussed later, let us now focus on correlation functions. The three plots in Fig. 18 show (7||(r) 
and g^i{r) of the singlet state and g(r) of the polarized state. In all three cases, the correlation functions 
of the CI state and the SRI state are quite similar. Most apparent differences appear at large distances. 
On a torus, the largest possible separation between two electrons is r = a/\/2. On the other hand, the 
correlation functions are very precisely identical for small r. This shows that, e.g. in the polarized GS, the 
wavefunction contains the factor {zi — Zj)^ for CI as well as for SRI. In other words, the Laughlin state (as 
the GS for SRI) describes exactly the short-range behaviour of an incompressible state of even long-range 
interacting electrons. Fig. 18 demonstrates that this is true (at least in a very good approximation) also for 
other ground states where the analytical wavefunction is not available (e.g. the singlet GS). 
In fact, for the singlet GS there is a tiny but perceptible difference in g-^i (0) for the two types of interaction. 
Since g^i (0) is almost zero, this observation suggests that a yet modified interaction might lead to analytical 
results, {Vb, Vi, . . .} = {oo, a, 0, 0, . . .} in terms of pseudopotentials. Sect. 1.3. Such an interaction 
enforces g^i{Q) = 0, which is anyway almost fulfilled for the current SRI, and on the other hand it retains 
the pleasant property of SRI in polarized systems, i.e. it is one-parametric. 

There is yet another significant difference between SRI and CI which is not obvious in Fig. 18 at first 
glance. The difference concerns the placement of zeroes in the wavefunction and we will concentrate on 
the v = 1/3 ground state now (see Sec. 1.3). 

In a general fermionic state, there must always be a zero bound to each electron in order to fulfil the Pauli 
exclusion principle: two electrons (of the same spin) cannot be at the same point in space simultaneously, 
ergo if zi = Zi then the wavefunction must vanish. Factors {zi — Zj)^ in the Laughlin state mean that there 
are two extra zeroes exactly at the position of each electron. That is why g{r) oc for small r's. Fig. 17, 
right and (51). For CI, the Laughlin WF is only an approximation to the ground state. In the real ground 
state, the one obligatory zero is still sitting on each electron and the two others are only near rather than 
exactly on the top of the electron. In Fig. 17 we can even see how far they are on average. These two 
extra zeroes are now mobile and their position depends on the position of all other electrons. Note that this 
distance depends on the system size [53]. 



Local power analysis 

A comment is due on the way how the plots in Figs. 17,16 were obtained. It is basically a section of g(r) 
along one straight line going through r = 0. This function was then transformed by 

9{r) ^ (51) 

which gives a local degree of the polynomial behaviour. Let us give an two examples. If g{r) were ar" 
then d ln3(r)/dlnr — n. If g(r) oc (r — tq)" then dlng{r)/d\nr — nr/(r — j-q) n for r ^ tq. In 
other words, if there is a dominant r" term in g{r), the quantity plotted in Fig. 12(b) gives the exponent. Of 
course, it is only approximate except for the case g{r) = ar" but it is quite easy to evaluate and moreover 
it gives a global property of the wavefunction as compared to fixing electron positions Z2, . . . , z„ and 
examining the WF as a function of zi where results depend on where we fix the electrons Z2, . . . , z„. 



2.1.3 Some excited states 

There is a rich variety of excitations to the incompressible FQH states. For instance quasiholes, excitons 
(quasihole-quasielectron pairs), charge density waves (CDW) or spin density waves (SDW), all of them 
can be described analytically (at least to some extent), and then of course all the rest of excitations which 
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Fig. 17: The incompressible ground state alv — 1/3 (with ten electrons), Coulomb interaction (left) and a short-range 
interaction (right). Section through the density-density correlation function g{r) along r — {x,x) is taken and the 
'local degree' of the polynomial behaviour is determined (see the text). While the local behaviour around r = is 
g{r) oc for the SRI, indicating that there is exactly a triple (6 = 2-3) zero of the wavefunction on on each electron, 
we can clearly see only one zero at each electron's position for the Coulomb interaction, g{r) cc r'^ and 2 = 21. 
However, going away from r = 0, the 'local degree' grows and beyond ~ 1.5i!o it approaches the curve of the SRI 
state. The conclusion is that one zero (the obligatory Pauli exclusion principle zero) is fixed to each electron (r = 0) 
in the Coulomb state and the other two zeroes are only loosely bound to the electron. First from distances > l.Si'o this 
compound object looks like an electron with two attached flux quanta. 

has not been understood up to now. Following the introduction given around (28), we will now demonstrate 
how to identify some of these excitations in spectra obtained by exact diagonalization at the example of 
V = 1/3. 

Charge density waves 

CDWs can be excited for example in the liquid GS at = 1/3. Disregarding the possibility of spin flips 
(as it may be reasonable when Zeeman energy is too high), it turns out that these are the lowest excitations. 
In Fig. 19 spectra of several short-range-interacting = 1/3 systems (tori of different sizes) are presented. 
The horizontal axis is modulus of k^, i.e. the 'crystallographic fc-vector' described in Subsec. 1.5.2. The 
Laughlin state has = and a CDW of wavevector Q excited from this state has /('' = Q. Beware 
however, that not every state which has /c*" 7^ must be a charge density wave! Apart from other possi- 
ble periodic excitations, there are also basically nonperiodic excitations (e.g. quasiholes) and such states 
are forced into periodicity only 'artificially' by the periodic boundary conditions imposed in our exact 
diagonahzation model. 

The lowest excitations in Fig. 19 form a well developed branch E{k^), which is usually called magnetoro- 
ton branch, and other excited states form a quasicontinuum. The dispersion of the magnetoroton branch 
can be calculated analytically in the single mode approximation. The original calculation by Girvin et al. 
[27] for Coulomb interacting systems ?A.v = 1/3 showed a well pronounced minimum in £'(|/( '"|) of the 
magnetoroton branch at k^ig w 1.4. In a short-range interacting system, shown in Fig. 19, the situation is 
slightly different. Having reached its minimum value, E{\k^'\) remains constant beyond k^io w 1.4. 
A point worth of emphasis is that the magnetoroton branch in Fig. 19 contains points (energies) from 
exactly diagonalized systems of different sizes. This confirms our hope that these states are not bound to 
some particular geometry of the elementary cell and that they appear also in an infinite system. 
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Fig. 18: Correlation functions of the singlet and polarized ground states at — 2/3: comparison between the Coulomb 
and short-range interaction. The curves are identical for small r and slight deviations occur at longer scales. This is 
another way to demonstrate that it is sufficient to consider short-range interaction in order to get (almost) correct 
ground states under FQHE conditions. 



Dealing with finite systems, we will always have only a finite, and usually quite small, number of allowed 
values for (27). On the other hand, the more points in /('-space we can access, the better we can 
recognise modes in exact diagonalization spectra, just like the magnetoroton branch in Fig. 19. Note also 
the large space between /c'' = and the next smallest j/c*"! « 0.5£q ^ in Fig. 19 which corresponds to the 
longest wavelength compatible with the periodic boundary conditions. 

The traditional way to improve these limits (few /(''-points, too large smallest \ k^\ > 0) is to study larger 
systems. This is however prohibitively difficult with exact diagonalization. An alternative approach may 
be to study systems with aspect ratios X ~ a : b slightly deviating from one. This allows us to deform 
the lattice of allowed /(''-points continuously (27 contains A), and on the other hand, we can expect that 
the states will not suffer from the slight asymmetry in a : 6 in line with the argument that these states 
are not bound to any particular geometry of the elementary cell. This method is demonstrated in Fig. 
19 by the blue points. The aspect ratio was varied from one up to 1.3. Since the energies of the CDW 
states still lie well on the magnetoroton branch, we can conclude that this variation was still only a small 
perturbation, i.e. acceptable for studying this branch. A more reliable critierion would be to check overlaps 
of wavefunctions at a : 6 = 1 and a : b > 1. 

CoiTelation functions of several states in the magnetoroton branch (Fig. 19) are shown in Fig. 20. The first 
look at g{r) (upper row in Fig. 20) may be sometimes not enough to distinguish their charge density wave 
nature. The CDW is superimposed on the structure of the mother Laughlin state, which these states are an 
excitation of. The periodic structure of g{r) is thus more clear if we subtract the corresponding correlation 
function of the LaughUn state first, Fig. 20 lower row. We can find three periods in y direction (horizontal 
waves at a; = 0, 0.3 and 0.6) in the state A or 4 periods in y and one period in x in the state C, in agreement 
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Fig. 19: The ground state (at fc = 0) and low excita- 
tions in SRI fully polarized v = 1/3 systems of differ- 
ent sizes (4-10 electrons). Energy plotted against the 
\k'^\, see Sec. 1 .5.2. To be able to study more points in 
k space, systems of different sizes are compared and 
also systems with aspect ratios slightly varying from 
one, see text. Note the well pronounced magnetoro- 
ton branch. Note that energies calculated in systems 
of different sizes lie on the same branch indicating that 
these states are not much system-size-dependent (and 
therefore relevant even in infinite systems). The cor- 
relation functions for three states lying on this branch 
(A,B,C) are depicted in Fig. 20. 



Laughlin 



state A, k7ko=(0,3) 



state B, k7ko=(0,4) 



state C, kVk(,=(1,3) 







x/a 




Fig. 20: Correlation functions of the ground state (Laughlin) and several CDW states which lie on the magnetoroton 
branch of a v = 1/3 system with short-range interaction (eight electrons). Upper row: correlation functions g{r), 
lower row: g{r) of the CDW states from which the Laughlin state g{r) has been subtracted. 



with their values of k^' /kg = k''' / {Nmi^ /Q). Note, that it is harder to distinguish the periodic structure in 
the /c*" = (0, tt) state (B), which may be partly because this is a point of high symmetry in /c'"-space. Fig. 6. 
In conclusion, we have shown how (the best known type of) charge density wave states on a torus can be 
identified in the exact diagonalization spectra and in correlation functions. Generally, we can expect that 
charge density waves excited from incompressible liquid states will form branches in provided 
of course that their energy is not hidden in a quasicontinuum of other excited states. Correlation functions 
show indeed the expected periodicity of a CDW superimposed on the structure of the ground state. 



2.1.4 Finite size effects 

Consider a = 1/3 system with its exact GS written as 4*/,, the LaughUn wavefunction (WF), see (7). 
Particle density in the state ^'l is very precisely constant provided we stay within the disc of radius 
^oV27r • 3A^. The first striking observation is that the density of the ground state obtained from exact 
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(a) Inhomogeneous density in an A^e 
4 system. 



(b) Comparison between systems of different sizes. 
Note the scaling of the x axis. 



Fig. 21: Density of the v = 1/3 incompressible ground state as obtained in torus geometry for different system sizes. 
The oscillations can be traced back to the centre-of-mass part of the wavefunction. As far as this effect is considered, 
differences between Coulomb interaction and short-range interaction are small. 

diagonalization varies quite strongly. Fig. 21. At the same time we notice that the ground state, which is 
claimed to be incompressible, hence non-degenerate, is actually triply degenerate. 

Fortunately, this does not mean that finite size calculations are completely wrong. Both facts can be 
attributed to the centre-of-mass part of the wavefunction (CMWF) which is not present in but is present 
in numerical calculations, Subsec. 1.5.2. As far as isotropic states are considered, this seems to be the most 
serious effect coming from the finite-sizedness of the system and in the following we will discuss its origin 
and how it can be eliminated. 



Centre-of-mass part of the wavefunction 

The complete WF of the Laughlin state at ly 
1.5.1) might be 



1/3 (for n particles) in the disc geometry (Subsec. 1.2, 



^-1/3(^1, . . . , z„) = F(Z)cxp(-|Z|V2^^) X exp ( - (jzip + . . . + \zX)/Ml) J|(z, - (52) 



i<j 



<I'i,(zi,...,z„) 



with Z ~ zi + . . . + Zn and for example F{Z) = or any other analytic function with three zeroes 
Zi, Z2, Z3. The CMWF '^cm has the form (Sec. 1.5.1) of a WF for one particle somewhere in the lowest 
Landau level to which a single variable Z is attributed. In torus geometry, the WF must be changed in order 
to comply with periodic boundary conditions (PBC) which amounts to replacing {Z — Zi) terms by theta 
functions of the same argument, Subsec. 1 .5.2. Example of "^cm obtained from the numerically calculated 
ground state '^gs in ^ system with four particles is shown in Fig. 22a. The CM part was extracted from 
the complete WFby the scheme *c'a/(4A) *(zi + A, . . . , 2:4 + A)/*(zi, . . . , Z4)'^cm{Q)- Note that 
this result fully matches what we expect from analytic considerations, a in Fig. 5a. 

If we calculate quantities like the density or correlation function in the state ^' = '^r'^cM, we evaluate 
integrals of the type 



(53) 



"-*.*cm(^)= / dzi...dz„|^'r(2;i,...,z„)| \^cm{zi + ... + Zn)\ 5{zi - z) . 
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Fig. 22: Left to right: (1) cm{Z)\^ , Z = zi + Z2 + zs + Z4 for one of the three degenerate ground states in a 
4 particle u = 1/3 system and (2) the density n(r) of this state. (3) |*cM(^)r + I*cm(-^)I^ + I*cm(^)P of 
those three states and (4) the sum of their densities (divided by three). Note that the last density is nearly constant (as 
it should be for the Laughlin state) and thus by adding up densities of the three states differing only in the CM part, we 
eliminated the effect of the CM part of WF. 



Recast in CM and relative variables, this integral is a multidimensional convolution of and ^/r- 

Assuming that is isotropic but non-constant which is true for the Laughlin WF, the function n{z /4) can 
be thus shown to have the same periodicity as Less exactly but in more illustrative terms: n{z) 

is basically a smeared |^'cj\f(4z)p. Note that this explains why n{z) varies much stronger along x than 
along y. 

These considerations can be summarized in the following way. Even though e.g. the Laughlin state is 
translationally invariant, the CM part of the wavefunction which is always present in the exact diagonal- 
ization studies, will cause the density to be inhomogeneous. Consider eigenstates of J. This determines 
the form of to be as in Fig. 22a, cf (30) the density n{z) of an A^e-electron state (a) will be 1/Ne 
periodic along x, Fig. 21 or Fig. 22b, (b) will be l/SiVg periodic along y, Fig. 22d, (c) will be modulated 
much stronger along x than along y (compare scales in Fig. 22b and 22d) and (d) will rapidly converge to 
a constant for A^e ^ oo. 

Similar ideas have first been presented by Haldane and Rezayi [32]. 
How to suppress the effect of the CM part of the WF 

To suppress the effect of the CMWF it would be ideal to calculate the density as 

n^iJ,{z) = y"d-2i ■ ■ ■dznl'^rizi, . . . , z„)p(5(zi - z) 

rather than by (53). In other words, it would be nice if we could replace ^'cm(^i + . . . + Zn) by a constant 
in the numerically calculated wavefunction ^'^^'cm- 

Even though we could numerically calculate "i/cM and then calculate the density in the state / 'if cm, this 
is technically quite labourious and requires numerical evaluation of (n — l)-fold integrals. Instead we can 
make a trick. Consider again the example of the ly = 1/3 GS. The state is triply degenerated in the CM part 

12 3 

and the three different "^(j^'j (as they come from ED in subspaces with sharp J) have the pleasant property 
that the sum of their squared moduli is nearly constant, or in a more restrained (and honest) terminology, 
its variations are much weaker than those of individual Fig. 22. 

With this in mind we expect that the sum ^^i.^^^*^ (z) + (z) + n^-^^^qr^ (z) will be a good approx- 

imation to n<s,^ (z). The reader may check with Fig. 22 how well this is fulfilled. 

Other finite size effects 

Here, we will try to abstract from the effects due to the CM part of the calculated wavefunctions. Since the 
operator for density-density correlation depends only on relative coordinates we expect that g{r) will be 
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free of the finite size effects described in previous paragraphs. Since the curves for g{r) obtained from the 
z/ = 1/3 ground state in systems of different sizes (Fig. 12b) match very well for \r\ going at least to one 
third of the elementary cell we may have good confidence in these results even within the scope of infinite 
systems. In clear terms, we may believe that g{r) of the infinite system is nearly the same as g{r) obtained 
in a finite system (with a : = 1) as far as up to r sa 0.35a. 

Another type of finite size effects which are 'finer' than those originating from the CM part of the WF is 
shown in Fig. 22, the rightmost plot. The density plotted should be constant after averaging over the three 
states degenerated in the CM part in the infinite system. The weak (1/A'^,„) -periodic structure {N,n = 12 in 
Fig. 22) which we still observe reflects the quantization of one particle momenta by the PBC. One particle 
can be localized only around one of A^„i discrete set of points in the x-direction. This effect is the same 
along X and y, since we have lost the quantum number J (30), by averaging over the three states, belonging 
to J = 2, 6, 10 in the present case. Note, how extremely small this finite size effect is. 

2.1.5 Conclusion: yet another comparison to composite fermion models 

For a large part we were concerned with the u = 1/3, 2/3 and 2/5 incompressible ground states in 
this section. All these states, including their possible spin polarizations, can be described in terms of 
Landau levels (LL) filled with composite fermions (CF), Fig. 9a and Sec. 1.4. In particular, wavefunctions 
suggested by Jain, Subsec. 1 .4.2, are very close to the many-electron ground states calculated by exact 
diagonalization, as it is demonstrated by comparing the wavefunctions calculated by the two approaches 
in terms of overlaps which approach unity [77] or of correlation functions shown in this Section, Figs. 11, 
12 and 14a. 

However, we have seen in this Section that this picture is not as intuitive as someone may believe. Correla- 
tion functions of states with p filled CF LLs are quite different from those of states with p filled electronic 
LLs. Changing orientation of the effective magnetic field following from the CF LL 'quantization' alters 
the correlation functions drastically. It is hard to establish a relation between the ground states at = 2/3 
and 2/5 on the level of comparing the electronic correlation functions. We should also mention a discrep- 
ancy in the CF model for the i/ — 2/3 polarized state. It is both a particle-hole conjugate to the i/ = 1/3 
Laughlin state and a state with two filled CF Landau levels and effective magnetic field antiparallel to the 
real magnetic field or attached flux quanta. As Wu, Dev and Jain [77] noted already in their original work 
about antiparallel flux attachment, these two approaches give two non-equivalent microscopic wavefunc- 
tions. Surprisingly enough, both wavefunctions have high overlaps (w 0.99) with the polarized ground 
state obtained by exact diagonalization [77]. Thus, either both models are in fact indeed equivalent or this 
result shows that even such high overlaps may be not enough to prove the correctness of a trial many-body 
wavefunction. 

Another point worth of notice is that the 'CF cyclotron energies' (sometimes denoted by Hwcf) extracted 
from exact diagonalization with electrons are not quite the same in 2/3 and 2/5 systems. Fig. 23, the 
scaling factor 5 : 3 makes Bf-g equal in both systems. In the picture of non-interacting CFs, only the 
direction of the effective field B^g is reversed. Thus, if B^g has the same modulus in both cases and 
Zeeman energy vanishes then Ep{Ne = 8) — Eu{Ne = 8), i.e. the difference of energies of the polarized 
and singlet GSs for 8-electron systems, should be equal to four times the CF cyclotron energy in the both 
systems (cf. Fig. 9b). In the exact diagonalization spectra of iVg = 8 systems are regarded, the difference 
of hoJcF for V = 2/3 and 2/5 is small, about 5%, Fig. 23. However, the quantitative agreement becomes 
worse when we attempt to extrapolate the energies to larger systems. 

Also comparing 2/5 to 2/3, differences in the lowest excitations from the polarized and singlet ground 
states (energies, quantum numbers) are quite apparent. Fig. 23. 

All these facts demonstrate that it can be misleading to think of the 2 /3 and 2 /5 states as of an exact copy 
of Landau levels completely filled with electrons. Composite fermion models must be taken seriously since 
they provide us with many very good predictions (explicite forms of wavefunctions, e.g.) but apart of that 



adp header will be provided by the publisher 



43 



Fig. 23: Spectra of the v — 2/3 and 2/5 systems with 
eight electrons, zero Zeeman energy (SRI). Energies have 
been shifted to set the singlet GS in both systems to 
zero energy. In the picture of non-interacting composite 
fermions, the energies of the polarized (S = 4) ground 
states (the 2/3 and 2/5 ones) should then be the same 
when Beff is correctly rescaled (the indicated scaling of 
5/3). For each spin, the levels are sorted according to the 
orbital quantum number J (30), which can take on three 
non-equivalent values (left to right): 0, 1, 2. 

they are not exact, they fail to describe some phenomena Uke e.g. position of zeroes in Coulomb interacting 
states, Subsect. 2.1.2, the analogy between electronic and CF Landau levels is sometimes weak. One of 
the inherent problems not mentioned so far is the question of mixing between CF Landau levels: whereas 
LL mixing can be neglected for electrons in the limit B oo, there is no such case for CFs. 
The nature of many incompressible FQH states is therefore still not completely clear, for example the 
ground states at = 2/3 and 2/5. Results in this Section indicate that the singlet states at these filling 
factors comprise of pairs of spin up and spin down electrons which we would not expect from the CF 
analogy - at least not at first glance. Furthermore, in the i/ — 2/3 singlet with electron density n, the t — i 
pairs seem to form a state which could be constructed by taking a system with the lowest LL completely 
filled with electrons of density n/2 and then replacing each electron by an | — j pair. This behaviour is 
not observed in the i/ ~ 2/5 singlet. We may again conclude, that even though the 2/5 and 2/3 ground 
states are very closely related on the level of composite-fermion theories, their electronic properties are 
different. It can thus be misleading to extend our intuition concerning the (completely filled) electronic 
Landau levels to states interpreted as (completely filled) composite fermion Landau levels. 

2.2 The half-polarized states at filling factors 2/3 and 2/5 

In the previous section we dealt with the spin singlet and polarized ground states at filling factors 2/3 
and 2/5 and it was mentioned that it is the Zeeman splitting (or better, Ez / Ec oc a/B) which determines 
which of them is the actual ground state. It is the singlet state for vanishing Zeeman splitting (low magnetic 
fields) or the polarized state if the Zeeman term dominates (limit B — > oo). All this can be understood 
within the composite fermion concept. Fig. 9a, where we even obtain the prediction that there is a direct 
transition (crossing) between these two ground states at some critical value of Ez /Ec or equivalently, at 
some critical magnetic field Be, if we sweep magnetic field and keep the filling factor constant, cf. also 
Sect. 3. 

However, experiments by Kukushkin et al. [46] indicate that this picture may be incomplete. They suggest 
that some exactly half-polarized state becomes a stable ground state in the vicinity of Be- In this Section 
we wiU describe one candidate for such a half-polarized state ground state and discuss its properties. 

2.2. 1 Ground state energies by exact diagonalization 

At first glance, spectra of homogeneous small finite systems with Coulomb interaction (Sect. 3, Fig. 48) 
do not suggest any intermediate state at the transition. The picture is quite different when short-range 
interaction is considered. In an interval of magnetic fields around Be the GS is a state with total spin equal 
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Fig. 24: Ground state energies alv = 2/3: the half-polarized state may become the absolute ground state in a narrow 
interval of magnetic fields. Left: Energies of SRI ground states in the subspaces 5* = (singlet), 5* = Ne/4 (half- 
polarized) and 5* = Ne/2 (fully polarized) as a function of Zeeman splitting (or magnetic field; note the energy 
units e^/e^o oc ^/B). In all cases, energies of the two largest systems available to calculations are shown. Middle: 
extrapolation of the GS energies to infinite systems (1/7V —> 0). Right: The energy- versus-magnetic field diagram for 
extrapolated ground state energies. This indicates that even then the HPS will be a ground state close to the transition. 



to Fig. 24a, i.e. a half-polarized state (HPS). This holds for all system sizes accessible to numerical 

calculation and, by extrapolating energies to ^ 0, Fig. 24b, it seems to hold also for infinite systems. 
It is probably only through the finiteness of the system that a half polarized ground state did not appear in 
Coulomb interacting systems (Fig. 48). The SRI systems may be less sensitive to this generical drawback 
of exact diagonalization models. On the other hand, SRI models predict wrong values of Be (see Subsec. 
3.2. 1) and thus the scheme presented in Fig. 24a must be checked in systems with Coulomb interaction. 
Considering Coulomb-interacting systems, the scheme suggested in Fig. 24a is supported by extrapolations 
of GS energies performed by Niemela, Pietilainen and Chakraborty [57] in spherical geometry. Fig. 25a, 
and it is not supported by analogous calculations on a torus presented here. Fig. 25b. We would like to 
stress that the extrapolation of the energy of the HPS is based only on two points, the third point {N^ = 4) 
in Fig. 25b, is not very reliable, Subsec. 2.2.3. Therefore the question of whether the HPS becomes the 
absolute GS or not remains basically open until exact diagonalizations of larger systems become possible. 
Nonetheless let us assume in this Section that a half -polarized state can indeed lower its energy sufficiently 
so as to become the absolute ground state. We will therefore focus on the S = sector of systems 

at filling factor 2/3, and also at 2/5 in Subsec. 2.2.5. Studies were mostly focused on the SRI states 
where it is easier to identify the best candidate for the half -polarized ground state. Its Coulomb-interacting 
counterpart is discussed later, in Subsec. 2.2.6. 

By convention a half-polarized state with 12 (8) electrons will consist of 9 (6) electrons with spin up 
(majority spin) and 3 (2) electrons with spin down (minority spin). 



2.2.2 Identifying the HPS in systems of different sizes 

Provided some particular physical half-polarized state GSoois the ground state in an infinite system, we 
may ask what its realizations in finite systems of different sizes are. Vice versa: given the half-polarized 
states calculated in a system of A^e ~ 12 (4, 8. . . ) electrons, which state corresponds to GSoo? In this way 
we can think of states which 'correspond to each other' in systems of different sizes. The trouble is, of 
course, that we do not know GSqo- 

Regarding the computational capacity available, we could study v = 2/3 systems with 4, 8 and 12 particles, 
the next larger system, = 16, would require diagonalization in spaces of dimension many hundred 
million. It seems likely that the analogues to GSooare the GSs in = 12 and = 8 systems (GS12, 
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(a) (b) 
Sphere Torus 




1/N 1/N 

Fig. 25: Extrapolation ^ of the GS energies for Coulomb interacting systems on a sphere and on a torus at 
V = 2/3. 

GSs) and that it is a low lying excited state (st03) in the smallest system, iVg = 4. In the following, 
reasons for this are proposed. 

(i) GSi2 and GSs belong to the same symmetry class defined by the 'crystallographic k"^' (27). They have 
both k"^ = (tTjTt), i.e. they lie in the 'corner of the Brillouin zone' (Fig. 6, Subsec. 1.5.2). This is also 
closely related to the fact that both GS12 and GSg are non-degenerate. 

(ii) The states GS12 and GSs are well separated from excitations within the S = Ne/i sector and the 
energy of the lowest excitation is similar, 0.01 (e^/eio), in systems of different size. Fig. 26. 

(iii) Though not completely identical, the inner structure of GS12 and GSs is very similar as seen by the 
correlation functions. Fig. 27. 

(iv) The GS of the iVe = 4 system has a lower symmetry than the formerly described states. Looking 
for a state of inner structure (correlation functions) similar to the one of GS12 and GSs within the sector 
k"^ = (tt, tt), we find remarkable similarities with the second excited state ('st03', marked in Fig. 26), 
Subsec. 2.2.3. However, we should bear in mind that for = 4 there is only a single electron with 
reversed spin in other words the system is indeed extremely small. A consequence is for example that 
gil{r) = 0. Relevance of such states with respect to infinite systems is thus doubtful. 



2.2.3 Inner structure of the half-polarized states 

Focus of this part will be the correlation functions of the states GS12 and GSs and a brief comment will be 
made on iVg = 4 states. As mentioned above and as the kind reader may verify in Fig. 27, GS12 and GSs 
look indeed similar. 

GS12 and GSs match in all three spin-resolved correlation functions, (7||(r), gii{r), gii{r), Fig. 27. The 
match is especially good (quantitative) on short distances, r < 3£q. This suggests that states GS12 and 
GSs are not bound to some particular system size and we can thus hope that if we could make the system 
larger, they would eventually develop into the GStxj- 

Differences between correlation functions of GS12 and GSs at longer distances r are understandable, given 
the normalization (44). The Nf, — 12 system is 'larger' than the Ne — 8 one, yet the integral J drg{r) must 
be the same. Perhaps the most apparent difference between various correlation functions is whether they 
have a maximum or a minimum 'in the middle' (6^0 or 8^0 in Fig- 27). In ideal case, a strong maximum 
occurs when g{r) is monotoneous in an infinite system while a minimum, or a weak maximum following 
a foregoing minimum, means that g{r) has some structure, one or more maxima for finite r^. In reality, 
however, the former behaviour occurs also when r; is larger than the finite system size. A manifestation of 
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Fig. 26: Low lying energy levels in the S = N^/i (half-polarized) sector of 2/3 systems with (left to right) — 4, 
8 and 12 particles. The states are sorted according to \ k^\ (Subsec. 1.5.2). The GS of an infinite system is likely to 
have = (tt, n) which is a point of very high symmetry in the /(' -space. This symmetry class is however distinct 
from the one of the singlet and polarized incompressible ground states. 



this is seen in g^iir), Fig. 27. The flat maximum at r.i « 5£o, foflowed by a minimum, observed for the 
Nm = 18 torus does not occur for the smaller system (Nm = 12). Looking only at the smaller system we 
could have wrongly concluded that the correlation function is almost structureless. 
Some further points are worth of notice. 

(i) g^iir) is suppressed nearly to zero at r = in spite of the missing Pauli principle, only on account of 
the repulsive interaction. It displays strong maxima around r w 3.4f o- 

(ii) Even though by far not identical, (r) and (r) are similar to each other. The clear shoulder around 
r « 2£o seems to stem from the 'exchange hole' (of the LLL) (7^=i(r) = 1 — cxp{~r'^ /2£q), see (45). 
After subtracting a suitably scaled function gy^i{r) the shoulder completely disappears and the remaining 
parts of both .gtt (r) and gii{r) are cx close to r = 0, Fig. 28 and discussion below. 

(iii) Up to a high precision the sum of g||(r), gii{r) and g]i{r) (with appropriate scaling, see Fig. 15 for 
explanation) is identical with gu=i {r), however with £o replaced by v^^o- Not shown here. 

Let us now turn to the smallest system where S = states may occur (at i' = 2/3), i.e. = 4. 

Figure 27 shows correlation functions of the lowest two states in the sector of = (tt, tt). Out of these, 
the second state (i.e. st03) seems to be analogous to 5* = Ne/4: OS's in the two larger systems {Ne = 8, 
12): is again a sum of the 'correlation hole' and a function cx r^, guir) shows a peaked structure 

with maximum around 2.8£q (both of these features are missing for the lower state st02). However, as 
mentioned above, the Ne = 4 system is too small for a reliable study of 5 = Ne/4: states {gii{r) = 0). 
Back to the GS12 (called HPS here), it is very interesting to study the 'cx part' (P6P) of the like-spins 
correlation functions, g||(r), gii{r). What we mean by 'P6P' is the rest after we subtract the 'lowest LL 
correlation hole', i.e. the 5^=1(7') part causing the shoulder in gaa{f) around r w 2£o, curve A in Fig. 
28a,b. 
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Fig. 27: The likely analogues of the half-polarized (S = Ne/4) GSoo(ground state in an infinite system) on a torus 
with Ne = S and Ne — 12 particles and filling factor 2/3. Left to right: density-density correlation for t| (majority 
spins), and H (reversed spins), upper row: in the whole primitive cell; lower row: sections along the diagonal. 
Note the isotropy of the state, i.e. visual manifestation of its high symmetry. 
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Fig. 28: The half-polarized state (GS in the S — Ne/4: sector) for A^e = 12 particles, analysis of the correlation 
functions (details in text). Left to right: (a) (minority spin), (b) g-^^ (majority spin) and (c) gi^^. Legend for 
particular curves (in all plots). A: lowest LL correlation hole, gv=i{r) scaled to fit the shoulder. B and D: gii{r) 
and 3tt(r) without shoulder. C and F: gii{r) and gn{r) of the singlet state (A^e = 10, i/ — 2/3), gii{r) without 
shoulder. E: g{r) of the u = 1/3 Laughlin state. 'Without shoulder' means, that the curve A was multiplied by a 
suitable constant to fit the shoulder and subtracted. 
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One of many facts we can extract from Fig. 28 is that P6P/ tt [refering to g-^^ (r)] and P6P/ H are similar 
but not identical. For example, they both exhibit a peaked structure but the first maxima do not coincide, 
they occur at S.O^o and 5. 8^0 for P6P/ tt and P6P/ H, respectively (curves B and D in Fig. 28a). 
Let us compare the P6P/ ii of the HPS with P6P/4 j of the singlet incompressible 2/3 GS, curves B and C 
in Fig. 28a). Match of these two is very good up to r w Mq, the absence of the peak at 5. 8^0 in the singlet 
state could be due to smallness of the system where the singlet state was determined {N^ = 10). It might 
appear in the next larger system, = 12, cf. similar situation in Fig. 27. 

On the other hand, P6P/tt of the HPS seems to resemble the singlet state less than P6P/i j of the HPS. The 
form of P6P/tT seems to be not very different from the one of the correlation function of the Laughlin 1/3 
state. Fig. 12b, whose first maximum occurs however already at r = 4.4^o (curves D and E in Fig. 28b). 
In any case, P6P/tt of the HPS matches better gv=i/3{f), ie. the Laughlin state, than P6P of the = 2/3 
singlet state. Here we mean especially behaviour on ranges < 3£o. 

Last but not least, the correlations between unlike spins are also very similar in the singlet state and in the 
HPS, Fig. 28c, in particular positions of the maxima differ by as little as O.l^o (both are around r w 3.4£o)- 

2.2.4 Discussion 

Findings presented above suggest that the = 2/3 half-polarized ground state in short-range interacting 
systems is a gapped state in which the singlet and polarized incompressible states coexist. Below, some 
key points regarding the HPS are summarized. 

Symmetry and energy 

Both in eight- and twelve-electron systems, the ground state has k'' = (tt, tt). This is one of two points 
of the highest symmetry in the /(''-space, another one is k"" = (0, 0), Fig. 6. In particular, the 'highest 
symmetry' means that this /c ''-point is not related to any other point by a symmetry operation in the k^- 
space corresponding to relative translations. Sec. 1.5.2. This in turn implies that states with /c'' — (tt, tt) or 
(0, 0) - and only such states - are non-degenerate, except for center-of-mass and incidental degeneracies. 
Together with the relatively large lowest excitation energy A{Ni, = 8, 12) from both GS12 and GSg (10% 
of the gap of the Laughlin state. Fig. 26) , this suggests that the ground state is gapped. Also the relation 
A{Ne = 8) < A(iVe = 12) speaks in favour of this hypothesis. If the gap were to vanish in an infinite 
system, we would expect the lowest excitation energy to decrease with system size. Naturally, we must 
be careful, since we can compare systems of only two different sizes and the function A{Ne) may be 
non-monotonous. On the other hand, A(A''e = 12) w 0.01 (e'^/elo) is much larger than a typical level 
separation between excited states. Fig. 26 and for a mere finite size effect, this gap seems too large. 

In spite of the similarities to the singlet and polarized incompressible ground states, k '' clearly distinguishes 
HPS from these two states, since they have both k^ = (0, 0). Also in spherical geometry, where |/c''| oc L 
(end of Subsec. 1.5.2), these incompressible states have L = while the HPS has L ^ S, where S is the 
total spin [57]. Thus, even though we showed that the HPS could be gapped, it is of different nature than 
the singlet and polarized ground states. Meaning of this different symmetry is however not clear. 
It would be interesting to study this state in a system with hexagonal elementary cell [31] which was 
unfortunatelly out of the scope of this work. This geometry is nearer to an isotropic 2D system than a 
torus (it has a six-fold rather than a four-fold rotational symmetry) while it is still compatible with plane 
waves (in CDWs). Most importantly, there is only one point of the highest symmetry in this geometry and 
a straightforward question is whether or not the HPS will maintain its high symmetry. 

Inner structure again 

Features of the HPS described by points (i-iii) in Subsec. 2.2.3 are actually strikingly similar to those of 
the incompressible singlet state at 1/ = 2/3. Investigation of the g^-^ (r) after the 'shoulder' was subtracted 
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(P6P/tt) suggests again some relation to the Laughlin state which is the particle-hole conjugate to the 
polarized incompressible state at = 2 /3. 

Especially manifest is the hint at pairing between unlike spins, the maximum around 3A£q in g^i{r). On 
the other hand, the shoulder in correlation functions of like spins seems to be rather a manifestation of 
filling factor > ^, since it occurs also for other states at filling ~ 2/3 (than just for the singlet, polarized 
and half-polarized GS) and it does not occur at filling v = 2/5 < ^, Subsec. 2.2.5. It suggests that some 
ly = 2/3 states with less than full polarization can be interpreted in terms of holes rather than electrons 
even though particle-hole symmetry applies only for fully polarized states, Subsec. 2.1.1. 
In the following Sections we will continue investigating the half-polarized states at filling factor 2/3 by 
other methods and continue discussing the hypothesis of coexisting singlet and polarized states. First, 
however, we look at two different minor issues. 



2.2.5 Half -polarized states at filling v = 2/b 

At filling 2/5, the situation is much less transparent than at filling 2/3. First, only systems with four and 
eight particles are accessible to exact diagonalization, the twelve particle system implies matrix dimensions 
in the order of hundreds of millions. Second, the spectrum of the eight particle system in the S = N^/^^ 
sector is quite different from that of a 2/3 system. Fig. 29: 

(i) the ground state lies at a different point in the /c ''-space, (0, 0), than the 2/3-HPS having /c*" = (tt, tt). 

(ii) The excitation energy from this GS is very small, less than a third of that one of the 2/3 HPS. 

(iii) The symmetry of the low excited states is lower than for A^g = 8, 2/3 system. 

Regarding the possibility that within the 8 electron calculations it is not the lowest energy half-polarized 
state at 1/ = 2/5 to be the counterpart of the HPS at = 2/3, there are two 2/5 states displayed in Fig. 30: 
(a) the one with the lowest energy in 5 = A*'e/4 sector and (b) the lowest state with the same symmetry as 
the 2/3 HPS, i.e. k = (7r,7r). 

Similarly, as for the = 2/3 states, the 2/5 HPS bear features of the polarized and singlet ground states. 
Let us regard the state in Fig. 30a: 

(i) Near r = the functions gn (minority spin), g-^i and (majority spin) are cx r^, and r^, respec- 
tively. In this respect, gn and g^i resemble the singlet state and g^j resembles the polarized state. 

(ii) Up to the first maximum, (/|| of the HPS is the same as in the polarized state, but shifted by about 0.2^o 
outwards. Positions of the first maxima mismatch slightly more (by 0.4£o)- The strong maximum in the 
centre of the cell is not present in the HPS. 

(iii) gii of the HPS and the singlet GS match very well even beyond the first maximum. Positions of the 
maxima are identical, r 3.5^o- On contrary to the previous point, there is another maximum in the centre 
of the cell in the HPS state and nothing in the singlet state, indicating that the similarity between the singlet 
and the HPS has certain limits. 

(iv) gii of the HPS and the singlet GS also match very well up to r « Then there is a deep minimum 
in the HPS which is absent in the singlet GS. 

Turning to the state (b), we might say that it is less alike to the singlet state. The minimum in gn is much 
deeper than for state (a), the first maximum in g-^i does not match the maximum seen in the singlet state. 
On the other hand, f/|| seems to be more similar to the polarized state. 

Lowest excitations in the high symmetry sectors show even less similarities to the singlet and polarized 
GSs, especially gn is quite dissimilar beyond the r w range and maxima in match less well. 
In conclusion, if there is a counterpart to the 2/3 HPS at filling 2/5 at all, we may expect it to be the state (a) 
(the absolute GS), even though hints for this are not very convincing. Again, this suggests that differences 
between filling factors 2/3 and 2/5 are not only of quantitative nature (gap energies, for instance) but 
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Fig. 30: Half-polarized states at filling 2/5 and their inner structure (density-density correlation functions); eight 
electron system, (a) the GS in the S — Ne/4: sector is non-degenerate but it has a different symmetry, i.e. (0, 0), than 
the HPS of v — 2/3. (b) the lowest half-polarized state (at 2/5) with the same symmetry, i.e. (|-, ^) as the HPS of 
v — 2/3 (marked by 4|k in Fig. 29). 



may be as substential as existence or non-existence of some particular ground state, which is nota bene 
completely unexpected on the level of non-interacting CF picture. 



2.2.6 Short-range versus Coulomb interaction 

Let us conclude with observations regarding the Coulomb- and short-range-interacting (SRI) systems in 
the sector of half -polarized states. 

The spectra do not look very similar. Fig. 2.2.6a. However, the absolute ground states have in both cases 
the same symmetry, they lie in the same point of the k space. 

The Coulomb and SRI ground states in the largest system available, N,, = 12, have very similar structure. 
The correlation functions g|| and g-^i match nicely while gn show some differences between the CI and 
SRI states. In spite of this, the overlap between the two states is as large as 95%. This allows for the 
following conclusions 
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Fig. 31: Half-polarized states in Coulomb interacting systems. 



(a) The two states 'correspond to each other', (b) The short-range part of the interaction seems to be 
essential for this state (very similar as for the Laughlin state), (c) Deviations in gn (minority spin) might 
come from the fact that spin-down electrons are very far separated from each other (they have an effective 
filling of only v ~ i). Thus the long-range part of the interaction substantially influences their motion. 
In TVe = 8 systems, the most likely analogue to the Nf, = 12 ground state is the state (), Fig. 2.2.6a. This is 
the lowest 8-electron state with the same symmetry (value of k^') as the = 12 ground state. Correlation 
functions of the two states (8- and 12-electron ones) match reasonably. Fig. 31(b). Compare also with 
differences between A^e = 8 and ~ 12 short-range ground states. Fig. 27. 

Among the excited states the level order is often modified, comparing the Nf. = \2 Coulomb and short- 
range systems. When trying to assign CI to corresponding SRI states, calculating overlap between two 
states seems to be a more reliable tool than comparing correlation functions. 

In summary, in spite of differences in the excitation spectrum, the half-polarized ground states of Coulomb 
and short-range systems seem to correspond to each other Differences in the excited states and in the 
correlations between the minority spin electrons indicate that the definition of the short-range interaction 
should be improved when we study the half-polarized states. Since the minority spin electrons are relatively 
far from each other, non-zero values of higher pseudopotentials {Vm, m > \, Subsect. 1.3.5) should 
probably be considered. 



2.3 In search of the inner structure of states: response to delta impurities 

Now that some candidates for the half -polarized ground state at filling 2/3 have been introduced we wish 
to look at them more closely and learn more about their inner structure. The ultimate goal of such efforts 
can be to propose trial wavefunctions just as the Laughlin wavefunction at filling v ~ 1/3. Even though 
this has not been accomplished, the results presented below shed some light on relations between the 
half-polarized state and the singlet and polarized incompressible states. 
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(a) Density response: incompressible- or com- 
pressible-like. Note that density integrated over 
the hatched area remains unchanged for the in- 
compressible system when the inhomogeneity is 
switched off. 
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the type 5 



(c) Two possible ways of how a response 
to an inhomogeneity can change with 
system size; different lines in one sketch 
refer to the same state in systems of dif- 
ferent sizes. 



Fig. 32: 5-line inhomogeneity and sketches of possible resulting effects. 



As a probing tool, the homogeneous states are subjected to a (5-line impurity and response in density and 
polarization is observed. In first quantization, 

Hin^purity ^^W{r,) , W {x , y) = 5{x ~ xq) (54) 

i=l 

This inhomogeneity profile (Fig. 32) was chosen since it is compatible with the torus symmetry. For studies 
of point-like impurities, spherical geometry is more suitable since it preserves the rotational symmetry, cf. 
references in Subsec. 2.3.1. The (5-line form is particularly apt to unveil a tendency of the state to build 
plane charge or spin density waves. We should keep in mind, that due to the restriction to the lowest 
Landau level, even a (5-like potential has an effective cross section of £o [61]. 
As we are dealing with spinful electrons, inhomogeneities can be principially of four distinct types; 

Hei = W{r) ■ ((5,-r + (5,J , Hmi,^ = W{r) ■ S^^ , 

Hmi = W{r)- ((5,-r Hmia =W{r)- 5,^ , (55) 

where the function W{r) describes the spatial form of the impurity. Fig. 32 shows the form of W{r) chosen 
in the present study. It is important to note that these impurities fail to conserve but they do conserve 
Sz- Also, owing to the form of W{r) — W{x), they conserve /(^ and thus also J (30) and they spoil only 
the /c ^-symmetry. This is very convenient from the computational point of view as matrix sizes remain 
tractable. From the physical point of view, this inhomogeneity is a soft tool which does not completely 
destroy the high symmetry of the studied states. For example, it allows us to stay in the Sz = iVe/4 sector 
when we study the half-polarized states. 

The first type {Hei, electric impurity) is an ordinary non-magnetic impurity or external electric potential. 
The magnetic impurity (Hmi) favours particles with correct spin (|, if W{r) > 0) and costs energy for 
particles with wrong spin (f in this case). The last two types describe an impurity which is seen only by 
one group of spins. In case that a system consists of two separated subsystems, one of spin up particles and 
another of spin down particles, these impurities allow to test only one of them without directly disturbing 
the other one. 

Note that some inhomogeneity types in (55) may be redundant, depending on the state we apply them to. 
For instance, the effect of Hmi,^ and Hmi.i must be the same up to a sign for all states in the 5r = 
sector. 

Before we turn to the exact diagonalization results, let us briefly think about what types of responses can be 
expected. Most importantly, consider the difference between compressible and incompressible states. As 
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a classical compressible system imagine a playground of fixed size filled with a gas of negatively charged 
footballs of density n{r), {n{r)) ^ N. A negative impurity at r = will repel the gas causing n(0) < N, 
Fig. 32a right. Beyond some distance r^, the density will reach a constant level again and this level will be 
slightly higher than the original density, N + S = n{r) > A^, |r| > r^, so that the constraint {n{r)) = N 
remains preserved. Some charge has been depleted away from the impurity, thereby compressing slightly 
the gas in the rest of the system. If the depleted charge equals the charge of the impurity, the charge 
distribution (charge of the footballs plus charge of the impurity) in the system will remain constant in spite 
of non-constant n{r) and the gas particles far away from r = will not 'see' the impurity anymore. This 
is the case of ideal screening. 

A classical incompressible liquid, for example again charged footballs, will not react at all, because it 
cannot change its density. Even though particles of the liquid feel repulsion from r = 0, the density will 
remain constant ??(r) = N. We can also encounter a different behaviour. Fig. 32a (left). Though the 
density decreases directly at r = 0, an oscillatory structure develops in n(r), so that the integral density in 
the region \r\ < remains as it was without the impurity. The density then also remains at its original 
value beyond rh- This is a non-ideal incompressible behaviour: at very short distances, the density can 
vary slightly, but averaged over distances of (or larger), the density remains constant. Also, since no net 
charge was depleted from the region |r| < r^, the impurity is completely unscreened on distances larger 
than Th. 

Compressible-like response as shown in Fig. 32a can be combined with quantum interferences (Friedel 
oscillations) and it is also possible to think of some overscreening effect which would lead to an oscillatory 
n{r). This means the sole fact that n(r) exhibits oscillations does not necessarily have to imply incom- 
pressibility. A more reliable criterion is that the integral density over | r | < remains the same with and 
without impurity. This procedure is delicate in finite systems where can be comparable to the system 
size. 

The last Figure, 32c, shows two possible ways of how responses change with system size. The right panel 
suggests that the state is not fixed to a particular size of the finite system and especially we could expect 
oscillations with period ri also in an infinite system. On the contrary, the left panel shows a state with no 
intrinsic length scale and e.g. the width of the peak is related to the (finite) size of the particular system. 
Now, let us proceed to fractional quantum Hall states. 

2.3.1 Electric (nonmagnetic) impurity 

The effect of electric impurities on incompressible ground states has been under investigation since the 
early times of the fractional quantum Hall effect. The main reason is that some disorder is needed for the 
integer quantum Hall effect to be observable, but on the other hand, too strong disorder will destroy the 
effect [79]. For the fractional quantum Hall effect, two of the basic questions were, (i) how strong impurity 
potentials may be so that they do not destroy the gap and (ii) how does it change the ground state. Basic 
studies with the Laughlin state were performed as early as in 1985 [61, 86, 27]. 

Since the exact diagonalization is limited to finite, and in fact quite small, systems, it is very delicate to put 
forward statements about the infinite 2D electron gas. Therefore, when we use the word 'incompressible' 
we mean rather 'incompressible-like' in terms of Fig. 32. In fact, the main purpose of the following 
Subsections is to see how the polarized and singlet state respond to impurities in afinite system and later to 
compare them to the half-polarized state again in afinite system. We will focus on short-range interacting 
systems here. 

The Laughlin state or the fully polarized 2/3 state 

The fully polarized i/ ~ 2/3 state is a particle-hole conjugate to the = 1 /3 Laughlin state in a homo- 
geneous system, Subsect. 1.5.4. In this part we will study the latter state. Strictly taken, the particle-hole 
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symmetry is lost when an arbitrary impurity is considered since the Hamiltonian is no longer translationally 
invariant. Differences between the v = Xji and 2/3 polarized states are however small if the impurity is 
weak. In particular, for inhomogeneities considered in this paragraph, it has been checked numerically that 
n[x) — Nf. are almost the same for the two states. Moreover, the larger N^, the smaller are the differences. 
The response of a = 1/3 system to an impurity of the form (54), a (5-line along y, is shown in Fig. 33. 
Different curves show the ground state density n(x) in systems of different sizes (iVg = 4 to 10 particles). 
The repulsive impurity is always located at a; = and it is weak, its strength is ^ 10% of the gap. These 
results agree very well with the densities presented by Zhang et al. [86], who considered a 5 rather than 
a (5-line impurity, though for A^'g = 4 systems only. Comparison between rectangular, spherical and disc 
geometry showed, in all cases, very similar behaviour [86]. Note also that findings in Fig. 33 assume 
short-range interaction whereas [86, 61] considered Coulomb interaction. 

Results in Fig. 33a support the conclusions of Zhang and Rezayi. The oscillatory response of n(x) is 
size-independent and it has a period r\ w Lhi^. The response, measured by n(0), does not vanish with 
increasing system size but it decays with distance from the impurity. Comparing n[x) in Fig. 33a to the 
model cases in Fig. 32a, we may tend to classify the Laughlin state as an incompressible one. Incompress- 
ibility of the Laughlin state is locally not perfect, otherwise n(x) would remain constant, at least in infinite 
systems. However, if some net charge were accumulated even in a larger region (of the order Vh) around 
a; = 0, we would expect n{x) at large distances to be consistently higher than the no-inhomogeneity value 
n{x) = Ne- Recall the difference N to N + 6 in Fig. 32a. If just a unit charge is depleted from the 
impurity, then 6 ~ 1/N . In infinite systems, the difference 5 will vanish, but data in Fig. 33a come from 
rather small systems iV < 10 where 5 is not negligible. This is not seen in Fig. 33a. 
Zhang et al. suggest that the observed response is a local charge density wave (28), a strong argument 
supporting this idea is given in point (iv) below. Under this view, it is not surprising that the response to a 
(5-line shown in Fig. 33(a) is very similar to the response to a (5-peak studied by Zhang. Only the envelope 
function, not the wavelength depends on the particular form of the exciting impurity. 
We should again add several comments: 

(i) oscillations observed in n(x). Fig. 33a, are not related to Friedel oscillations which appear in the Fermi 
gas. This is where a sharp Fermi surface exists giving rise to interferences, just as in correlation functions 
of a free Fermi gas, Subsec. 2.1.1. 

(ii) small wiggles on the A'^e = 4 density in Fig. 33a are due to the center-of-mass (CM) part of the 
wavefunction. Being a finite size effect, they fall off rapidly with system size as we indeed see in Fig. 33a, 
Subsec. 2.1.4. 

(iii) the ground state of the homogeneous system is triply degenerate in the CM part, Subsec. 2.1.4. This 
degeneracy is lifted by the inhomogeneity, but energy differences between these three states remain much 
smaller than their separation from the lowest excited states for the inhomogeneity strength considered [86]. 
The response n{x) of any of the three states depends slightly on the position of the impurity within the 
elementary cell, but this dependence and also differences among the three states in energy and in n{x) 
quickly vanish with increasing system size. In Fig. 33a always the impurity giving the strongest response 
in n{x) was chosen. 

(iv) Period of oscillations: As Rezayi and Haldane [61] note, numerical calculations as in Fig. 33 agree 
with results of the single mode approximation proposed by Girvin et al. [27]. The linear response function 
xiq) (in the v = l/i Laughlin state) is dominated by the magnetoroton collective mode around (/o^o ~ 1 -4. 
Would it be = — qo)^ the density response to a point impurity potential would be n{r) cx Jo((Zo'')- 
This density profile looks like damped oscillations with the first node at r = 1.7£o- 

Regarding a more realistic profile of x{(l)^ this estimate for n(r) is a very good approximation to n{x) in 
Fig. 33a. 
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(a) The polarized state (i/ = 1/3 considered, see text) 
in systems with Ne = 4-10 electrons. Note that the 
response does not decay with increasing the system 
size. 



(b) The singlet state. The thin line shows the CM os- 
cillations in a homogeneous system, thick lines show 
responses to an impurity with CM oscillations sub- 
tracted. 



Fig. 33: Polarized and singlet 2/3 state and (non-magnetic) impurity in the form of a 5-line (along y). Normalized 
density along x is plotted. 



As the purpose of the present work was to study systems with spin, we will now continue to spin sin- 
glet states at 1/ — 2/3. Some quite new results for the v = 1/3 Laughlin state have been achieved 
byMiiller [53]. 



The singlet state 

The 2/3 singlet ground state shows basically the same signs of incompressibility as the polarized state. The 
period of the density oscillations incurred by a (5-line impurity is almost the same (ri « "ito), and also in 
terms of classification of Fig. 32a, the singlet state shows an incompressible-like behaviour, cf. discussion 
of the polarized state. The striking feature of the singlet state is, that the strength of the response is about 
an order of magnitude less than in the polarized state. Thus in an 8-electron system, the density response 
is 'hidden' under the center-of-mass oscillations. Fig. 33b. 

This strong difference between the singlet and polarized ground states is unexpected since 'incompress- 
ibility' gaps of both states are similar. 

This hints at unusual stability of the singlet state with respect to charged inhomogeneities. In terms of 
perturbation theory, this is not due to energetic reasons but rather owing to small matrix elements of Hei 
between the ground state and excited states. Energy of the first excited state, however, decreases when 
impurities are present and thus, in spite of the quite stable density of the GS, the gap will eventually 
collapse. 

Regarding the response in systems of different size, we find a considerable attenuation when going from 
eight to ten-electron systems. Fig. 33b. Nevertheless we assume that the response remains finite even in the 
thermodynamic limit. To support this hypothesis, a fact worth of emphasis is that the = 8(10) singlet 
state occurs in systems with Nm = 12 (15) flux quanta, i.e. with system area A = (1)- These are 

the two smallest systems considered in Fig. 33a. For the these two systems we also observe a considerable 
attenuation of the n{x) response when going from the iVg = 4 to iVg = 5 state, Fig. 33a, and this reduction 
in response is definitely only a finite size effect. As close as this analogy is, observations presented in Fig. 
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33b are not conclusive and an investigation of the singlet state in a larger system (Ne ~ 12) would be 
needed. 

Let us just briefly mention, that non-magnetic impurities have no effect on the polarization of the singlet 
ground state. 

2.3.2 Magnetic impurity in incompressible 2/3 states 

As far as spin polarized states are considered, magnetic impurities (55) cannot have any other effect than 
the electric impurities do. Therefore only the 2/3 singlet ground state will be discussed here as the half- 
polarized states deserve to be considered separately, Subsec. 2.3.4. 

Considering the density. Fig. 34a, we find a yet weaker response than for non-magnetic impurities. Fig. 
33b. The response reminds of an incompressible system, in terms of Fig. 32a, and may remain finite in the 
thermodynamic limit, cf. discussion of non-magnetic impurities. 

Polarization ni{x)/n{x) behaves quite differently. Fig. 34b: the response is large and it looks compress- 
ible. Again in terms of Fig. 32a. In particular, note that the polarization (x) / n{x) in Fig. 34b approaches 
« 0.51 as we go 'far away' from the impurity, i.e. a different value than the polarization in the homoge- 
neous case, O.5.. Electrons with 'correct spin' (f) accumulate around the impurity, n ^{O) / n{0) drops from 
the homogeneous value (0.5) by as much as by 5%, whereas the average polarization off the impurity 
slightly increases so as to keep the overall average value 0.5 as required by Sz = 0. This behaviour differs 
strongly from the density response. Fig. 34b. 

It should also be noted that both density and polarization are here much less system-size dependent than in 
the case of non-magnetic impurities. 

These are quite remarkable findings. It seems that the singlet state is locally much more 'incompressible' 
than the polarized state. On the other hand, the singlet state is relatively easily polarizable which is partic- 
ularly striking when compared to the weak response in the density. If we assume the density in Fig. 33b 
to be the response of two independent liquids, then the polarization in Fig. 34b should be (i) smaller by 
a factor of five for = 8 than what is observed and (ii) considerably smaller for = 10 compared to 
the A^e = 8 case. This again contradicts the picture of two uncorrected 1/3 Laughlin liquids, one spin up, 
another spin down, which we could wrongly infer from the view of filled composite fermion LLs. Remind, 
however, that it is in fact not the claim of CF theories, that particles of n = 0, t and n = 0, J, CF LLs are 
uncorrected. 

2.3.3 Integer quantum Hall ferromagnets 

A brief introduction to integer quantum Hall ferromagnets (QHF) was given in Subsec. 1.6. It is instructive 
to keep in mind the scheme of Landau levels. Fig. 8. 

Here we will focus on the Sz ~ sector in prototypes of Ising and Heisenberg QHFs with neglected LL 
mixing. These states {Sz = 0) are analogues of the half -polarized states at filling 2/3, the explanation 
follows. Disciples of CF teachings deem the v = 2/3 ground states to have i/qf ~ completely filled 
CF LLs, Fig. 9b. Transitions between the singlet and polarized GSs occur, when the {n, a) ~ (0, j) CF 
LL crosses the (1, f) CF LL. It is then plausible to neglect the low lying (0, f) CF LL and look only at the 
two crossing CF Landau levels. The two ferromagnetic Ising states - the singlet, and polarized electronic 
GS at = 2/3 - correspond to all CFs placed in the (0, |), and (1, f) CF LL, respectively, compare to 
Fig. 7b. Hence the half-polarized state ~ 2/3) corresponds to half-filled (0, |) and half-filled (1, f). 
Disregarding the fully occupied (0, |) CF LL, i.e. counting only particles in the two crossing CF LLs (in 
total Ne CFs), the ferromagnetic Ising states are Sz — ±iVe/2 and the 'half-half state is Sz = 0. 
In this Subsection we study the same situation as the one occuring at the 1^ — 2/3 ground state transition 
(within the picture of crossing CF LLs) but for electronic Landau levels, i.e. with electrons instead of 
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Fig. 34: Singlet 2/3 ground state for an attractive magnetic impurity Hmj^-^ (55) in the form of a 5-line (along y). 
Normalized density and polarization along x are plotted. 



composite fermions. We therefore study a v ^ 1 system with spin degree of freedom, where spin down 
(spin up) electrons lie in the n = (n ~ 1) Landau level, respectively, and we disregard the fully occupied 
(0, t) level. The physical system we model herewith has thus ly = 2. Technically, this requires only 
implementing modified values of pseudopotentials. Fig. 2. Without electron-electron interaction, these 
two Landau levels are set to equal energy so as to model the LL crossing. Mixing to the fully occupied 
(0, t) LL as well as to all higher LLs is neglected, since all these levels are well separated from the two 
crossing levels. 

Heisenberg QHFs are not related to v = 2/3 and we investigate them just for the sake of comparison 
between Ising- and some other type of QHF. In the integer QHE regime, Heisenberg QHF occurs e.g. 
when (0, t) and (0, i) LLs cross (and v = 1) as it is the case for instance at vanishing Zeeman splitting. 
With CFs, this happens atiy= 1/3, i.e. i/cF = 1, Fig. 9b. 

We will first briefly discuss homogeneous states in these QHF systems and then we will turn to their 
response to magnetic inhomogeneities ((5-lines). 

Ising quantum Hall ferromagnet 

There are two degenerate ground states of an Ising ferromagnet: both with S = Ne/2, one Sz = Ne/2 and 
another Sz = —Ne/2. In general, excited states are no eigenstates to as a consequence of the omission 
of the fully occupied (0, t) level and may only be classified according to 5*2. They are all situated well 
above the ground states. Fig. 7b, and their energy grows with Ne/2 — \Sz\- In the following we will only 
speak about Sz — states. The whole Sz — sector is quite high in the complete spectrum. Unlike for a 
Heisenberg ferromagnet there is nothing like a S ~ Ne/2, Sz — ground state for an Ising ferromagnet. 
Low lying Sz = states of the considered Ising QHF are apparently arranged into a flat dispersion branch. 
Fig. 35a. For a fully occupied Landau level, J 'coincides' with ky. Precisely, ky ~ {Ne/2 — J)\j2'K /Nm 
for Ne even in the sense of (27). Centre-of-mass degeneracy is absent. The anomalous form of the branch 
in a iVe = 8 system, seems to be of finite-size origin, since Ne = 10, 12 and 14 spectra are all similar. 
States of the lowest branch have /c"" of the form {2-Kn/Ne, 0), n = 0, ±1, . . . , Ne/2, or (0, 2'Kn/Ne). This 
is in agreement with the symmetry between x and y (square elementary cell). It shows that rotational 
symmetry is absent in the low energy sector - otherwise we would observe also states with k"^ ^ (k^, ky), 
kx,ky 7^ 0. The lowest branch flattens and becomes well separated from excited states with increasing 
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(b) Left: Density in the ground state, 12 particle system; right: change in density in response to a <5-line impurity, 
cf Fig. 37b. 



Fig. 35: Half-polarized states (Sz — 0) of an Ising quantum Hall ferromagnet. ty — 2 



system size, and the minimum energy remains at k '' = (0, 0). Also, other branches develop, the second 
lowest branch is described by = {Trn/Ne, ±27r/A''e) (plus the x-y symmetric partner) and minimum 
energy at points (tt, ±27r/A^e), see the = 12 spectrum in Fig. 35a. Apart from these branches an 
isolated k '' — (0, 0) state is present (shown in gray in Fig. 35a) and it is hidden within the branch. We can 
hypothesise that this state becomes the absolute ground state i.e. gets separated from the lowest branch in 
sufficiently large systems. 
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(a) Spectra. (b) Density, 12 particle system. 

Fig. 36: Half-polarized states (5, = 0) of a Heisenberg quantum Hall ferromagnet. v — 1 



The flat branch is reminiscent of results of Rezayi^ et al. [63] and could correspond to domain states, i.e. 
stripes along x or y with alternating spin polarization, in a system which does not prefer any particular 
domain size. The origin of the highly symmetric isolated (gray in Fig. 35a) state is unknown. 
As it can be expected, low lying states have homogeneous density and it is true even for the whole lowest 
branch, as an example we show density in the ground state (/c*" = 0), Fig. 35b. The anticipated domains 
would probably be visible first in coiTelation functions. States in the second lowest branch show unidkec- 
tional charge density waves. 

In summary, in a homogeneous v = 2 Ising QHF we observe 

(i) a flat branch of low lying states, which could become degenerate in infinite systems. Fig. 35a. This first 
branch probably consists of stripe domains - or spin density waves - of all possible wavelengths A = a/n, 
n = 0, 1, . . . iVe/2 just as in the system studied in [63]. Contrary to isotropic states (like Laughlin liquid), 
the wave must be parallel to one side of the square elementary cell. 

(ii) second branch with pronounced dispersion, which could be a charge density wave 

(iii) continuum of excited states above the two branches and 

(iv) another state, with high symmetry, k"^ = (0, 0), which lies among the states of the lowest branch. 
Heisenberg quantum Hall ferromagnet 

The situation here is quite different from the Ising ferromagnets. The Hamiltonian (Coulomb interaction 
projected to the lowest Landau level) conserves the total spin S and it even commutes with 5"+ and 5*" 
which change 5*2 while keeping the length of the total spin. The ground state is fully polarized, S = Ne/2, 
but its z-component of spin is arbitrary. Fig. 7b. 

Looking at the sector Sz = 0, Fig. 36b, the lowest state is thus the ferromagnetic S = state. Other 

low-energy states form again a branch, k'' = {±mi/Ne,0) and (0, ±7rri/A^e), n = 0, . . . ,Ne/2 {x-y 
symmetry present, rotational symmetry absent). Contrary to the Ising QHF, this branch does not seem to 
flatten. States in the branch fulfil S ^ Np/2 — n: the ferromagnetic (ground) state is polarized and going 
up the branch, the polarization decreases. In this respect, the excitations of the lowest branch markedly 



The cited work concerns the situation when the lowest and the third Landau levels of different subbands cross. Rezayi et al. 
had first to show that this system is an Ising QHF. See subsection 1.6 for more details. In the Sz = sector of his system 

Rezayi et al. found a multiply (almost) degenerate ground state with k'~ just of the sequence {2TTn/Ne,0), similai' as we see 
in Fig. 35a for Ne = 12. 
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differ from spin density waves. What we observe in the Heisenberg QHF are most likely states with n 
weakly interacting spin waves which were observed under the same conditions on a sphere by Wojs and 
Quinn [76]. 

Half-polarized QHF states and magnetic impurity 

If a homogeneous state cannot be established and domains formation is more favourable, then no particular 
domain size is preferred. This is the central message of the following paragraph and it applies to both Ising 
and Heisenberg QHFs described at the beginning of Subsection 2.3.3. 

The two systems were subjected to a (5-line magnetic inhomogeneity, just as the incompressible singlet 
ground state in Subsec. 2.3.2. However, QHF systems and incompressible liquid states at u = 1/3 or 2/3 
behave quite differently. Looking at a QHF and comparing the response in systems of different sizes, we 
observe no intrinsic length scale. Fig. 37. Rather, the form of the response reflects the size of the system 
(cf. the left panel of Fig. 32c). This statement applies both to the Ising, Fig. 37b, and the Heisenberg 
QHF, Fig. 37a, where we show the polarization of the energetically lowest state in a system subject to the 
inhomogeneity. 

It is also interesting to study the density of the disturbed QHF states. The density of the Heisenberg QHF 
remains almost unchanged (it is constant) unlike the density of the Ising QHF state. Fig. 35b right. This 
is understandable. Whereas in the Heisenberg QHF spin up and spin down one-particle states have exactly 
the same wavefunction, both spin up and spin down states are from the lowest Landau level, this is not the 
case for the Ising QHF. In that case, spin up and spin down states come from different Landau levels. Thus, 
even when the magnetic impurity shuffles the spin up and spin down particles somehow in the Heisenberg 
QHF, the density does not change. 

Finally, we comment on densities in the inhomogeneous states in the Ising QHF. Results shown in Fig. 35b 
belong to quite small systems (12 particles at most). In the largest system studied, we observe a maximum 
in the density direct at the position of the impurity (x ~ 0) and the maximum approaches the value of 
density in a homogeneous system. With some imagination this allows for a hypothesis that - if domains 
are formed in an infinite system - the density will be inhomogeneous close to the domain boundary while 
remaining homogeneous inside a domain. However, we would have to study larger systems to confirm this 
speculation. 

2.3.4 The half-polarized states 

The inner structure of the half-polarized (S ~ -/Ve/4) ground state at filling 2/3 is investigated in this 
Subsection. We will argue that this state (assuming short-range interaction) resembles rather the incom- 
pressible singlet and polarized ground states at u = 2/3 than the Ising quantum Hall ferromagnet in the 
5*^ = sector as described in Subsection 2.3.3. 

In this Subsection, by half-polarized ground states we mean the 8- and 12-electron S = Ne/4 states GSs 
and GSi2 as introduced in Sec. 2.2, cf. correlation functions in Fig. 28. 

The ground state in a homogeneous system has a nearly constant density (oscillations due to the center-of- 
mass part wavefunction are less than 0.1% in the 12-electron system). This changes when a weak (5-line 
magnetic impurity along y is applied. Not only the polarization but also the density becomes inhomoge- 
neous. Fig. 38. The first minima of n{x) are at the same position ri w 2.2£o in the two system sizes 
considered and decaying oscillations are likely to follow at larger distances. Comparing the two system 
sizes in Fig. 38a, we find a much weaker response in the larger system, but this still does not have to imply 
a vanishing response in an infinite system, cf. discussion of the singlet state in Subsec. 2.3.1. 
Unlike the Ising quantum Hall ferromagnet discussed in Subsect. 2.3.3, the half-polarized states seem to 
have an intrinsic length scale in n{x) (of the order of ri). Fig. 38a. It is remarkable that this ri matches 



adp header will be provided by the publisher 



61 




-5-4-3-2-1 1 2 3 4 5 

x/ln 



1 

0.75 
0.5 
0.25 




Ml=0 
8/8, Ml=0.01 
\l2/12, Ml=0.01 



r 



X/Ir 



(a) Heisenberg QHF (v = 1). 
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Fig. 37: Different quantum Hall ferromagnets (QHF), half-polarized states: polarization response to a magnetic 5-line 
impurity in systems of 8 and 12 particles. 



quite well the position of the first maximum in the density of the LaughUn state (1/3) responding to an 
impurity. Fig. 33. 

Contrary to the density, the polarization does not show an intrinsic length scale as positions of the first 
minima in iVe = 8 and = 12 systems mismatch considerably. Fig. 38b. However, the polarization 
response here differs from the behaviour of the singlet state. Fig. 34b. Rather, Fig. 38b suggests that 
ni{x)/n{x) 0.75 as we go away from the impurity for the half-polarized states. This behaviour was 
classified as 'incompressible' in Fig. 32a. 

These observations suggest that the presence of an impurity will not lead to a splitting of the state into two 
domains (one with spin up, second with spin down), which we could expect for Ising QHF, Fig. 37. It 
seems that an impurity will rather change the polarization of the system only locally, in an 'incompress- 
ible manner'. Fig. 32a. The density response has the same characteristic length scale as the singlet and 
polarized ;/ = 2/3 ground states and such a length scale is absent in the polarization in agreement with 
behaviour of the singlet state. Fig. 34b. 

A state with Sz = iVp/4 comprises of jN^ electrons with spin down ('minority spins') and jN^ electrons 
with spin up ('majority spins'). Since the two populations are not balanced, we may gain extra information 
by speaking to them separately. The simplest concept, assuming non-interacting electrons, would be that 
Hmi,i, Hmi.i and Hmi (55) give rise to responses in ratio i : | : 1. Very roughly, this is indeed the case. 
Heights of the central peak {x = 0) for these three types of inhomogeneities are indeed approximately 
in this ratio, both for the density and for the polarization. Fig. 39. In the following we will discuss 
investigations with spin-dependent perturbations in more detail. 

Let us separate the density of majority and minority spins. Fig. 40. We will argue that the half-polarized 
state with electrons consists of two coexisting and weakly interacting liquids: Ne/2 electrons in a fully 
polarized liquid (with = ^e/4) and Ne/2 electrons in a S*" = state. Minority spins are thus present 
only in the S*" = liquid whereas majority spins occur in both of them. Concentrate on Fig. 40c. 

(i) Minority spins (|) react almost equally to Hmi.] and ~Hmi,i- They reflect only changes in the 
~ Q liquid and there are as many up as down spins in it. In fact, the Hmi,^ impurity influences also the 
polarized liquid component, but we cannot see it in the density of minority spins provided the two liquids 
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Fig. 38: Half-polarized ground state {S — Ne/4:) responding to a 5-line magnetic impurity, Hmi,^ (55). 
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Fig. 39: Half-polarized ground state (S = Ne/4:) responding to a (5-line magnetic impurity. Different impurity types 
are considered: Hmi, Hmi^] and Hmi,i (55). 



do not interact appreciably. The combined effect of H^j ij — Hmi,i causes a response of about the sum 
of these two. 

(ii) Majority spins (f) react differently to Hmi,] and —Hmi.i- We should keep in mind that rt] reflects 
changes in both (polarized and S" = 0) liquids. The latter impurity inflicts changes only on the S"" = 
part, whereas the former impurity acts on both liquids. If both liquids would have the same sensitivity to 
the considered impurities, we could expect responses in ratio 4:3:1 (Hmi to Hmi.] to Hmi.])- The fact 
that responses observed in Fig. 40c (measured by the height of the central maximum) are in ratio 3:2:1 
could be an indication that the polarized liquid is less sensitive than the S"" = liquid. 
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(a) 8 electrons, repulsive. (b) 8 electrons, attractive. (c) 12 electrons, attractive. 

Fig. 40: The same as Fig. 39, but the density is decomposed into the density of majority spins (n^) and minority spins 
(nj.). By attractive (repulsive) is meant that the 5-impurity at x = is attractive (repulsive) for the majority spin. 

(ill) Note also, that responses are the same (up to an inversion) for attractive and repulsive impurities, Fig. 
40a and 40b, provided the impurities are weak. 

Studies of the eight electron half-polarized state. Fig. 40b is not in conflict with this interpretation, re- 
sponses in densities are quantitatively different though. However, we should be cautious in drawing strong 
conclusions as these systems with primitive cell of size 12 flux quanta correspond to the smallest system 
{Ne = 4) considered in Fig. 33 (i^ = 1/3 state plus an impurity) and in that case finite size effects are 
already very strongly pronounced. Thus, the twelve electron system can be considered as the smallest 
system with finite size effects not playing a major role. 

Conclusion 

The hypothesis of the coexistence of the spin singlet and polarized liquids in the half-polarized states 
(HPS) seems to be supported. We have pointed out some similarities between the HPS and the former two 
incompressible states. In contrast, response to magnetic impurities seems to be different for the HPS and 
the Ising quantum Hall ferromagnet (in the Sz — sector) which would be the direct counterpart of the 
HPS if composite fermions are substituted by electrons. 

In general, it is not very surprising that electronic systems (i^ = 2 Ising QHF) differ strongly from the CF- 
counterparts. We have already seen this in correlation functions in Subsec. 2.1.1. However, the observed 
differences seem to be too deep to allow us to establish a relation between QHF states and the half -polarized 
states introduced in Sec. 2.2. 

2.4 Deforming the elementary cell 

In this Section we discuss another way of how to investigate fractional quantum Hall states. We will exactly 
diagonalize = 1/3 and u = 2/3 systems in elongated rectangular elementary cells. The dimensions are 
a by b, the aspect ratio is thus a : b > 1. The area of the rectangle is always kept constant, ab = 
(1), and therefore 

ab=2TTllN„,, =^ a = in\j2'KN,n\ , b = 4v/27riV„j/A , \ = a:b. (56) 

What can we expect? In the first approximation, we would say (i) nothing happens for an isotropic state 
such as the = 1/3 Laughlin liquid and (ii) crystalline or wave-like states will change both in energy 
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(a) Coulomb interaction (1/3). 



(b) Short-range interaction (2/3). Blue 
crosses mark the energy of a full Landau 
level. See the comment [3] for details. 



(c) Polarized, singlet and 
half-polarized states aX u = 
2/3. Overview. 



Fig. 41: Spin polarized eight electrons at filling factor 1/3 and 2/3. Energy of the lowest states versus aspect ratio of 
the primitive cell. An overview of the polarized, singlet and half-polarized states in the same scale is presented in the 
last panel. 



and in density. The reason is, that structures in homogeneous liquid states (as for example in correlation 
functions in Subsec. 2.1.1) are intrinsic and not incurred by the finite system size. Consequently, we 
expect the liquid state to change neither their energy nor their correlation functions, at least not on short 
distances, if a : 6 is slightly varied. On the other hand, an integer multiple of the period of a wave-like or 
crystalline state must be necessarily equal to a and/or b, hence by varying the aspect ratio we force it to 
change its period. In a classical crystal this means compression, or better deformation, since total 'volume' 
ab remains constant, and we expect it to cost energy. 

This investigation of = 2/3 systems was partly motivated by the work of Rezayi et al. [63] who 
investigated one particular type integer quantum Hall ferromagnet. Their exact diagonalization on a torus 
showed an iV„i-fold nearly degenerate ground state and the authors argued that these states comprised of 
stripes of alternating spin polarization (Subsec. 2.3.3) oriented parallel to one side of the rectangle, for 
example a. As they varied the aspect ratio, the states still remained degenerate, and their energy E{X) 
changed proportional to b. In fact, the degeneracy even improved: the small energy differences between 
the Nm states dropped. This was a strong argument for the stripe order, since then dE{\)/db can be 
interpreted as energy per unit length of an interface between a spin up and spin down stripe. 

2.4.1 Incompressible ground states 

As usual, we will start with v ~ 1/3, being probably the best understood system. This will also be the 
only case where we will discuss Coulomb interacting systems, in the rest we will stay with short-range 
interacting systems. 

Coulomb versus short-range interaction: u = 1/3 

The spectrum of a Coulomb-interacting system has a quite rich structure. Fig. 41a. The ground state energy 
exhibits several minima as a function of the aspect ratio of the elementary cell. In fact even more structure 
seems to appear in larger systems, as far as it could be inferred from comparing 6, 8 and 10 electron 
systems. In the following, we will show that this structure occurs mainly due to the long-range part of the 
Coulomb potential, it should be possible to describe it mainly by the Hartree part of the total energy or 
simply that it is due to formation of charge density waves (CDW) resembling Wigner crystals. Differences 
between Wigner crystals and CDWs are discussed below. In fact, energy of the states in question. Fig. 
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(a) The 'square' crystal state (aspect ratio 2). 
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(b) The hexagonal crystal state, aspect ratio 4/\/3 ~ 2.31. 

Fig. 42: Charge density waves resembling Wigner crystal states are among the lowest excitations in a Coulomb in- 
teracting 1 /3 system. Their energy is minimized (as a function of aspect ratio) when the elementary cell matches the 
crystal geometry. Correlation functions in eight-electron systems are shown, length of x- and i/-sides corresponds to 
the particular aspect ratio. 



42, will contain strong exchange contributions. Nevertheless, these states are very similar to the classical 
states which minimize the Coulomb energy. In a second step, we will discuss how correlations (and energy 
due to correlations) depend on the aspect ratio. Fig. 41b. 

In order to understand the the aspect-ratio dependence of energy of the Coulomb-interacting ground state, 
let us focus on two low excited states marked by arrows in Fig. 41a. These two states are just the CDWs 
mentioned above and they look almost like Wigner crystals: one hexagonal, another square, as density- 
density correlation shows. Fig. 42. It is then not surprising that the energy of such states is minimal, when 
the aspect ratio matches its geometry. For eight electrons considered here, this happens for^ 4rf : 2(i = 2 and 
4<i : (2V3/2d) = 4\/3/3 for the square and hexagonal crystal, respectively. Perhaps the most apparent 
difference between a CDW and an (unpinned) Wigner crystal is that for the latter state we expect the 
correlation function to drop almost to zero between the 'lattice sites' and this is not the case here. Fig. 
42. The reason is that at filling factor v ~ 1/3, the system is too densely populated, or mean interparticle 
distance is too small, rmean/^o = ■\/27r/i/ « 4.35 (1) to allow the electron density (or correlation function) 
to vanish between two sites. Remember that an electron within the lowest Landau level cannot be localized 

^ d is the 'lattice constant'. 
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(a) Aspect ratio 1.00. 



(b) Aspect ratio 2.00. 



(c) Aspect ratio 2^3 ^ 3.46. 



Fig. 43: Evolution of the Laughlin state with aspect ratio of the elementary cell (Coulomb interaction). Correlation 
functions are shown. 



more strongly than on a length scale of the order of unity (magnetic length £q). Even if we assembled a 
hexagonal Wigner crystal at = 1 /3, the wavefunctions at neighbouring sites would strongly overlap and 
it is then more favourable for the electrons to retain some features of the Laughlin correlations. As a result 
we obtain a CDW (or a 'strongly correlated crystal' [47]) like the state in Fig. 42b. At lower filling factors, 
rmean/^o is larger and Wigner crystal states become possible. This can be interpreted as a quantum phase 
transition from liquid to solid as the filling factor is decreased and the extensive studies in this field suggest 
the critical value ~ see Sec. 5.7 in Chakraborty [14] for a review. 

The ground state (GS) energy reflects these geometrical conditions. This state also minimizes its energy 
when the square crystal can easily be formed, but at short distances it strictly preserves the liquid-like 
correlations. Fig. 43. It is isotropic at short distances, in Fig. 43b, the ring corresponding to the first 
maximum is circular and not deformed into an oval for instance. Fig. 42, and also g{r) cx (not obvious 
in Fig. 43). It seems plausible that the increase of GS energy around a : w 3, Fig. 41a, is due to the loss 
of isotropy at shorter distances. The ring of the first maximum in g{r) disappears. Fig. 43c, the feature 
g{r) (X however remains. It is important to know, that unlike the energy, the structure (correlation 
functions) of the ground state is quite insensitive to the type of interaction (Coulomb or short-range). 
Let us now proceed to short-range interacting states. The crystalline states disappear from the realm of 
low-energy excitations. The ground state energy is completely independent on aspect ratio, it is zero. Fig. 
41b. Energies of fully polarized 2/3 states displayed therein are equal to those of 1/3-systems up to a 
constant shift. This constant depends on aspect ratio, but the dependence is imperceptible up to a : 6 « 4 
(for 4 electron system). In fact, the ground state is rigid in the following sense: a state can have zero energy 
only if there are three zeroes on the position of each electron in the wavefunction. The wavefunction is 
completely determined by this condition together with the confinement to the lowest Landau level. It is 
even surprising, that given how strictly determined the ground state is, it once resembles a liquid (for 
a : & w 1) and another time a CDW state (larger aspect ratios). Fig. 43c. 

Assuming fully spin polarized electrons, 2/3 and 1/3 systems (e.g. 8/12 and 4/12) are particle-hole conju- 
gated. Thus, spectra of these systems are identical up to a constant energy shift, which is just the Coulomb 
(or short-range interaction) energy of a completely filled lowest Landau level. Note that this energy varies 
with aspect ratio (both for Coulomb and for short-range interaction). The common statement that interac- 
tion energy of a full LL is a constant is valid in a broad range of aspect ratios, but not everywhere. In Fig. 
41b, this holds up to a : 6 < 4, Subsect. 1.5.4. This is shown in Fig. 41b. The = 1/3 Laughlin state 
has zero energy for any aspect ratio (not shown), the 2/3 ground state energy is then just the Hartree-Fock 
energy of a completely filled Landau level. Beyond a : 6 w 4, this energy is no longer constant, indicating 
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(a) Eight electrons. 



(b) Ten electrons. 



Fig. 44: Low lying states at filling 2/3 under vanishing Zeeman splitting versus aspect ratio of the primitive cell. Note 
that energy of the singlet ground state remains about constant for aspect ratios < 1.4 in the larger system. 



that the deformation of the elementary cell becomes pathologic and beyond this point (at latest), the model 
no longer describes a 2D system but rather an effective ID system. 

Consider a/b 3> 1. Then the A^e electrons are located on a very thin cylinder [62] of length oc \fafb 
(area of the cylinder is fixed by filling factor, ah = 2TiN,n) and single electron states resemble 'rings on a 
pole'. The mean distance between electrons is then oc ^Ja/b/N^ and Coulomb energy is then proportional 
to The increase of the ground state energy for very large aspect ratios. Fig. 41a, is due to the 

repulsion between an electron and its own periodic image in the 'short-direction' . 

Excited states are sensitive to deformation of the elementary cell even more. Energy levels group into 
branches beyond a/b « 2, Fig. 41b. Keeping in mind the transition towards an effective ID model, 
these branches can be attributed to 0, 1, 2, etc. pairs of 'rings on a pole' sitting at neighbouring sites. In 
illustrative terms, there is no longer enough room for two electrons to be positioned in 'vertical' direction 
(along y axis, i.e. the shorter side of the elementary cell) except when they freeze into a crystal. 
In conclusion, going beyond aspect ratio ^ 2 (in a N,n = 12 system) the system cannot be taken as a 
faithful model for an isotropic infinite system. 

The singlet state 

The singlet ground state is apparently more sensitive to varying the aspect ratio. Its energy changes at much 
smaller deformation than that of the polarized state. Fig. 41c. However, comparison between systems of 
different sizes shows that its energy is also constant, provided that the aspect ratio is not much larger than 
one and the system is large enough. Fig. 44. This is another hint at isotropy of the state. A crystalline state 
responds more strongly to a change of a/b, since this is in principle an attempt to compress the lattice in 
one direction while expanding it in the other direction. Recall just the CDW states in = 1/3 systems 
marked by arrows in Fig. 41a. 

This is in agreement with a direct observation of correlation functions. Fig. 45. In particular, the ring 
structure in g^i (r) (or maximum at rg ~ 3.4£o) remains preserved even for aspect ratios a : 6 « 3, Fig. 45 
right. This is similar to how the ring structure of the first maximum was preserved in the deformed = 1/3 
Laughlin state. Fig. 43b. Also, looking at .gtT (•^) ™d .9Ti (^) the deformed singlet state, the sum of these 
two seems to remain constant beyond ro even in deformed systems, in spite of 9|j(a;) decreasing beyond 
X = tq. This was just the conclusion in a : = 1 systems. Fig. 15, and it suggests that the singlet state 
did not change much even in a quite strongly deformed system {a : b < 3). Moreover, this finding allows 
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Fig. 45: The singlet 2/3 state in elementary cells of different aspect ratios: density-density correlation between unlike 
spins (the two aspect ratios shown in 2D plots are a : b = 1 and 2). 
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Fig. 46: Half-polarized states (S = at filling 2/3 versus aspect ratio. States with next larger spin are well above 

(out of scale here). 



us to use deformed systems to study what happens on slightly larger distances than distances accessible in 
a square cell of a fixed area because maximum distance between two electrons in a deformed elementary 
cell, ^£o^/2^^Nm{X + 1/A), grows with increasing A. 

Regarding the energy, which seems to react more sensitively to deformations than the correlation functions, 
the following speculation seems plausible. If the singlet state is a liquid of f — | pairs of characteristic 
size ro ~ 3.4^0. Subsec. 2.1.1, it ought to be more sensitive to aspect ratio variations than the Laughlin 
state just because such a pair in the jy = 2/3 singlet state is larger than a single electron in the = 1/3 
Laughlin state. 



2.4.2 Half-polarized states 

The half-polarized states can be expected to suffer severely under the finite size of the system. A system 
with eight electrons contains only two electrons with minority spin. Contrary to fully polarized systems 
(where eight particles is already fair enough), it is thus the smallest system with S = Ne/4: where many- 
body effects can be studied. 

Let us compare how systems of two different sizes respond to varying aspect ratio. In an eight-electron 
system. Fig. 46a, there are four low lying states: the ground state at a/b = 1 with = (tt, tt), a (0, 0) 
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state which becomes the ground state at a/6 > 1.5 and a pair of degenerate states, (0, tt) and (tt, 0) [k"^ 
is defined in Subsec. 1.5.2). The former two states are isotropic (and lie in high symmetry points of 
the Brillouin zone), the other two are spin-density waves in x and y direction, judging by the correlation 
functions (not shown). Moving away from aspect ratio one, degeneracy of the latter two is lifted - just as 
the 90 deg rotational symmetry of the elementary cell is broken - and the wave along x (the longer side) 
becomes energetically more favourable. It is quite conspicious that this state evolves parallel to the (tt, tt) 
state for aspect ratios above w 1.4. For these values of a : 6, the inner structure of these two states seems 
to be very similar, too. 

In the low-energy sector, a (0, 0) state is absent in a 12-electron system. Fig. 46b. In other respects, 
however, the situation is quite similar to the smaller system. There is a well-separated (tt, tt) ground state 
in a square cell and this state becomes nearly degenerate with a (0, tt) state for aspect ratios > 1.4. Also, 
the energy of these two states decreases with increasing aspect ratio and eventually reaches its minimum. 
In contrast to the smaller system, the minimum occurs later, at a : & w 2.4 (Fig. 46b) compared to w 1.6 
in Fig. 46a, but this occurs also for the incompressible states, e.g. the singlet at i/ = 2/3 (Fig. 44a vs Fig. 
44b). The correlation functions of these two states, (0, tt) and (tt, tt), are similar to those of the (0, tt) and 
(tt, tt) states in the eight-electron system (not shown). 

Now turn to the correlation functions of the (0, tt) and (tt, tt) states in a 12-electron system, Fig. 47. 
Both states are quite isotropic, for a square elementary cell, first at a very close look, we find a slight x 
versus y anisotropy in the (0, tt) state. However, already under slight variation of the aspect ratio, stripe 
structures parallel to the shorter side evolve {a : b ~ 1.2, Fig. 47a). In this respect, both states look quite 
similar. Fig. 47 (or compare Figs. 47a and 47b), and we should stress that the differences in the correlation 
functions between the isotropic (at a : = 1) and the wave-like state {a : b ~ 1.8) are very large, both in 
isotropy/anisotropy and in the short-range behaviour. Fig. 47c. This is in a strong contrast to the behaviour 
of the incompressible states, e.g. the Laughlin state which preserves lot of its original isotropy even at 
a : & « 2, Fig. 43. 

These observations suggest the following interpretation. The half -polarized ground state at = 2/3 is an 
isotropic state which however inclines to the formation of a spin-density wave. The wave has the shortest 
period allowed by the number of electrons, i.e. it resembles an antiferromagnetic ordering (titi • ■ • rather 
than ttii ■ • ■> for instance) as the correlation functions in the rightmost column in Fig. 47a suggest. 
Since there are just three j-electrons in the system, we expect two stripes (the third j-electron is just at the 
origin) in gii{x, 0) in the case of titi • ■ • ordering. In more detail, see Fig. 47d: the minima/maxima 
in g^iix, 0) match well with the maximum/minima in g^i{x, 0). In other words, spin up is followed by 
spin down. However, the amplitude of oscillations in {x, 0) is moderate. Fig. 47c, and hence we should 
rather classify the state as a 'spin density wave' than e.g. a state with stripe domains of alternating spin 
polarization. 

On the basis of the present investigation, it is not clear whether in a large enough system, this spin wave 
state is the ground state, a low-energy excitation or it is degenerate with the isotropic ground state. Even 
though the GS at a : 6 > 1 (spin wave) has a lower energy than the isotropic state at a : 6 = 1, Fig. 46b, 
this does not say much about which state would be the ground state in a larger system. We saw a similar 
situation for the v ~ 1/3 Laughlin state. Fig. 41a, or the singlet ly = 2/3 state. Fig. 44. The energy of the 
ground state was not at its minimum at a : b ~ 1, yet the isotropic {a : b ~ 1) state is probably the real 
ground state in the thermodynamic limit. The question how to decide which state - isotropic or anisotropic 
- will be preferred in infinite systems remains open, but comparison between systems of more different 
sizes could be very helpful. 

2.4.3 Conclusions 

It has been demonstrated that isotropic states like the fully polarized or singlet incompressible = 2/3 
ones tend to be insensitive to slight deformations. The response was observed in the energy of the state 
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Fig. 47: Evolution of two half-polarized states lowest in energy with growing aspect ratio of the elementary cell (12 
electrons, v — 2/3, short-range interaction). Correlation functions are shown. 



and in its correlation functions, where we saw that especially the short-range behaviour remains basically 
unchanged. The insensitivity improves with increasing the system size (number of particles). We also 
registered some differences between the singlet and polarized state. In systems of equal size (area) the 
former state was disturbed by smaller deformations. This agrees with our previously mentioned hypothesis 
(Subsec. 2.1) that the singlet ground state consists of pairs of electrons with unlike spin. Since the typical 
size of such a pair was rather large (3.4^o), the singlet state will suffer under the finite size of the system 
more than the polarized state where the 'relevant particles' are still electrons whose size is about £q. Imag- 
ine filling a container once with ten tennis balls (~ polarized state) or with five footballs singlet state). 
Slightly deforming the container will probably affect the latter system stronger. 

Investigation of the half -polarized state revealed that while the state is isotropic in a square cell, it tends to 
build a unidirectional spin density wave for aspect ratios not far from one. In this regime, it also becomes 
degenerate with one other state. Correlation functions of the both states (in deformed elementary cells) are 
quite similar to each other We suggested that these states have an antiferromagnetic ordering in agreement 
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with both correlation functions and wavevectors of these two states, = (0, tt) and (tt, tt). The question 
which state (isotropic or spin wave) is the real ground state in an infinite system remained unanswered. 

2.5 Summary and comparison to other studies 

2.5. 1 The incompressible states: the polarized and the singlet ones 

We studied various properties of the fractional quantum Hall states with spin degree of freedom at filling 
factors 1/3, 2/3 and 2/5: correlation functions, response to magnetic and non-magnetic (5-line impurities 
or to deformation of the elementary cell. Briefly summarized: 

(i) The results are in agreement with the concept of incompressibility of these states and also (in the case 
of = 1/3) with some earlier studies, e.g. [86]. 

(ii) Even though these states can be imagined as composite fermion systems with integer filling, the analogy 
to Landau levels completely filled with electrons can often be misleading. For instance electrons of unlike 
spin are strongly correlated in the v ^ 2/3 singlet state while they are completely uncorrected in a = 2 
singlet state. 

(iii) We inferred pairing of spin up and spin down electrons in the 1/^2/3 singlet state. In the spin- 
unresolved density-density correlations, this state looks as if the two electrons in each pair were located 
exactly at the same position and the pairs then formed a ly = 1 state. This conclusion was not possible 
for the 1/ = 2/5 singlet state thereby highlighting differences between fillings 2/5 and 2/3 which are very 
closely related within composite fermion theories. 

2.5.2 Half-polarized states 

We identified a highly symmetric half-polarized state at filling factor 2/3 which comW become the absolute 
ground state in a narrow range of Zeeman energies (or magnetic fields). Such a state is completely unex- 
pected in mean-field composite fermion theories. Extending earlier studies with exact diagonalization on 
a sphere we showed that extrapolating the energy of this state from finite size exact diagonalizations to the 
thermodynamic limit is problematic and the question whether the half-polarized state really becomes the 
absolute ground state remains open. 

Investigations on this state both for short-range and Coulomb interacting systems showed strong similar- 
ities to the incompressible singlet and polarized states at = 2/3. Consequently, we suggested that the 
singlet and the polarized state coexist within the half-polarized state. The state might be gapped for short- 
range interacting electrons but even if yes, it is probably not gapped for Coulomb interacting systems. 
These differences in spectra accentuate the fact that extrapolations to infinite systems should be taken with 
extreme caution. It also means, that the definition of the short-range interaction should be reconsidered. 
The model may be an oversimplified if we study the half-polarized states since the mean distance between 
two minority spin electrons is rather large, higher pseudopotentials should also be taken into account. 
The half-polarized state forms a pronounced spin-density wave, or antiferromagnetic order, when anisotropy 
is introduced from outside (deformation of the elementary cell) but we could not conclude whether this 
spin-wave will be more energetically favourable than the isotropic form in much larger systems. 

2.5.3 Half-polarized states: other studies 

Let us first briefly recall other suggestions which appeared since Kukushkin et al. presented their exper- 
iment showing a plateau of the polarization at the value of one half. All works mentioned below can be 
apphed both to filling factor 2/3 and 2/5 in principle. Unless necessary, we will not distinguish between 
these two cases. 



72 



K. Vybomy: Spin in fractional quantum Hall systems 



Ganpathy Murthy [54] was attracted by the idea that correlations favour either the spin singlet or the fully 
polarized state. At the point where the two ground states cross (recall Figures 24 and 9a), electrons could 
prefer to form a translationally non-invariant state consisting of regularly alternating areas of (locally) 
singlet and (locally) polarized states arranged into a partially polarized density wave (PPDW). He argues 
that this structure ought to have square rather than a hexagonal symmetry. The energy of the PPDW state is 
evaluated within the Hamiltonian theory of composite fermions [56] and it is shown that the PPDW state is 
stable (against one-particle excitations) and lower in energy than the (homogeneous) singlet and polarized 
states. The period of the density wave should be 2^Jtt£* (19) which is 7.93^o for filling 2/5 and 6.14^o for 
2/3. Charge modulation in the wave should be quite weak (in the order of 1%). 

Apal'kov, Chakraborty, Niemela and Pietilainen [10] object that the energy of the PPDW is too high and 
claim that a homogeneous Halperin state in the two crossing CF Landau levels (see below) should have 
a lower energy. Without invalidating the following results, this estimation seems to be however incorrect 
[55]. As the mentioned Halperin state cannot account for the half-polarized states, Apal'kov et al. suggest 
another candidate for the half-polarized state, a non-symmetric excitonic liquid. They consider only the 
'active levels' meaning the two CF Landau levels which cross. These {two) levels have total filling of 
one, i.e., there are only Nm electrons for Nm places in the t level and Nm places in the [ level. By 
convention, they define a f-particle as an 'electron' and a missing J, -particle as a 'hole'; an 'electron'- 
'hole' pair is an 'exciton' and a pair of a ^-particle and a missing f-particle is 'vacuum'. Owing to the 
constraint N-\ + = N„i, one-particle states can be mapped onto a system consisting solely of 'vacua' 
and 'excitons'. The partial filling factor Ni/Nm € [0; 1] then gives simultaneously the polarization and 
the number of 'excitons' (by iV,„)- Note, that 'excitons' are bosons by virtue of an integer spin. 
From this viewpoint, the 1^ = 1 quantum Hall ferromagnet (being described by the Halperin (1,1,1) state) 
is a Bose condensate of excitons. In that case, all the excitons are non-interacting and have zero angular 
momentum L. This is most easily seen by the fact that .g|| (0) = 0. On an 'electron' (f particle), there is no 
I particle, i.e. there is a 'hole'. In an exciton (hydrogen atom), the only wavefunctions with tp{r = 0) ^ 
are those with L = 0. On the other hand, g^i{0) ~ follows from the fact that the Halperin state 
has maximum polarization and thus the spatial part of the wavefunction must be totally antisymmetric. 
Apal'kov et al. suggest that the half -polarized state at i/ = 2/3 or 2/5 could be a condensate of excitons 
with L = 1 for which they call it nonsytnmetric. 

To support this idea, they perform exact diagonalizations in a = 1 system with several model interactions 
which are meant to describe the two - active - crossing CF Landau levels. These interactions are derived 
from the Coulomb potential with suppressed short-range component, probably (without justification) with 
the intention to describe interacting composite fermions. Stability of the half-polarized state is substanti- 
ated by showing that the energy versus polarization curve has a downward cusp at half-polarization. On 
the other hand, g^i{0) 7^ in the half -polarized state indicates that the 'excitons' do not have L = 0. The 
particular value of L = 1 is demonstrated by other means. 

Finally, the idea of Eros Mariani [51] should be presented. Parallel to the previous two works, the two 
'active' crossing CF Landau levels are considered. An assumption is made that they both have a partial 
filling of 1/2 rendering (after a second Chern-Simons transformation) two Fermi seas of 'free' composite 
fermions (of second generation). Mariani et al. show that interaction of these objects with fluctuations 
of the gauge field leads to an attractive effective interaction between particles with opposite spin and 
momentum. In analogy to superconductive pairing, this implies a gapped ground state. An estimation of 
the gap is given. 
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2.5.4 What are the half-polarized states then? 

Presently, it is not clear which (if any) of the candidates proposed in the previous subsection describes 
the half-polarized reality. As Murthy correctly mentions, the final answer should be given by exact diag- 
onalization of a large enough system. Unfortunatelly, we dispose of systems not larger than 12 particles. 
Nonetheless let us compare the candidates with what was presented earlier in this chapter. 
The downward cusp in energy-versus-polarization dependence cannot be assured by the calculations pre- 
sented here. However, if the lowest half-polarized state indeed becomes the absolute ground state at the 
transition between the singlet and polarized state (see extrapolations in Fig. 25), the cusp is likely to be 
present. In the other case, it will turn into an upward cusp, as the calculated spectra suggest. 

Results presented here indicate, that the half -polarized ground state (2/3) has (k^, ky) ~ (tt, tt) and that 
it shows similarities to the singlet and polarized ground states (Fig. 27). In particular .gtx(O) « 0, which 
is in contrast with the model of a nonsymmetric exciton liquid (cf. the correlation functions in [10]). 
Comparison between short-range interaction and Coulomb half-polarized states (Fig. 2.2.6) suggest that, 
similar to the Laughlin state, the short-range part of the interaction plays the major role. From this point, 
the model discussed by Apal'kov et al. [10] seems to be more appropriate rather for some other systems. 
Positioning of the half-polarized state out of the centre of the Brillouin zone could be an indication that 
it is indeed a standing wave. This is also supported by spectral properties when the elementary cell is 
deformed. Fig. 46. The two lowest states becoming degenerate at aspect ratios larger than 1.4 could be 
a charge/spin-density wave (note also the correlation functions. Fig. 47). The fact, that the energy of the 
ground state lowers with increasing aspect ratio could indicate that this state is more stable than an isotropic 
one. However, caution is advised here, since the singlet incompressible ground state does the same. Fig. 
44, while its isotropic form is the true ground state. 

Theory of the 'superconductive' pairing was not addressed so far Comparisons on the level of correlation 
functions, possibly in /c ''-space, are in principle possible, but quite complicated because of the two Chern- 
Simons transformations involved. 



3 Quantum Hall Ferromagnetism at z/ = 2/3? 

Like the previous Chapter, this Chapter also starts from the fact that there are two distinct ground states at 
filling factor 2/3: a spin-singlet and a fully polarized one. Their structure was studied in Chapter 2 and 
we also interpreted them in terms of composite fermions. Fig. 9a. Whichever of these two becomes the 
absolute ground state depends on the Zeeman splitting which favours spins aligned parallel to magnetic 
field. The singlet state is the lowest in energy for vanishing Zeeman splitting. However, increasing the 
Zeeman splitting, its energy remains unchanged while the energy of the fully spin polarized state decreases 
and eventually this other state becomes the absolute ground state. This simplest scenario, sweeping the 
Zeeman energy while magnetic field is kept constant, is not very usual, albeit it is experimentally possible 
[50]. However, even if we simply sweep the magnetic field (and keep constant filling = 2/3 which 
requires a simultaneous change of the electron density), the Coulomb energy of the singlet state changes 
c>c ^fB and that is slower than the Zeeman energy of the polarized state in the limit of large B. Therefore, 
the qualitative discussion above is still valid. The total energy bilance of the two ground states (in SI units) 
is thus 

polarized: Ep{B) 
singlet: Es{B) 



- giiBN.B = -\ap\VB - \D.p\B , 



|CJx/5, 
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Fig. 48: Energies of low lying states at v = 2/3 in a homogeneous Coulomb-interacting system with Zeeman field: 
transition from an incompressible singlet ground state to a fully polarized incompressible ground state. Different 
numbers of particles in a square with periodic boundary conditions are considered, the scenario is however the same 
in all cases. 



where is the number of particles and > are the total Coulomb energies in units e'^ /Aireio (as 
calculated by exact diagonalization, for example; not per particle). Obviously, Ep{B) < Es{B) for B 
large enough. What the critical field Be is, where both energies are equal, depends on {E'^ — E^)/Ne- 
This quantity is accessible only numerically and depends on Ne although we may hope that it stays nearly 
constant for large Ne- 

Figure 48 demonstrates this singlet to polarized transition for 4, 6, 8 and 10 Coulomb-interacting electrons 
on a torus. Note that the energy units in Fig. 48, e^/ (47re£o) oc VB, change with magnetic field. In these 
units, the potential (Coulomb) energy of all states stays constant (singlet state) and Zeeman energy scales 

as cx e2/(47re4) ■ 

A close look at Fig. 48 shows that the magnetic field Be, at which the ground state transition takes place, 
varies non-monotonically. However, an extrapolation of energies of the two ground states to 
allows for a rough estimate of Be ~ 7 T in an infinite system, Subsect. 2.2. This is in quite good 
agreement with experiments [45], even though in some samples Be as low as « 2 T was observed [46, 70]. 
This could be due to deviations from an ideal 2D system, Subsect. 3.1. 

In the following we want to show that the existence of the spin structures and the formation of domains are 
of central interest for the understanding of the ground state transitions. 

The ground state is always either a singlet or fully polarized in a homogeneous system. The energies 
of these two states are equal at the transition. This is similar to an Ising ferromagnet, if we label the 
polarized state by pseudospin up and the singlet state by pseudospin down. In an infinite system at non- 
zero temperature, the Ising ferromagnet prefers a state with domains, some with (pseudo)spin up, others 
with spin down, to the two homogeneous states. First because entropy of the former is higher [41] and 
second because the total magnetization of a domain state is approximately zero while, locally, most spins 
are parallel to their neighbours thereby minimizing the energy of magnetic stray fields [11]. None of these 
two mechanisms was included in the studied model of a = 2/3 system. Nevertheless, the question 
was addressed how the system responds if such a domain-inducing mechanism is modeled by a magnetic 
inhomogeneity. Will the ground state split into regions of different spin polarization? With this question in 
mind, the inhomogeneity should prefer the singlet ground state in one part and the polarized ground state 
in another part of the system. 

From the experimental side, there are quite strong hints at ferromagnetism, mentioned in the introduction. 
Hysteresis, saturation (in time) of magnetoresistance, Barkhausen jumps etc. hint at ferromagnetic states 
with domain structure near the transition point. Intention of the present study is to support this interpreta- 
tion. 
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(a) Homogeneous system. 
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(b) System with 'domain-enforcing' inliomogeneity. 



Fig. 49: An idea of how to enforce domains at the crossing of singlet and polarized ground states of = 2/3. An 
average Zeeman field is chosen so that the both homogeneous states have the same energy. Modulation of the Zeeman 
field prefers the singlet state 'in the middle' and the polarized state 'at the edges' (note however the periodic boundary 
conditions). 

3.1 Attempting to enforce domains by applying a suitable magnetic inhomogeneity 

This and the following sections will be concerned with various attempts to induce the formation of domains 
close to the transition point. At the beginning we must discuss (i) how to enforce domains, what to add 
to the Hamiltonian, which form of inhomogeneity and (ii) how to detect them, which quantities should be 
observed. 

3.1.1 First attempt: the simplest scenario 

The simplest scenario is sketched in Fig. 49. In the homogeneous case, the Hamiltonian consists of two 
terms 



Coul 



H. 



Zeeman 



— T- 



1 



(57) 



i<j 



the Coulomb interaction and the Zeeman term. If the Coulomb energy is fixed, the energies of the two 
incompressible ground states can be shifted with respect to each other by varying the Zeeman term. If B 
is fixed at B = Be (i.e. the two ground state have the same energy), the Zeeman energy can be still varied 
by means of the g factor. Decreasing g slightly, the singlet state will become the absolute ground state, 
increasing g the polarized state will prevail. 

The idea of a 'domain-enforcing' inhomogeneity is to the turn the constant g into g{xj) = 5o + dii^j) in 
(57) and g{x) > go in one part of the system whereas g{x) < go in another. 

Or, speaking in terms of Fig. 48: we slightly modulate the magnetic field B and in one part of the system 
we consider B > B^ while in another i? < B^- By slightly we mean that only the spin degree of freedom 
is affected, not the orbital. This is an approximation. 
The full Hamiltonian to consider is thus 



H — Hcoul + H Zeeman + HmI , 

Hmi = ^gi{xj)fiBBcrl , {(pi\HMi\ipj) 



SijEMi 



0,1,. 



-1,... 
-1,... 



(58) 
1 

-1 
1 
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(a) An ovemew: spectrum and (inset) polar- 
izations of the singlet, transition and polarized 
ground state (left to right). 



(b) Polai'ization of the three states 
in detail (n = 10). Mean polar- 
ization subtracted. 



(c) The transition in sys- 
tems of different sizes. 



Fig. 50: Response of a ten-electron system to a weak (Emi = 0.002) magnetic inhomogeneity of the form given in 
(58). No sings of domain formation observed: the transition state does not respond stronger than the incompressible 
states. 



where \ipj) is a one-particle state localized around x = {j /Nm)a, cf. (21). This roughly corresponds to 
gi{x) having a 'rectangular wave' form {gi = 1 for < x < and |a < x < a and gi = —1 for 
\a<x< fa). 

The basic results of this model are: the ground states slightly change in accord with the inhomogeneity 
and nothing peculiar happens near the transition. As we sweep the magnetic field through B = Be, even 
in the presence of a weak inhomogeneity, the singlet state evolves 'smoothly and monotonously' into the 
polarized state, without any remarkable intermediate states. 

Typical results are shown in Fig. 50. A magnetic inhomogeneity (58) was applied to a ten-electron 
Coulomb-interacting system and its strength Emi was chosen to be ~ 10% of the incompressibility gap. 
Regarding the ground states and the gap, the spectrum remains virtually unchanged. Fig. 56a shows a 
comparison of the spectra between homogeneous and inhomogeneous systems. Looking now at the singlet 
and polarized ground states, we find a spatially varying spin polarization' n-\ {x) /n{x). Fig. 50a. However, 
the mean values of the polarization still remains at 0.5 (1) as it was in the homogeneous singlet (polarized) 
state. Fig. 50a, leftmost (rightmost) inset. The polarization of the 'transition state' has a mean value of 
0.75, i.e. just in the middle between the polarized and the singlet state. This is not surprising, since the 
'transition state' was taken to be a symmetric linear combination of the two crossing states (see Subsect. 
3.1.2). What is more interesting, is the variation of the polarization around the mean value. Fig. 50b: in 
this point, the 'transition state' lies just between the singlet and polarized states. Contrary to what we ob- 
serve in Fig. 50a (middle inset), formation of domains near the transition would mean that the polarization 
of the transition state should vary between 0.5 and 1. 

It could be that the system is simply too small for domains to evolve near the transition. On the other hand, 
this does not seem to be the case, since the response to the inhomogeneity does not grow with increasing 
system size but rather stays about the same. Fig. 50c. 

The particular parameters of the model presented in Fig. 50 could have been chosen unluckily so that 
domains could not evolve. Let us therefore discuss the inhomogeneous = 2/3 systems more thoroughly. 



Thi'oughout this Chapter, we will refer to p{x) = n| (x) /n{x) as to polarization. In the literature, another detinition is more 
common, P{x) = [n-^{x) — n^{x)\/n{x), both quantities are, however, equivalent; P(x) = 2p{x) — 1. 
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3.1.2 Turning crossing into anticrossing: inhomogeneous inplane field 

AiB = Be there is actually a crossing between the singlet and polarized ground states. Fig. 50a, rendering 
the transition jump-like just as in a homogeneous system. For the transition state (the middle curve in Fig. 
50b), we took a fifty-fifty linear combination of these two ground states. One could say, the transition 
occurs in an infinitesimally small interval of magnetic field around Be. 

In a realistic system, the transition is unlikely to happen all simultaneously in the whole system. Two 
mechanisms causing a more continuous transition in a finite interval of B are conceivable: 

(i) weak inhomogeneities: The spectrum (as a function of B) looks basically the same as in Fig. 50a, but 
there is an anticrossing at B w Be- 

(ii) strong inhomogeneities: The energy gap between the pair of the crossing ground states and the ex- 
cited states at B = Be (Fig. 50a) is reduced compared to the incompressibility gaps of the singlet and 
polarized ground states far away from Be, i.e. for B ^ and B ^ oo. Under influence of stronger in- 
homogeneities, it could be that some originally excited state (or more states) become ground state around 
B w Be while singlet and polarized incompressible states remain lowest in energy only far away from 
Be- If this turns out to be the case, it could be that more states (possibly of different Sz) can be mixed by 
the inhomogeneity, eventually rendering the ground state compressible. 

In this Subsection we will discuss the former possibility, the latter will be the topic of Subsect. 3.1.3. 
The ground state transition in Fig. 50a is a crossing even in the presence of the magnetic inhomogeneity 
Hmi (58) because the symmetry of Hmi is too high and it does not mix the two crossing ground states. In 
particular, [Hmi, S'^] ^ but [Hmi, S^] = and the singlet state \S) has Sz = Q whereas the polarized 
state \P) is Sz = Ne/2. Consequently Sz\S) = and Sz\P) = [Ne/2)\P), therefore {P\Hmi\S) = 
{Ne/2)-\P\SzHMi\S) = {Ne/2)-^{P\HMiSz\S) = 0. The inhomogeneity iJM/ mixes states with 
different S but only those with equal Sz- The fully polarized ground state (S = Ne/2, Sz = Ne/2) has 
also an S'^ = counterpart, since the homogeneous Hamiltonian commutes with spin lowering operator. 
This state, however, is a highly excited state at _B w Be, since its Zeeman energy is zero. 
Which terms added to the Hamiltonian can break this symmetry and what will then be the response of the 
ground state? 

Weak (inhomogeneous) inplane magnetic fields will have the desired effect. This scenario is not unlikely 
to occur in a realistic system. It merely means, that the extra fluctuating magnetic field which interacts 
only with the spins, is not pointing exactly in the direction of the (strong) external magnetic field causing 
the Landau level quantization. The existence of such symmetry-breaking inhomogeneities is very likely in 
realistic systems, although they might be very weak e.g. hyperfine interaction with nuclear spins. 
Let us consider a Hamiltonian with inplane magnetic inhomogeneities (IMI) of the form 

H = Heoul + H Zeeman + Hmj + HjMI , 

HiMi = ^go^J.BBxXxj)al , ((pi\HiMi\^j) 

3 

The main claim of this Subsection is that weak Hjmi only opens an anticrossing at the ground state 
transition. In other words, the relevant states still basically form a two-level system comprising of the 
(slightly disturbed) singlet and polarized ground states. The width of the anticrossing grows with increasing 
strength of the inplane field inhomogeneity, Ejmi- As to the width we refer either by the level splitting. 
Fig. 51b, or by the range of magnetic field where (Sz) noticeably changes. Fig. 56. 
This fact is best demonstrated in Fig. 5L Inhomogeneities are weak there (compared to both incompress- 
ibility gaps Eg Eg), i.e. Emi,Eimi ^ Eg, and the spectrum remains almost unchanged. Fig. 51a. 
The energy levels of the inhomogeneous system (points) are almost equal to the energies in a system free 
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(a) Spectrum. Solid lines: homogeneous system; 
points: inhomogeneities switched on. 



(b) Energy of the lowest two states at the anticrossing as IMI is being turned 
on. Right: level splitting 



Fig. 51: Response (of an eight-electron system) to a weak inhomogeneity in perpendicular and inplane direction (59). 
Perpendicular component Emi is the same as in Fig. 50. Exponent four comes from Ne/2, see explanations in the 
text. 
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(a) Weak IMI. (Eimi = 0.001) 



(b) Intermediate IMI. 
(Eimi = 0.005) 



(c) Crossing at different values 
of IMI. 



Fig. 52: Inplane magnetic inhomogeneity (IMI) turns the crossing between the singlet ground state and the polarized 
ground state, Fig. 51a, into an anticrossing. The cross-over between the two ground states can be observed either in 
the spectrum or in {S^) of the ground state as B is swept through Be- 



of impurities (lines). Only directly at i? w Be an anticrossing opens and the level separation AE grows 
with increasing Ejmi {Emi is kept constant). Fig. 51b. Even for a quite strong inplane inhomogeneity (of 
the order of Eg), the level splitting remains small (^ Eg). The reason for this is simple: Himi couples 
only states which differ by ±1 in Sz, since it is a one-particle operator (allowing for only one spin flip at 
once). Thus, a coupling of the two ground states occurs for a A^'e = 8 system first in the fourth order of 
perturbation theory {N^ = 8 implies 5*2 = 4 for fully polarized system). This interpretation fully agrees 
with the finding AE cx (-E/ji//)'*, Fig. 51b. We can therefore expect that, for an inhomogeneity of constant 
strength, the level splitting will vanish exponentially at ^ oo as long as Ejmi is much smaller than the 
gap at B sa Be- 

Another view at the crossing for Eimi 7^ is presented in Fig. 52. If we focus on the ground state and 
sweep B through Be, we may observe how {Sz) (or (5)) of the ground state smoothly passes from to 
iVe/2 = 4. The transition observed in this way (i.e. (5*2) w 2 = N^/A) coincides with the transition 
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(a) Weak MI. {Emi = 0.002) 



(b) Intermediate MI. (Emi = 0.01) 



(c) Strong MI. {Emi = 0.02) 



Fig. 53: Stronger MPs bring another ground state into play (Sz = 1) and the transition from the singlet to the polarized 
ground state becomes more gradual. 
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Fig. 54: Strong perpendicular mag- 
netic inhomogeneity, just as in Fig. 
53c, combined with inplane in- 
homogeneity (IMI). (Sz) of the 
ground state. 



observed in the spectrum (the 'anticrossing region'). Fig. 52a,b. The larger Eimi, the smoother the 
transition and the broader the range of B in which the transition occurs. Fig. 52c. 

So as to conclude: most importantly, an inplane magnetic inhomogeneity (IMI) transforms the ground state 
transition into an anticrossing. This effect should fade away for larger systems [N^ S> 1). We also remark, 
that the IMI shifts the transition point Be to lower fields. Fig. 52c, but this effect seems to be rather small 
for inhomogeneity strength not exceeding the incompressibility gap. 



3.1.3 Strong inhomogeneities 

In the following we suppress the inplane inhomogeneities again and let us study stronger perpendicular 
inhomogeneities of the form (58). 

If the strength of the 'rectangular wave' impurity becomes comparable to the gap at i? « Be, Emi ~ Eg, 

the situation at the 'singlet-to-polarized' transition changes dramatically. The excitation gap closes and 

many states of different spin polarizations crowd around the ground state. Even at zero temperature and 

in spite of lack of anticrossings of states with different Sz (i-e. Sz is a good quantum number again), the 

transition becomes more gradual, when measured by Sz of the ground state. Fig. 53c. 

Primarily, this is owing to the S = 1 state which profits best from the inhomogeneity. Keeping in mind 

its value of /c'' = (1.07, 0)^g"^, this state seems to be a spin density wave in x-direction pinned by the 

inhomogeneity potential. It is also important that states with other spins are very near to it. 

This r ~ transition can be again smoothened by an inplane inhomogeneity, as shown in Fig. 54. Here, 

the transition 5*2 = ^ 1 becomes much more gradual than the transition 5*2 = 1 ^ 4. Reason for this is 
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n| (x)/n{x). 

Fig. 55: Destruction of the singlet state by very strong magnetic inhomogeneities: the system splits into two domains, 
one with spin up, another with spin down and the electrons avoid the 'interface region' (minima in the density). 

again that the inplane inhomogeneity couples directly only states with AS*. = ±1. Other quantities than 
just Sz (e.g. polarization) are shown in Fig. 56. 

A strong magnetic inhomogeneity has also another quite pronounced feature: the singlet-polarized ground 
state transition Be shifts to higher magnetic fields. Fig. 53. This effect is considerably stronger than the 
shift to lower fields in case of the inplane inhomogeneity (Fig. 52). Origin of this shift to higher B is the 
decreasing energy of the singlet ground state. Fig. 53 or Fig. 55d. 

Let us look at this issue more closely. Increasing Emi, there is no apparent transition (crossing) in the 
ground state of the Sz = sector (not shown). The total spin of the ground state increases smoothly from 
zero and saturates around S* w 1.6 for Emi ~ 0.02, Fig. 55d. Beyond this point, the label 'singlet ground 
state' becomes inappropriate. At such values of Emi, the polarization achieves the maximum variation 
between zero and one. Fig. 55a. The eight electrons, four with spin up, four with spin down, split into 
two nearly independent groups: the spin up (down) electrons gather in the region where gi{x) is positive 
(negative), see (58). Such a state where e.g. no spin up electrons occur in the 'wrong region' (Fig. 55c, 
Emi = 0.02) is no longer even remotely related to the homogeneous incompressible state, even though 
it has Sz = 0. Rather, we could interpret it as two ly = 1/3 systems living next to each other: one with 
spin up, another with spin down. The strong spatial variation of density in this system. Fig. 55b, indicates 
that electrons try to avoid the 'interface region'; an alternative point of view is to compare the 'spin-down 
domain region' (seen in the polarization. Fig. 55a) with the density of spin down electrons. Fig. 55c. 
However, we must always be aware that we investigate only a finite system which is too small to observe 
the 'inside' of a domain where we expect the density to be constant. In a sufficiently large system, the 
maximum in Fig. 55c should spread into a plateau. Therefore, also conclusions about the interface region 
must be interpreted with caution. 



3.1.4 Quantities to observe 

Polarization is the most natural quantity to study when looking for domains of polarized and singlet states. 
Nevertheless, it could be useful to search for other observables as they might bring some more information 
on what is happening in the states. 

Here, we suggest to study the local expectation values (or densities) of otherwise 'global' operators such 
as Sz or S'^. These are defined by 



ATe 

5:,,^(r) = S'^,2(g)n(r), S^{r) ^ S^ ® n{r) , where n(r) = ^ 5(r - r,) 
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Fig. 56: Response to strong inhomogeneity of the form in (59). Emi = 0.02, Eimi — 0.005. 



and they should be plotted in the form Sz{r)/n{r). Their meaning is the following: Imagine an n-electron 
state which is an equal-weight superposition of two states: one localized in the region < x < a/2 which 
is fully spin polarized (Si ~ n/2) and another localized in a/2 < x < a which is a spin singlet. This state 
is 5*2 = n/4, yet its Sz{x)/n{x) is equal to ri/2 or in the two respective regions. 

What these quantities reveal is demonstrated in the case of a strong magnetic inhomogeneity, both per- 
pendicular and inplane. Fig. 56. The low-energy part of the spectrum does not change considerably. Fig. 
56a, even though the singlet state is separated almost completely into a spin-up and a spin-down domain 
by the inhomogeneity, seen in the polarization. Fig. 56b. Note that again the states near the transition have 
smaller variations in the polarization than the 'inhomogeneous singlet-state'. 

The lower two plots of Fig. 56b show Sz{x)/n{x) and Sx{x)/n{x). Obviously, 5*^ stays quite constant 
with X, at least on the scale ranging from Sz ~ Q (singlet) to 6*2 = 4 (fully polarized). Albeit polarization 
(or relative density of spin down electrons) varies strongly, Sz {x) remains nearly constant. This indicates 
that the state does not really separate into domains of locally different Sz- Observation of the quantity 
5^ (a;) (not shown) points in the same direction. 

Local expectation values of Sx indicate that states near the transition are more susceptible to inplane 
inhomogeneities. At any B, this response is much stronger than for perpendicular inhomogeneities. The 
following picture explains this behaviour: if we imagine a 'classical' spin vector pointing in z-direction 
and accept that it fluctuates by a small angle A(y9, then 5*2 cx cos A(/3 w 1 whereas Sx oc sin A(/3 w A(y9. 



3.1.5 Different geometries of the inhomogeneity 



Disregarding entropy, it is unlikely that a domain state will be the ground state in a homogeneous system. 
If it is an excitation we can hope to encourage it energetically by including a suitable inhomogeneity like 
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Fig. 57: Response of the singlet to various weak magnetic inhomogeneities. The inhomogeneity is similar to the one 
in (58) but it divides the system into two areas (stripes) of various ratios (6 : 6 through 2 : 10) or it consists of two 
stripes with size 3:3:3:3. 



H]\ii in (58). However, we do not a priori know what 'suitable' means. So far, we divided the system into 
two equal parts by Hmi- 

How the singlet state responds to inhomogeneities of different forms is shown in Fig. 57. Different lines 
correspond to the 'rectangular wave' inhomogeneities with different ratios of the 'plus' and 'minus' parts. 
All these inhomogeneities are thus a single stripe of various width (per elementary cell) parallel to y . 
Also, response to Hmi consisting of two stripes is shown (i.e. 'rectangular wave' with half period). 
Responses are basically very similar to each other and it seems by having focused on Hmi of the form of 
(58) we did not choose a particularly clumsy one. One particularly interesting information which can be 
extracted from Figure 57 is that the polarization response is always at least an order of magnitude larger 
than in the density, 10% against 0.3% in the present case. This confirms the conclusion of Subsection 
2.3.2: even though singlet incompressible states try to maintain constant density, they can be fairly easily 
polarized. 

A good example of the influence of the form of the inhomogeneity are the half-polarized states. Fig. 58. 
The lowest level in the S = 2 sector is six-fold degenerate (factor of three from the centre-of-mass and 
factor of two from the relative part). The two states are mirror images of each other with respect to the 
diagonal of the elementary cell. We split them into two groups J = 2, 6, 10 and J = 0, 4, 8 (within each 
group the states differ only by the center-of-mass part) and subject each group to one-stripe and two-stripe 
inhomogeneities. Fig. 57c. 

One group (J = 2, 6, 10) responds strongly to the two-stripe H^jj and is left almost unchanged by the 
one-stripe Hmi, upper row in Fig. 58b. Nearly the opposite is true for the other group. It gives a clear 
picture of the structure of these states. They are spin density waves with two periods in one direction 
and one period in another direction. This is in full agreement with the spin-spin correlation functions, 
not shown here. The conclusion is also underlined by the markedly lowered energy of the J ~ 2 state 
when it is addressed by the two-stripe inhomogeneity. Fig. 58a. This is a practical demonstration of one 
spin-density-wave state selected by an impurity from a degenerate manifold. 



adp header will be provided by the publisher 



83 



Ground state and gap 



Ml: one stripe (along x) 



-4.180 



-4.185 




Ml: two stripes (along x) 



-4.210 



hmg one- two- stripes 

(a) Energies. 



0.2 0.4 0.6 0.8 1 

x/a 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 



A 



A 



WW' 



0.2 0.4 0.6 0.8 1 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 



0.2 0.4 0.6 0.8 1 

x/a 



(b) Polarizations. 



Fig. 58: The half-polarized states (S = 2) and their response to a magnetic inhomogeneity of the form of one or two 
stripes. 



3.1.6 Transition at nonzero temperature 

Regardless of how intensely we try to help an eventual domain state to become the ground state, it may 
still be, that it is hidden among the excitations. Therefore we may try to take the excited states into account 
by means of thermal averaging. 

The strong impurity mode (Emi = 0.02) was chosen for this study. Three-fold degeneracy in center- 
of-mass of the incompressible ground states is lifted. The level splitting is however still smaller than the 
incompressibility gap, compare the black and grey points in the upper plot of Fig. 56a around B ^ 10 T. 
Various temperatures were considered: kT <C Eg means than we do not average even over all states of 
the originally degenerate triple. Knowing that Eg means the gap energy at i? or S oo, the other 
temperatures shown in Fig. 59 are self-explaining. 

Judging by polarization ni{x)/n{x), the state at the transition approaches a situation which we could call 
'domain'. In the middle (x/a ~ 0.5), the polarization drops to zero and only spin up electrons are present. 
In the other region (x/a w = 1), polarization is about 0.5, meaning that the number of spin up as spin 
down electrons in this area is the same. 

We should note though that an inhomogeneity which is strong enough to produce such nice 'domains' is 
also strong enough to change the originally incompressible singlet state completely. Fig. 59a. In other 
words, the response of the system at the transition is still weaker than the response of the singlet state. This 
manifests that spontaneous build-up of domains is not very likely within this model. 
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Fig. 59: Density and polarization at different temperatures: filling factor 2/3 with magnetic inhomogeneity (corre- 
sponding to Fig. 53c). Left to right: before, at and after transition, i.e. B —> 0, B Be and B oo. 



In the following Section we will try to suggest slightly different models which may put us on the trace of 
states which exhibit nontrivial behaviour at the transition between the incompressible singlet and polarized 
ground states. 



3.2 Systems with short range interaction 

As far as the transition between singlet and polarized ground state is concerned, the most obvious feature of 
the 1/ = 2/3 Coulomb-interacting systems is the energy 'gap' which separates the two degenerate ground 
states from excited ones even at the very crossing. Figs. 60a, 48. In the previous section we demonstrated 
that this picture may change when fairly strong magnetic inhomogeneities are applied. Fig. 53. We can 
cause a similar drastic change by replacing the Coulomb by the short-range interaction. Fig. 60b. 
Let us first concentrate on the calculated spectrum of the homogeneous system with short-range interaction. 
Fig. 60b. Again, we observe a gapped ground state with maximum spin and zero spin in the limit of 
B oo and B ^ 0, respectively. In between, states with different spins become the absolute ground 
states. Aforemost, it is the half-polarized state {S = 2), although states with other spins (S = 1 and 3) are 
not very far. Alternatively, this can be expressed by the S-dependence of the spin of the ground state. Fig. 
61a. 

The half-polarized states have been extensively discussed in Sec. 2.1 where they were studied as 'zero- 
temperature candidates' for the ground state in homogeneous systems. However, since inhomogeneities 
couple the ground state to the excited states, the properties of the lowest-lying state will not be determined 
solely by the those of the ground state. 

Spectral properties of the short-range interacting (SRI) system subjected to a 'perpendicular' magnetic in- 
homogeneity (58) are summarized in Fig. 61. In a similar fashion as for the Coulomb-interacting systems, 
states with other spins become the absolute ground state in some range of the magnetic field. Fig. 53. This 
is also manifested in the expectation value of spin (or 5"^) of the system even at nonzero temperatures. 
Most significant is still the 'half-polarized' plateau Sz{B) « 2, Fig. 61d. 

After this introduction let us look at the inhomogeneous states themselves. Their properties are highlighted 
especially in comparison to the Coulomb interacting (CI) states. When subjected to a 'rectangular cosine' 
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Fig. 60: Spectrum of a homogeneous system with Zeeman splitting (8 electrons, v = 2/3). 
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Fig. 61: Spectrum and the expectation value of the spin in the ground state in short-range interacting systems (eight 
electrons). 



magnetic impurity, those CI states showed a smooth monotonous transition from the singlet to the polarized 
state. The singlet was most strongly affected by the inhomogeneity, the polarized state was not affected 
at all, it was frozen by its symmetry. The polarized state has all spins up, Sz ~ Since magnetic 

inhomogeneity of the form in (58) preserves Sz, it does not couple the polarized state with any states which 
contain spin down electrons, since such a state must have Sz < Ne/2. The transition state was just in the 
middle. This is the finding both at T = 0, Fig. 50b, and at temperature low enough to average only over 
the three degenerate states of the homogeneous system. Fig. 62a. 

The SRI systems give a different view. The response to the inhomogeneity is slightly stronger at the 
transition than in the singlet state. Fig. 62a, inset. This is not very surprising given that there are quite 
many states near the ground state in the transition region. At slightly higher temperature where we average 
over about 10 states in the singlet and polarized limit, the distinction between Coulomb and short-range 
interacting systems weakens. Fig. 62b. 

Regarding Figure 62, it should be stressed once again, that the 'transition states' for the Coulomb and the 
short-range interaction have completely different character. In the former case, this state is basically a 
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(a) 'Low' temperature. (b) Slightly higher temperature. 

Fig. 62: Short-range interacting system with strong inhomogeneity, (58), iJj\/7 — 0.02: polarization in the singlet 
sector, around the transition and in the polarized sector, cf. spectrum in Fig. 61c. Insets show how polarizations 
fluctuate around anticipated mean values (0.5, 0.25 and for the singlet, transition and polarized state). 

superposition of the singlet and the polarized states, whereas it is a half -polarized state [Sz = for 
the SRI. 

It seems we are on the track of the domain build-up here. In an ideal case, we would expect negligibly af- 
fected singlet and polarized states while the polarization of the transition state varies between (polarized 
domain) and 0.5 (singlet domain). In real systems, we are still very far from such behaviour as the dif- 
ference between polarizations of the singlet and transition state is quite small. Nevertheless, the direction 
seems correct, in contrast to the Coulomb interacting systems. We may therefore conclude: 
(i) If nontrivial effects at the transition are expected, there must be more states involved than just the 
singlet and polarized ground states; (ii) it is likely that the half-polarized states play a major role; (iii) at 
low temperatures inhomogeneous states as in Fig. 62b can be observed simultaneously with a plateau in 
5,(5), Fig. 61d. 

The last point is a consequence of the fact that not only the ground state but also the lowest excited states 
have 5, = 2 in a part of the transition region. Fig. 60b. 

3.2. 1 Comments on the form of the short-range interaction 

For a short-range interaction, the form described in Subsect. 1.3.5, Fig. 2d was chosen. The basic idea 
there was to keep the pseudopotentials Vq and Vi at their normal values while setting the others to zero. 
Vo and Vi give the energy of two interacting particles in the state with angular momentum and 1. These 
are the states with smallest and second smallest interparticle separation possible and only the latter one is 
accessible if the particles have the same spin. 

As far as incompressible liquid states are concerned, not much happens during such 'pseudopotential en- 
gineering'. The best measure for this are directly the density-density correlation functions. Fig. 18. The 
good match between correlation functions of Coulomb- and short-range-interacting states agrees with the 
common claim that their energy is determined mostly by effects occurring at short distances. Also excita- 
tion energies remain essentially unchanged as long as 'zero momentum' states are considered (as opposed 
to charge or spin density waves). 

What strongly changes is the energy difference between the polarized and singlet state: it is 0.0632 for 
Coulomb and 0.3693 for short-range interaction in an eight-electron system, with zero Zeeman energy. 
This happens because the average Coulomb potentials felt by electrons in a singlet state and in a spin 
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polarized state differ. Roughly, we take the average over set {Vq, Vi,0,0, . . .} in the former case (all 
m's allowed) and over {Vi,0,0, . . .} in the latter case (only odd 7ti's allowed). This is a fundamental 
problem. The requirement of equal averages is not compatible with preserving the ratio of Vo and Vi as 
in a Coulomb interacting system. We would have to use higher Haldane potentials to achieve this, losing 
thereby the simplicity of the definition of short-range interaction. Therefore, with short-range interactions, 
we must be very cautious whenever we compare absolute energies of states with different spins (and thus 
parity properties of the wavefunction). Namely, position of the singlet-polarized transition depends directly 
on the energy difference of the singlet and polarized ground state. Sec. 3. 

This difficulties apply to spectra in this subsection. Fig. 61, 60b. Fortunatelly, the polarizations in Fig. 62 
do not suffer from this, provided that half-polarized states indeed become the absolute ground state some- 
where around the transition. 



3.3 Systems with an oblong elementary cell 

So far in this Chapter, we have only considered square elementary cells a by a so far If we somehow 
e.g. by means of a magnetic inhomogeneity, manage to split such a system into two domains of the same 
size, this will be a/2 by a, cf. Fig. 57c. Consequently, the spin singlet and spin polarized states which we 
expect to appear in these domains would necessarily have to be deformed as in a cell of aspect ratio 1 : 2. 
In principle, this could even suppress the formation of such domains or at least shift them to higher excited 
states. The energy of any of the two incompressible ground states depends on aspect ratio (the stronger 
the smaller the system is), Subsect. 2.4.1. There is no reason to expect that the energy of a domain wall 
between two such states is constant. In the following Section we will investigate systems in a rectangular 
cell with aspect ratio 2 : 1 which have the possibility of splitting into two square domains. All results in 
this Section refer to Coulomb interacting systems. 



3.3.1 Overview of the transition: which states play a part 

Going from square elementary cell to aspect ratio 1 : 2, the overall view of the transition changes. The 
crossing between singlet and polarized incompressible states is no longer well separated from excited 
states. Fig. 63. 

Similarly, as for short-range interaction, states with different spin appear near the transition: most promi- 
nently 5 = 1 and 5 = 2. Again, Figs. 53, 61, these states are promoted by the perpendicular magnetic 
inhomogeneities in the form of a 'rectangular wave'. Fig. 64. A consequence of this is a gradual change 
of the spin in the ground state as we sweep magnetic field (or simply increase Zeeman energy). Here, we 
should point out the difference between the present case and the Coulomb interacting system in a square 
elementary cell. Fig. 53. For an oblique elementary cell, (i) the 5 = 1 state becomes the absolute ground 
state near the transition even in homogeneous systems, (ii) A much weaker inhomogeneity is needed to 
make the 5 = 2 state the absolute ground state in some range of the magnetic field. Fig. 64c shows that 
Emi = 0.004 is sufficient for this to happen in a 2 : 1 system, while Emi = 0.02 is not strong enough for 
a square elementary cell. Fig. 53c. 

By changing the elementary cell geometry we support possible domain states, but it is adequate to ask how 
much the incompressible singlet and polarized states are affected by this procedure. The inner structure of 
these states under elementary cell variations was addressed in Subsect. 2.4.1 and we saw indications that 
the states are liquid like (and very similar to the original states from square elementary cell) even at aspect 
ratio 1 : 2. However, overlaps between the square-cell and deformed states are noticeably below unity and 
hence their behaviour is not representative if we are interested in infinite homogeneous systems. Recall, 
that the square-cell polarized state is extremely close to the Laughlin state (overlaps w 99%). 
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Fig. 63: Spectra of homogeneous eight electron systems (with Zeeman splitting) for square and oblong elementary 
cell. 
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Table 3: Incompressible ground states (polarized and singlet) in an eight-electron system. Overlaps between states in a 
square elementary cell, oblong elementary cell and oblong elementary cell with intermediate magnetic inhomogeneity. 
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Fig. 64: Spectra and Sz of the ground state in a system with oblique rectangular elementary cell (aspect ratio 2 : 1, 
eight electrons). Magnetic inhomogeneities (58) of different strengths are considered. 
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Fig. 65: 'Domains' imprinted by a magnetic inhomogeneity of type 'rectangular wave' into a system with oblong 
rectangular elementary cell (aspect ratio 2:1). The strength of the inhomogeneity is about 20 % of the gap in the limit 
B — > (in particular Emi — 0.004). Plotted quantities are averaged over the three states which were degenerate in 
the homogeneous system (in the center-of-mass part). 



3.3.2 States at the transition 

The following paragraph deals with he central result of the investigations on systems with aspect ratio 
2 : 1. The low-energy states near the transition (5 = 1 and S* = 2 in Fig. 64) respond very strongly to 
a 'rectangular wave' magnetic inhomogeneity. Fig. 65 (the middle two lines). Already for intermediate 
strength of the inhomogeneity like 15% of the singlet incompressibility gap in a square cell, polarization 
varies between « 0.5 and « 0.05, Fig. 65a (values of 0.5 and would mean a state with Sz ^ and 
Sz = Ne/2, respectively). Equivalently, Fig. 65b shows that (a) the density of spin down electrons drops 
below 25% of its average value in the spin polarized region and (b) spin up and spin down densities are 
balanced up to 10% variations in the 'spin singlet region'. At the same time, variations of the total density 
remain small (less than 5%), but there is a clear deficit of electrons in the 'polarized region'. Fig. 65c. 
In order to check that the inhomogeneity is not too strong ('destructive') compared to the Coulomb inter- 
action responsible for the formation of the incompressible ground states (far away from Be), we should 
observe the incompressible S = and 5* — Ne/2 states. Fig. 65. For both of them, responses are much 
weaker than for the transition states. 

Let us now concentrate exclusively to the half-polarized states and try to analyze their nature. Observe first 
the homogeneous system near the transition. Fig. 66 and focus on the half -polarized sector {Sz = N^/A = 
2) with one particular value of J, Subsect. 1.5.2. The low lying states show pronounced spin structures 
and, moreover, several distinct types of spin structures appear in the low energy part of the spectrum. 
This is heralded by different values of k '' which are (0, 0), (±1, 0) and (2, 0) for the lowest three states 
(st01,st02H-st03, st04, the middle pair is degenerate) and the different spin structures can be seen best 
in the density-density correlation of minority spin, gii{r). Fig. 67. Half-polarized states contain six spins 
up and two spins down here, which we choose to call majority and minority spins respectively. 
The lowest state looks isotropic as far as the rectangular elementary cell allows, the other two (st02H-st03 
and st04) are different kinds of spin density waves in the 'long direction' {x). Keeping in mind that these 
states are energetically close to each other as compared to incompressibility gaps at = 2/3 in a square 
elementary cell, for example, we can indeed expect strongly modulated polarization in response to suitable 
not very strong inhomogeneities. Polarizations in Fig. 65 were a good demonstration of this prediction. 
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Fig. 66: Low lying states near the transition (in absence of inhomogeneity, eight electrons). The J = sector is where 
(a) both incompressible ground states (singlet and polarized) occur and (b) the lowest half-polarized state occurs. 
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Fig. 67: Density-density correlations (gn, i.e. minority spin) in the lowest half-polarized states. Coulomb interaction, 
a homogeneous system, aspect ratio 2:1. 



This leads us to the question what types of spin structures can be imprinted into these states. Are they com- 
pletely 'soft' or are some particular structures preferred? An answer is given by polarizations in response 
to various types of inhomogeneities. Fig. 68. Briefly summarized: a variety of spin structures is possible 
but 'periodic' structures are preferred. Among the 'periodic' structures, the largest period available is pre- 
ferred. This means one stripe or just the 'domains' as in Fig. 65. By periodic we mean commensurate with 
the elementary cell period, for instance a 'rectangular wave' in contrast to a delta peak since otherwise, any 
structure is periodic in our system due to periodic boundary conditions. 

Looking only at polarizations. Fig. 68a, responses to all types of inhomogeneities considered here seem to 
be the same (in strength) within a factor of two. However, a closer look reveals some differences between 
those which are 'periodic' and the others. Fig. 68b. The one-stripe and two-stripe inhomogeneities mix 
mostly only the lowest four states: (sum of squares of) projections of the inhomogeneous state to states 
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same as in Fig. 57. 



(b) Projections of the inhomogeneous ground states to the lowest four homoge- 
neous states (in Fig. 67). 



Fig. 68: Half-polarized state and different forms of inhomogeneity. Strength of the inhomogeneity is the same in all 
cases, Eaii = 0.004. 



st01-st04 give in these cases > 90%. It seems that a one-stripe structure, or domain state in Fig. 65, 
stems from the /c'' = (±1, 0) states (st02H-st03) and the two-stripe structure comes from the = {2, 0) 
state (st04). In both cases, however, projections to the lowest state (stOl) remain high. 
A different situation occurs for 'non-periodic' structures like a delta peak. Inhomogeneous states are then 
'constructed' in the main from states which were originally energetically higher in a homogeneous system. 
Such states (e.g. with a delta peak in the polarization) only have a strong projection to the k '' ~ (2, 0) state 
(st04), but still more than 50% of weight comes from higher states. Fig. 68b. 

This scheme, 'periodic -welcome, others-less welcome', is confirmed also in terms of energy. While the 
'periodic' states (one- and two-stripes) profit energetically from the inhomogeneity, the delta-peak state is 
shifted to higher energy. Fig. 68b. 

Finally, the following conclusion about the v = 2/3 system near the transition seems to be possible. The 
softening against magnetic inhomogeneities of different forms, as observed in Fig. 65a, stems not only 
from the spectral properties of the system (small level spacing. Fig. 66) but also from the fact that more 
different (inner) spin structures occur among the low lying states. States belonging to a single value of S 
(e.g. S = iVe/4) are capable of generating a response as shown in Fig. 65a. 



3.3.3 What is inside the domains? 

We will now only consider the 'one-stripe' inhomogeneity, in sense of Fig. 57, which has lead us to the 
states with polarization varying almost between zero and one half. Fig. 65. In other words, we could 
distinguish between two domains of about equal areas in that state: one, with only spin up electrons 
and another with as many spin up as spin down electrons, whereby the total density is spatially nearly 
constant. Now we are interested in the inner structure of these domains. One of the aims of this thesis 
was to find side-by-side domains comprising of the incompressible singlet and incompressible polarized 
states. Unfortunately, the results presented in this Subsection cannot conclusively answer whether the 
states discussed in the previous Subsection are of this type or not. Also, it would be surprising if they 
could in the view of the small systems (eight electrons) we study. One of the reasons why studies of finite 
systems on a torus or on a sphere were so successful was that these models contain no edges. On contrary, 
there are 'edges' in the state with 'domains': the domain walls. Nevertheless, these results provide at least 
some basis for comparing the inside of the domains to the incompressible states and, in particular, highlight 
some differences between these two. 
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As a probing tool we chose density-density correlation functions. As we are dealing with inhomogeneous 
states, we must use g{r, ro) oc {5{ri — r)6{r2 — ro)) rather than g{r) cx {S{ri — r2 — r)). The former 
quantity is the conditional probability to find an electron at r given there is an electron at ro and we will 
separately address the cases when both electrons have spin up or when they both have spin down. By 
convention, majority electrons are spin up (expected to be present in both domains) and minority electrons 
are spin down (absent in the fully polarized domain). 

Roughly, we can say that the eight electrons are organized in four vertical stripes: two in the polarized and 
two in the unpolarized domain. For instance, if we catch a majority spin electron in the left stripe in the 
polarized domain, we will see another quite sharply localized (majority spin) electron in the same stripe 
and two delocalized electrons in the other stripe of the polarized domain. Fig. 69c. In the unpolarized 
domain, we will see the two majority electrons distributed nearly equally among the two stripes. 
Similarly, if we pin a majority spin electron in one stripe of the unpolarized domain. Fig. 69a, we find 
another (majority spin) electron in the same stripe. Four electrons in the polarized domain are distributed 
homogeneously to the two stripes. We will see almost the same picture with minority spin electrons, if 
we catch a minority spin electron at the same place. Naturally, we will see almost nothing in the polarized 
domain. Fig. 69b. 

Summary. In eight electron systems, the domain states comprise of four vertical stripes (i.e. parallel to 
the short side of elementary cell), each occupied by two electrons. In the polarized domain, each stripe 
contains two electrons separated by 6/2, and the two stripes can 'freely slide' besides each other Stripes of 
the unpolarized domain are preferentially occupied by electrons of the same spin and both spins (majority 
and minority) seem to be equivalent: schematically ("ft |L(ii |_r + (ii |l(TT I-R- The domains seem to be 
rather independent. For instance, regardless of where, within the unpolarized domain, we pin the majority 
spin electron, the density of electrons seen in the polarized domain does not change much. 
It should be emphasised that although the stripe structure is well pronounced in conditional probabilities, 
the density varies only weakly. Fig. 65c. But, even so, it contains indications of the four stripes. This 
structure suggests that the interior of any of the domains is rather anisotropic and this is quite distinct 
from the liquid states at u = 2/3 (polarized and singlet. Fig. 13) where at least the first maximum in 
g{r) occurs for all r with \r\ = ri (Subsect. 2.1.1) and not only in the x- or y-direction. The study of 
finite size effects comparing the averaged and non-averaged correlation functions suggests that for liquid 
states, the anisotropy of non-averaged correlation functions should be much smaller than what we observe 
in Fig. 69. On the other hand, the results shown in Fig. 69 refer to a state which is inhomogeneous and 
strongly influenced by the aspect ratio being far from unity. A more thorough study of the non-averaged 
correlation functions in systems of various aspect ratios and comparison to systems of different sizes is 
therefore necessary to allow more definite conclusions. 



3.3.4 Comment on homogeneous half-polarized states 

It should be mentioned that the half-polarized state we study here is not the same as the half-polarized 
states discussed in Sect. 2.2. 

The lowest state here, in an elongated elementary cell, has k '' = (0, 0), whereas the half-polarized ground 
state in a square cell has k "" ~ (2, 2), Sect. 2.2. These are the two inequivalent points of the highest 
symmetry. 

In the present system (aspect ratio 2 : 1) all low lying half -polarized states belong to the ky = sector. 
States with other values of ky lie well above the four discussed states st01-st04, the lowest of these other 
states has an energy of —4.240, cf. the spectrum on the left in Fig. 66. Being interested in the states low in 
energy we may therefore stay restricted to the sector ky = 0. This is a pleasant fact since there is no need 
to consider inhomogeneities of very low symmetry (implying handling larger Hilbert spaces). 
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Fig. 69: Half-polarized state with intermediate magnetic inhomogeneity ('rectangular wave', Emi 
averaged density-density correlation functions. 



0.004), non- 



Also note that a spatially fixed inhomogeneity couples the relative and center-of-mass (CM) coordinates. 
At filling ly = 2/3 there are three possible CM states on a torus, which however remain mutually decoupled 
owing to the high symmetry of the inhomogeneity. The chosen sector J = corresponds to a combination 
of ky = sector and the one CM state, which leads to the lowest energy. Differences to other CM states 
are, however, not too large. 



3.4 Summary of studies on the inhomogeneous systems 

Perhaps the most important conclusion of this Chapter is that the two incompressible ground states at 
v = 2/3, the polarized and the singlet one, alone are not enough to create a state with 'domains', i.e. 
regions of polarization zero, corresponding to Sz = 0, existing side-by-side with regions of polarization 
one. We have demonstrated this in Subsection 3.1. When a 'domain-inducing' magnetic inhomogeneity 
is applied, the singlet ground state is more strongly affected than states near the transition. This claim 
remained true for different types of magnetic inhomogeneities, for different quantities used to detect the 
domains (apart from the polarization, also for Sz{x), Sx{x), etc.) and also for non-zero temperature. 

Different conclusions apply when more than just the polarized and the singlet ground states are present 
in the low-energy sector We have demonstrated, that near the transition, the gap could actually close in 
several different situations. In the present study, this happens for very strong magnetic inhomogeneities 
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(Subsect. 3.1.3), for systems with an elongated elementary cell (Sect. 3.3) and for short-range interacting 
systems (Sect. 3.2). The states which closed the gap always belong to an intermediate value of spin, 
most prominent are those with S = 1 and 5 = 2 and since we considered only eight-electron systems, 
the latter value of spin corresponds to the half-polarized sector S = Ne/4:. These states are considerably 
softer against magnetic inhomogeneities than the incompressible singlet and polarized states. On one hand, 
this fact follows from a small level spacing in the low energy sector when the gap closes. However, the 
magnetic inhomogeneities were also found to have large (~ 0.1) matrix elements between most of the low 
lying states. 

The states with the strongest tendency to form domains i.e. the 'softest' states, were found in systems 
with Coulomb interaction and an elongated elementary cell. Near the transition, even a moderately strong 
magnetic inhomogeneity (weaker than the incompressibility gap) was enough to make the polarization 
approach the values corresponding to the singlet and polarized states inside the domains. For these domain 
states, we have investigated the 'inside' of the domains by means of non-averaged correlation functions, 
Subsect. 3.3.3. We could not yet confirm that the domains comprise of an incompressible liquid but this 
system defintely deserves a more detailed study. Especially in this case, a comparison with larger systems 
would be very helpful. 



4 Conclusions 

Fractional quantum Hall systems at filling factors = 2/3 and 2/5 have been studied numerically by 
means of exact diagonalization techniques on a torus. In both systems, the existence of two different 
ground states is well established: one is fully spin polarized, another is a spin singlet and they are both 
strongly correlated. All four states can be visualised as composite fermion systems at integer filling factor 
(i^CF = 2). A transition between these two ground states can be induced by changing the Zeeman energy 
while keeping the filling factor constant. Chapter 3. 

At the beginning of Chapter 2, we investigated the polarized and the singlet incompressible ground states 
in terms of their density-density correlation functions. First, we highlighted the fact that - even if these 
states were exactly described by some composite fermion model - the inner structures of the ground states 
at 1/ = 2/3 and 2/5 differ strongly from the inner structure of a state comprising of two fully occupied 
Landau levels. In other words, in a composite-fermion state (e.g. lycF = 2), the correlations between the 
electrons are different than in a corresponding electronic state (i/ = 2). A more important result is, however, 
that the electronic correlations differ strongly also between the i/ = 2/3 and 2/5 states themselves. This is 
surprising, since both filling factors map to the same filling factor of composite fermions {I'cF = 2) and 
only the orientation of the effective field is different. Study of the correlation functions allowed to suggest 
a new interpretation of the singlet 1/^2/3 ground state. The electrons move along in pairs of opposite 
spins and the pairs form a state equivalent to a fully occupied lowest Landau level. This conclusion does 
nof apply to the = 2/5 singlet ground state. 

The central focus of the present work was on the low-energy states occurring near the transition between 
the singlet and the polarized ground states. Some experimental results indicate that another ground state 
distinct from the two ground states already mentioned could exist near the transition [46]. In Sections 3.4 
and 2.2 we found several arguments in favour of a half-polarized state (S = iVe/4) becoming the absolute 
ground state in a narrow range of the magnetic field. The systems available to exact diagonalization were 
however too small to allow for an unswerving prediction. Two candidates for the half-polarized ground 
state were identified. In Section 2.2 we concentrated on the 'isotropic candidate'. A study of its inner 
structure (correlation functions) combined with an investigation of the response to probing magnetic inho- 
mogeneities (Sect. 2.3) produced results resembling both the singlet and polarized incompressible ground 
state. A hypothesis that both these states coexist within the half -polarized state has been presented. 
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Calculations with elongated rectangular elementary cells (Sect. 2.4) suggested another candidate for the 
half-polarized ground state: a spin-density wave along the longer side of the elementary cell. A comparison 
between two systems of different size indicated that this state has the shortest period allowed by the finite 
size of the considered system (e.g. one third of the length of the cell for a state which contains three 
minority spins). Based on the present calculations it is not possible to decide which of the two candidates 
(if any) evolves into the ground state of an infinite system. 

At = 2/5, no obvious analogue to the half-polarized state at i/ = 2/3 was found. 

Employing magnetic inhomogeneities to enforce domains of different spin polarization near the transition 
at J/ = 2/3 (Chapter 3) we found that no signs of domain formation occur unless the energy gap closes. 
The loss of incompressibility could however still be compatible with the experimental observation of a 
plateau of polarization one half during the transition. It is enough if there are many states with S = Ne/4: 
and no (or only few) states with other values of S in the low-energy sector (Sect. 3.2). 
The 'best' candidates for domain states were found to appear in systems with an elongated rectangular 
cell. The fundamental idea here was that the elementary cell with aspect ratio 2 : 1 is divided by the 
inhomogeneity into two square parts which could be more convenient for the formation of isotropic states 
(the singlet and the polarized incompressible liquid). Examination of the domain state however showed that 
the inside of the domains does not resemble the incompressible ground states alv = 2/3. Nevertheless, a 
more detailed study is necessary here, since systems with aspect ratio far from unity can suffer more from 
finite size effects than what was demonstrated in Sect. 2.1. 

At the very end, I would like to acknowledge Daniela Pfannkuche for her support during my PhD stud- 
ies and many fruitful discussions, Benjamin Kriiger, Philip von Ende and Ondfej Certik for some of the 
calculations presented here and Matti Manninen for his hospitality in the period when this work was be- 
ing completed. Finally, I am much obliged to a group of my colleaugues who carefully read the very 
long manuscript and helped to substantially improve both its contents and language: Frank Hellmuth, Rob 
Knapik, Katrin Malessa, Christian Miiller, Michael Prouza, Arek Wojs and Jan Zemen. 
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